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MECHANISM OF ENERGY LOSS 


FLUID FRICTION 
‘By Boris A. BAKHMETEFF, 1M. “AM. Soc B. 


4 es “account of the inner processes, by which the ener energy r of fluid flow is 

lost. through friction, is presented in this paper. Previous investigations | have 7 

dealt, largely with bulk appraisals, ‘relating to over-all losses of head. _ The 


first expressed the nal energy ‘balance for motion i in viscous 
formulated. the 80-¢: -callled ‘dissipation function.’ Horace Lamb‘ broadened ithe 


Stokes reason reasoning, applying it to laminar motion in general. Move recently, an 


energy the and its final dissipation ‘molecular no 


| cide anid and that they | constitute the initial and the ultimate phase « of a a 
manifold proces 
2 The pap paper reveals these consecutive phases, explaining their dynamic sig- 7 


‘itil ificance and the nature of the losses involved in the different stages.  Thefacts 

ae disclosed are 2 especially pertinent for for gaining a more comprehensive nderstand- iS a! 

* ng of the mechanism of turbulent flow. The analysis begins with 

‘mniform established motion and i is extended to layer = 


_ Nors.—Written comments are invited for immediate publication; to insure publication the last sail ae 7 
cussion should be submitted by July 1, 1945. 


Prof., Civ. Eng., Columbia Univ., New York, N. Y. 

Associate Prof., Civ. Eng., College of the City of New York, New York, N. Y. 

-_ 3“On the Effect of Internal Friction of Fluids on the Motion of Pendulums,” by Sir G. G. Stokes, 

Bere Cambridge Philosophical Soc., Vol. [X, 1850, p. 8; also the collected ‘papers in Vol. III, p. 7. 
‘**Hydrodynamics,” by Horace Lamb, 6th Ed., Univ. Press, Cambridge, England, 1932, p. 579. es, os - 
5“Statistical Theory of Turbulence,” by G. I. Taylor, Proceedings, Royal Soc. of London, Series A, _ Sy 


§**Modern Developments in Fluid Dynamics edited by 8 
tngland, 1938, ‘Vol. p. 395. 
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LOSS. 


‘ee - The fundamental relations are catelilideed first for the simpler case of 


laminar motion, and subsequently are a are applied to to the more complex turbulent 


flow. 


The letter symbols i in this paper, defined in in the t text t where they first appear, _ 
essentially to American Standard Letter Sym mbols for Hydraulics = 
(ASA—Z10. .2—1942), prepared by a Committee of | the American Standards \e 
Association, w ith Society representation, an and | approved by the Association i in y 

= 


1. Intropuction Morton 


in psy p= = pressure intensity; / = distance ‘measured in the direction of 

. flow; y = = specific weight (weight per unit: volume) ; and h= head, measured | 
froma line. In this instance Eq. 1 measures the rate of energy per 
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Fiow 1n Smoots Pipe 

volunie | lost a by each and every element of & Cross section, in its dis- 


over a unit length of the conduit. The s same quantity (Eq. 
ser rves to determine the mean wall friction 


with that, the internal shearing | stress structure. The lat latter, in uniform 
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Values of 
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radius; 8; and y’, in two-dimensional motion, designate specifically 

_ the distance from the free surface (Fig. 2) or from the central axis as ag against Y, 

hich exemplifies a transverse coordinate i in general 
Another familiar concept is the bulk loss of power: 


dp o( B24). ft 
sec per Fe, 


a ele asa 43 

Eqs. 2 and 3, from the axis to an inner filament; 

sectional area; Py = the wetted perimeter; R = the 

eg Unit, 

Length ) 
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ENERGY 

or the total w wank done by the: flow in overcoming frictional resistances during 

= time and over a unit heal of the conduit. In Fig. 3, the discharge Q 


enters the unit length block and sub- 
hy=S1 sequently leaves it, in the interim unit 
ll time Suffering the loss W.; or, one may 
a state that within the unit length block, 
during each” second, mechanical 
yy energy at the rate of W, is withdrawn 
- from the flow and is spent to “overcome. 
resistances, by way of ultima ate dissipa- a 
‘tion into molecular heat. . Since the 

Velocities: do not: change, all of W, is 


s ithdrawn from the potential energy store carried by the flow, thereby caus- 
ing a reduction of the piezometric head head by hy = The foregoing 


 - 
Ww 


ith the subscript b designating the process of withdrawal or the ‘ ee ' 
of the energy from t the flow, and the subscript s indicating t the spending of the def 
energy onresistances.; 
The equality sk shown in Eq. 5 holds for bulk re relations only. ea In fact, con- i. bo! 
1 prt to Kq. 5, the local values of Wp and w, are generally not equal and do not for 
balance each other. Indeed, the cardinal fact which dominates the entire 
energy loss — mechanism is that the cross-sectional distribution of the | local Tes 
borrowing quantity wy and the local spending rate w, We ‘manifest altogether. a 
‘different, and in a sense opposite, outlines, and that for such reason the borrow- . - 
a = and the spending phases a are necessarily linked by an intermediate process, ( 
4 the function of which is t to transmit energy withdrawn ¢ at one part of the cross: pr 
Boo to other par arts, where the mechanism of energy exchange. calls for ‘its ca 
Tue Enercy W ITHDRAWING OR BorrowinG Function 
The distribution of the | local borrowing rs rate We is is determined by noting 
re. that the energy - withdrawn ; per unit time within an elementary tube dA (Fig. 3) 
is dis S dQ = Y S u dA, w with | u the local velocity. — Dividing by the volume a 


é dA X 1, one e obtains the unit rate “rege is, the work over a unit time interv: al, 


which | shows that the spatial distribution of of we follows the local velocity pro- 
“file, copying the latter to the scale ¥ S (curves. BA and B;Ai in Figs. 2(a) and 


+ In two- dimensional flow, for a block of unit width normal to the plane 


(W)o = ¥ 
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and the average cross-sectional unit rate will be 


Furthermore, the e aggregate energy amount ‘ ‘withdrawn” during a unit time 
from a partial « cross-sectional block, such : as that shown in Fig. 2(a) by the 


constituting the * “cumulative” borrowing rate between the transverse limits, as- 


indicated. The quantity “measured by the area OBba in in ‘Fig. will 
be termed the “cumulative discharge. the 
corresponding quantities [Waly are defined in 
Fig. 2(c) by curve with the coordinate 
be representing the area under the 
Bib; of the local ws-curve in Fig. 2(b). 
Obviously, for y’ = y’., equals as 
defined in Eq. 


‘Differ ential Form —To express the energy 

borrowing process ‘mathematically, note 

“for a dimensional element. (Fig. 4), the 


‘Tesulting friction force is dy dz — ‘In uniform 


“established m motion frictional traction is in balance with the pres result: ant 


dx ) dy which, apportioned to a volume unit, defines t the ‘equilibrium 


_ The sign of the stress gradient dr/dy y is necessarily negative, for the friction 


id resultant « opposes the motion and com- 

“gravity force, actuating the flow. The 
"stress is essentially : an internal agency. 


on the interfaces n and m of 
%+) joining filaments (Fig. 5(b)) is equal in 


size, but opposite in direction. In de- 


(a) 7 b) 


quently, one is guided by the direction of 
the traction on the face of a fluid block (Fig. _ 5(a)), terminated by the positive — 
ordinate y. As is evident, the sign of 7 will coincide appropriately with th: at 


0 Of du/dy, 
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kK. he unit time work performed by friction on the element (Fig. 4)is Fy ae dx dy u 


which is also a negative quantity to be compensated by the positive wniceti of the equa 
pressure resultant ((— dp/dx) dx dy uj. The latter, in turn, is equal to the valu 
positive unit rate at which energy i is withdrawn locally from from ‘the flow flow. Zz leis | Neit 


_ Note that, _— far, no no limitations have been imposed on n the particular nature hon 


; of the stress- -generating agencies. _ Therefore, Kq. 11 is quite general and will F- 
apply. to all cases of two-dimensional flow, whether laminar or turbulent. 
3. THE _ENeRcy ‘SPENDING PHASE gu 
Iti is expedient to first. purely viscous (laminar) pa patterns, in w whisk 
thes mechanism is relatively simple. Indeed, in accord with Stokes, 3 the spend- 


ing of energy has been identified with 
work involved viscous deformation. For 
© two- dimensional motion (Fig. 6), when the 
element is deformed into ab’e’ d, ‘during 
the time dt, the done is the "product 

‘the shearing force dx and 1 the displace. cen 

ment du dt = dy dt, the local unit time 

rate. of deformation work per unit volume is equal to 


the 

‘80 

or for laminar motion, with = axl 


_ Energy spending is lisse concentrated in the 1 region of high velocity gra- 


ee. which is principally near the walls of a conduit. In fact, at the axis, 


the local loss is nil. outline of the e local energy § spending 
curve w, (curve 0; ‘si, Fig. . 2(b)) i is thus ‘ “opposite” to that « of Wee 


analogy to 9, the cumulative spending, rate— 


a. is the aggregate work spent ¢ on resistances per per t unit time in the shaded | block ni 
de of Fig. 2(a), being represented i in Fig. 2(c) by curve 01828, i in which the coordi- eC 
‘measures the area a in ‘Fig. 2(b). For al total reach 0— pe 
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Naturally the areas O,B,A; and O,S8:Ai under curves Ws and Ws in Fig. 2(b) are 
equal. a any intermediary position, however, such as. , the cumulative 
values a and |W, |v" shown in Fig. 2(c) by and a2s2, do not coincide. 
Neither do the local values ws and w, represented i in Fig. 2(b) by arbi and ; asi. 
Indeed in the central region of Fig. 2(b), between the axis and the intersection 
point M, the borrowed energy We exceeds the rate w, at which energy is spent 


locally, shen, nearer to the wall bey yond M, the local ene energy rates, required 
for compensating the deformation work, exceed by far the comparatively small 
quantities we, available locally by borrowing. . It follows, accordingly, that, 


between the ‘ ‘borrowing” and the ‘ “spending” phase, there must necessarily be 
an intermediate ‘ “transmittance” function, by means of which the e excess of 


borrowed ¢ energy, available in the central re region o- -y' M; is conveyed tc to the wall 


zones to. cover | the requirements for spending. 
4. Taz Enercy Bauance Equation 
The local energy balance can be expressed by 


the cumulative curves W, and W. (Fig. 2(c)) terminating at the same ‘point Be. 


with we = wy — Ws designating the local “transmittance e energy rate.” 

central reach, with wy > Ws, We is positive, exemplifying the eXCess | of energy 
subject to transmittance. hh the wall region, where Ws e ceeds the local 
supply of of ‘we, We becomes negative, signifying that the deficiency in locally 


borrow ed. energy is supplied from the store accumulated by transmittance Stone 


Fig. 2(0), the (BT) )-eurve is obtained by subtracting O, from 
80 that ‘the local value for 1 y’ is sib) = = ait The resulting curve crosses the 2" 
axis (w, = = 0) at the intersection m, reaching its maximum negative 
A iT, at the solid boundary. 


In Fig. 2(c), the curve OstarAs | features the 


obtained by deducting 0.8282 from OsbsBs. the 
= asbe measures the area 0,B;bis: O.Bitia: in Fig. 2(6), indicating 
the aggregate excess energy, transmitted from the block toward the w all re- 
gion. The maximum cumulative value medion in Fig. 2(c), corresponding in 
position to the intersection point M of the local ‘We- curve and w.-curve, is the 
total excess energy, accumulated over ‘the reach and ‘made available 
through transmittance for spending i in the wall zone. * Obviously the coordi- 
nate in Fig. 2(c) represents ‘the positive excess area O i in Fig. 2(b), 
equal to the negative deficiency area MS.AiM. The decline of beyond 
point M signifies the gradual consumption in the wall region of the energy 
accumulated by transmittance. At the solid boundary the transmittance 
‘energy exhausts its functions. Accordingly, the curve goes 
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o grasp the “ transmittance” ’ concept consider, in two- dimensional flow, 
‘an element A (Fig. 7) i in its relation to the adjacent layers B and C (for clarity 
the interfaces are shown separated). In its mo. 

| tion 3 relative to layer B, element A pares 


positive traction, performing work on layer B at 


| the | unit time rate of + dz ue In turn, element 
A -Teceives work» from layer at the rate 
(@r re . dr) dx (u - + - du). The resulting positive 
performed by the element | on the adjoining 
fluid, over what it receive es, is rdru — (r + dr) 
ut du) dx which, after omitting 


higher order of 


Dividing by the volume dz dy, the unit local rate of work | transmitted to 


the adjacent fluid or the transmittance unit rate oan 


Eqs. 11 and 12a, energy balance, Eq. 14, for two- 
laminar flow, takes the basic differential fom 

_ 


Formally, by transposing his nisin Eq. 17 may be peri into: a complete 


- differential of 74 u, which furnishes the starting point in the Stokes* analysis. 
However, to give Eq. 17 the proper physical significance, the members must 
be presented i in the algebraic sequence as written. Note that, for established 
uniform motion, with necessarily negative (see Eq. 10), the quantity 
ier dy = is positive. . The quantity w, is also positive, for in T a both 


- members are of the same sign (see paragraph following Eqs. 10) . The only 


quantity that changes in sign is the local transmittance rate w, (Eq. 16). ~ 


-_* Since w; is the derivative of — (r u), the latter. must feature the ae 


transmittance rate W i In fact, for the block 0-y’ in Fig. 2(a), by Eq. | 


lw, wedy ly’ = be tu) dy’ = 


quantity u T is negative, being the product of a positive velocity and a stress, 


which in terms of Fig. 8 is is negative. Hence Tu) As ‘positive, indict ated 
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by curve in Fig. 2(c). The sign the local value Wy O on the other 
hand, conditioned by the tangent to the Weeurve, « changes from positive to 
negative at point tyr. 
As a matter of procedure, when “dealing unt Length 

with the transmittance phase, it expedient: 0 abi 
to start with the 

according to Eq. 18, is the product of “the: 

velocity curve (BA i in Fig. 2(a)) and the stress _— ee J 

line (OC in Fig. 2(a)). The maximum v: value 

of —ur is reached at tl the intersection ‘point ¥ 

Mi in Fig. 2(b) corresponding to we = Wy — Ws 

0. Point M is to be referred to 

“reversal point.” 

Cumulative Energy Balance Equation. — 

The local equation. (Eq. 17) is complemented 


for two- dimensional uniform. flow by the ‘cumulative expression 


or, in relation. to Fig. 2: 


“all the members (see Fig. 2(c)) being positive. ad or ot * 
be Again, as in the case of Eq. 11 and Eq. 9, no re were imposed on 


nature of the stress-forming agency. The “transmittance” phase for mulas, 
Eqs. 16 and 18, apply to all forms of flow, laminar or turbulent. nai 
6. MECHANICAL ANALOGY Jee 
‘The essence of the “transmittance | function” may be icieialiaa by a 


mechanical model. In Fig. 9 a sequence of wheels is mounted on a common 
frame, with adequate friction between the contacting surfaces. Because of this 
‘friction, energy derived from the “driving” agencies P,, Pe, ..., is transmitted 
¥ to where it is spent on the resistances R. Slipping between on wheels 
will that, although the friction components such as ihe = 2=P 

and ITrler C (Fig. 9(6)) are > equal, the velocities u: and Us differ, and thus only 
a part of the power |Tvs|us derived from P; (namely, |Poslu) will be trans- 


Tepresents a local mechanical ‘ ‘energy loss,” an amount dissipated directly 
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agencies (Pit to P, Pi) through wheel wheel tis 


into molecular heat. The work done cumulatively by all the preceding driving 


Of this work, however, only By u; is passed to the next w heel (7), the amount 


FP; (ui — uj) being dissipated locally in the slip betw een and (j). A As 4 
— 


result, the “cumulative” energy balance for the wheel is s expressed 

a 

resulta 


stress 7 in Fig. 8 is is equivalent, to the borrowed energy AL 


4 the local “slip” loss (Eq. 22) compares with w, in Eq. 12a; and DFA Au) in 
4 


= 
=>, 
© 
| 


Eq. 24 is the analogue of |W. ie _ Finally, the ‘ “transmittance” ri member 1; Fi 
is the both in and of the product Ir u| =| 
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Enerey “ ‘SPENDING”’ 


i Loss of energy signifies ultimate dissipation into molecular heat. Energy 
on ed from the flow changes its mechanical aspects and i is transformed into 
thermal form. . Such transformation is caused finally by viscous attrition, 
which enhances intermolecular activity. The ultimate dissipation develops _ 
thus on an ultramicroscopic ‘molecular 
In laminar patterns the Process | of dissipation 1 is assumed to coincide with 
viscous deformation (see Section 3). 3 _ Mechanical energy s spent i in oven 
turbulent flow, “spending? and “ portly. The 
term ‘ ‘spending”’ will be used to signify only the idea that energy borrowed from 
the flow has lost, irretrievably, its original flow form; and “dissipation,” under — 
all: circumstances, will designate the ultimate “degradation”: of mechanical 


ENERGY BALANCES In LAMINAR 


energy exchange process by purely procedure. 
Two-Dimensional M otion, Between Parallel Plates (Fig. 2). —wWith 


»= Uy’, as the unit width discharge | between the axis and the wall, the local 


is 


in Eqs. 7 and 8, 


By elementary algebraic procedure. (w hich is soil oo 
~ following dimensionless expressions fc for the local and cumulative unit rates that 
appear in Eqs. 17 and 20: 
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| Pipe: the Poiseuille re relation, 


with R = ~ and up = 2 U, leads to the following oe 
80a) 


In establishing the local rates, it is expedient to operate with an elementary 
cylindrical shell. For example, the w ork involved in the deformation shown 


5} 


ig. 10(a) results in the formula: 


~ dr 


The cumulative quantities es should refer to an internal cylinder of radius r and 


length. The: final inter are: 
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The cumulative cu curves reproduced for with the ev entual 


4 


Turbulent Flow 
—o— Laminar Flow 


Values of 


Since n no restrictions pe introduced to limit the expressions for the bor- 
rowing rate | Wp » (Eq. 11) and the transmittance rate w, (Eq. 16), the equa- 
tions are directly applicable to turbulent flow. z Naturally, 7 will be the local 
Reynolds stress, and a the mean steady local velocity, whereas the quantities - 


Tu) in 18 and Eq. 21 the tem- 


average unit borrowing a complex 
nature of turbulent motion, the rate Ws be a priori 
as in 1 the case of Eq. 126 for laminar flow. J On the other hand, 7m the stand- : 


SMOOTH PIPE | 
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—must « conserve the mathematical form of Eq: 12a. The: w,-quantity indicates 
4 the unit rate at which flow energy disappears under the action of agencicn 


pertinent to the turbulent process, i” 


Furthermore, ‘the incumbent stress and velocity quantities | cannot be 
Tived from analytical expression as in Section 3, and must be taken directly 
from observation. Empirical data vai vary Ww ridely, wall roughness, 


energy y balances for turbulent flow display certain 
which happen differ substantially from lan laminar outlines, and reflect the 


essence of turbulence dynamics. _ A typical example i is given in Fig. 1, based on 7 
 Géttingen data’ relating to turbulent flow in a 10 cm smooth pipe, w ith U = 876 -t 
em per sec, and — dp/dz = 480g percm® 


: a By us using in Eq. 34 the values for 7, a, and da/dy contained in the Nikuradse | 
report, one computes the consecutive local unit | rates D, De, and @, as plotted :_ 
in Fig. ‘1. The cumulative curves, obtained from the local outlines through tt. 
step-by-step integration are presented, in dimensionless form, in Fig. 
together with the comparative plottings for laminar flow. 
The . striking feature, characteristic of all turbulent patterns, is the intensive 
concentration of energy spending in a narrow zone next to the wall. For 

- laminar flow the reversal point (M) i is located at a distance 0. 293 ro from the 
wall, whereas in turbulent motion it is shifted to within 0.053 r, from the solid 
—. Furthermore, over the Teach constituting about 95% of 
the radial distance, relatively small amounts of energy ‘are spent locally, the 
major part of the locally borrowed energy wz being transmitted to the wall zone 

ak 
(ru —to) where the accumulated energy is used to offset the unusually high 
requirements for s spending. In fact, in Fig. 1, the calculated local values of 
‘w, and We close to the solid boundary outreach by far the bounds of the figure. 

Accordingly, to appraise the phenomena, one must refer to the dimensionless 
cumulative diagram, ‘Fig. 11, in which 1 the specific characteristics of the energy 
exchange i in turbulent become particularly lucid. ‘Through the central 
-Tegion, the cumulative transmittance curve e |W,|’ follows closely the borrowing 
outline In fact, over the reach 0-—M, in which more than 9 90% of 


‘the total energy Wo i is s withdrawn (borrowed) from the flow, the cumulative 


05 re) wall ‘zone, where the q quite insignificant local borrow ing accounts for 
less than 8% of W.. 7 In other words, the concentrated “spending” in ‘the wall 
zone takes effect principally at the expense of the transmitted energy. be os 


One should infer, : from the shape of the turbulent curves, that the energy 
“exchange i in the central zone and in the border zone must be actuated by en- 


tirely different Conventionally, i in fact, one may accept the 

premise 2 that the “reversal radius” ry separates the fluid body into two parts, ; 
of which the central cylindrical block (0 —rm) will be called the “transmittance 
region,” whereas the annular space near the wall (ru- M—To) will be referred al 


**Gesetzmiissigkeiten der turbulenten Strémung in Glatten Rohren,” by J. Nikuradse, 


‘Verein Deutscher Ingenioure, No. 356, 1932, pp. 20 and 31. iets 
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Turbulent entent is inane by the pre presence, in the flow, of an array 
of eddies that “sw arm” in an apparently unpredictable manner and at each 
_ instant add their rotary - components to the e respective local velocities of the 
main motion. * ~The ever-fluctuating turbulent pattern, superimposed over the 
steady flow s structure, devolves thus from the composite action shal a ‘multitude — 
of individual vorticing units, 
os The energy of a fluid in turbulent motion is cor compounded from the energy yo 


of the main flow, characterized by the energy 


29 


= value of (w’s)2/2, where 


is the temporal mean square of in ‘the three 


al statistically permanent state of ‘turbulence implies ° the maintenance of : a 


steady average eddying structure, which n means that new vorticing units must 


course of their | swarming the | flow. generation of 

eddies infers that: they are endowed with energy, W hich must ; necessarily be 
supplied from the energy store of the main flow. Indeed, the primary charac- 
teristic of the turbulence mechanism is that energy, at an apportioned rate, is 

withdrawn continuously from the main flow, and i is “invested” into turbulent 
eddies, such | “investment” involving the conversion of the borrowed flow energy < 
into. vortical eddying form. The generation of of turbulence i is thus synonymous 
with the producing of the turbulent eddies. — Turbulent e energy, manifest in. 
turbulent fluctuations, is originally flow energy, eoniornal into vortical form. 
By its very nature the conversion is an irreversible process, because energy, 
- once converted into eddies, can by no means s be restored to the flow. For such 
reason, so far as the main flow i is concerned, when energy withdrawn from the — 
latter is “spent” in the creation of turbulent eddies, it be be considered 
Actually “lost” flow energy reappears, at least partially, in the 


converted structure. The final 


in the course of swarming, the mechanical energy with which the eddies were 


initially endowed is transformed (degraded) through viscous attrition into 


According to views, the generation of turbulent eddies in estab-— 
lished “friction” flow i is largely concentrated i in a narrow zone near the solid 
boundaries of a conduit. Prandtl? most appropriately qualifies his definition 


* Proceedings, 5th International Cong. for Applied Mechanics, John & Sons, 
Ine., New York, ‘N. Y., 1939, p. 367,00 


4 
, _and by the kinetic en — 
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of such zones as “the at actual eddy producing mill,” whereas the 


passively.” Indeed, the of generating ed eddies seems to be a specific 

‘property, inherent to discontinuity ‘ ‘vortex sheets” that are » subjected to ex- 
excessive st strain. Whenever they are overstrained, such vortex sheets, or layers 
of concentrated vorticity, become unstable; and after passing through a cycle 

ever-increasing self-induced pendulations, they terminate eventually 
“curling,” or being “rolled” up, into ‘ “chains” of eddies or a sequence of inter- 
‘mittent lumps of vorticity. The process is comparatively well de- 
fined in cst case all “free” discontinuity sheets caused in the interior of a fluid 
As illustr: ated by Fig. 12 the: formation of the eddies, in 


Fre 12.- “ 


= instances, may be observed visually. bt (In these observations by L. 
"Schiller and his associates," the original vortex Sheet ws was made vi visible by the 
a injection of. dye.) In friction flow, past: solid surfaces, the reach of excessive 
_ strain is is supposedly associated with the | boundary zone near t the wall) where, 
: over an . exceedingly thin ‘ ‘sublayer,” the velocity drops abruptly tozero. /Thus— 
far, existing techniques have not permitted the direct observation of turbulence 
“generation near solid boundaries. There i is ample indirect evidence, neverthe-_ 
less, pointing to the essential similarity of the case with that of “free” separation — 


the reader is referred to the general literature in this field.)"1* 


q sheets. (For further information in this generally still obscure and mysterious — 


10 Proceedings, 5th International Cong. for Applied Mechanics, Mass., John & Sons, 
 Ine:, New York, N. 1939, p. 315. 
1**Fundamentals of Hydro- and Aeromechanics,” based on of L ‘Prandtl, by O. G. Tietjens, 
4 ‘translated by L. Rosenhead, Ist Ed., McGraw-Hill Book Co., Inc., New York and London, 1934, p. 222. ; 
a 12‘*The Formation of Vortices from a Surface of Discontinuity,” by L. Rosenhead, Proceedings, Roy al 
of London, Series A, Vol. CXXXIV, 1931-1932, pp. 170-192. 


= 


— 
im 
— 
iii 
4 
| 
— 
— | 
: 
| | 
| 
| 
4 
| 
| 
| 
j 


The pivotal phase in the sequence of events is conversion, in ne course od 
ow wide turbulent eddies are generated incessantly i in the conversion Zo! zone. — 
preconversion phases (Section 9) deal with energy of the main flow. ‘As illus- 
trated by Fig. 11, flow energy is borrowed largely in the middle regions of a 
_CTOSS section and through transmittance is concentrated near the wall, being ; 
made available for conversion. In the post-conversion phase, the remainder © 
- of the originally borrowed mechanical energy is in turbulent oddying form 
: Individual vorticing units, assumedly of the Fig. 12 type, are ‘“‘cast off” or 
a the generating zone into the adjacent flow, and as the eddies 
‘spread across the conduit. ‘on generally oblique | paths, their initial energy. is 
gradually dissipated into molecular heat. The’ possible interreaction between 
a vortical unit and the e surrounding | 
flow is illustrated by 13. On ot 
top of the eddy, the "rotary 
= adds to the local 


velocity ux, Whereas at the bottom 
subtracts from Us. The resulting 
difference in v velocity 
ina transverse Pressure differential, 
“prompting the eddy y to move 
crosswise ‘direction away from the 
wall, The e “shedding’ ” of the e eddies: 
from: the seat where they : are origi- 
en generated may b be thus con- 
sidered as caused by a “Magnus effect,” 


= similar i in n nature to the “ lift” of an 
in or to the propelling force in a rotor ship. The transverse motion, 
prompted by such induced pressure forces, combines with the drag of the 
current which tends to wash the eddy downstream. 
_ The transformation of the mechanical vortical energy int into 


by viscous attrition. Indeed, the | local rate. be expressed in 


terms: of Taylor f formula?*- 


which is a a numerical coefficient, and ) is a dimension of the 

small-size eddying group in the vorticing composite. ‘Inan eddying « composite, 

4 energy is dissipated i in large part through the instrumentality of the smallest o 

eddies, or the so-called micro-turbulence fringe. Large eddies, in turn, 
are responsible for producing the turbulent stress and more generally for actuat- oe 
ing convective transfer. The large eddies, accordingly, enhance turbulent 

diffusion, and their direct contribution to dissipation is comparatively small. 
It is surmised that, in the course of swarming, the larger eddies are subject to 

a & process of “ ‘grinding,” by v which they are split continuously into smaller _ : 


1888, Developments in Fluid Dynamics,” edited by 8. Goldstein, Clarendon Press, Oxford, 


va "Statistical Theory of Turbulence,” by G. I. Taylor, Proceedings, Royal Soc. of London, Series 
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4 units. ‘Taylor has treated this ‘general su 


: ‘The sequence involved in the mechanism of turbulent friction loss i is pre- 7 
sented ‘schematically i in Fig. 14. Of the total borrowed flow ¢ energy Wo, the 
amount W, is made available in the @ conversion zone. Energy W, is less than 
t 
| 


0 


; the cumulative loss |W,|+ in the transmittance zone. In the ab-_ 
a ‘Gute of adequate experimental knowledge, « one cannot judge as ~ what part 

of W, is dissipated locally i in the course of conversion. No o doubt, however 

_ (see Section 12), the conversion process is accompanied by a a considerable in- 
“herent loss. The latter, in Fig. 14, is symbolized by the ordinate CB. The 
-Temaining W (turbulent) = BA represents the part of the original W. which s sur 


vives conversion and which, liberated in eddying form, furnishes the ‘ ‘animating’ 


by the dissipation of mechanical energy into heat (Bq. 37 
ae An exact , quantitative appraisal of the sequence e would requ require experimental al 


15 shows, in part, the distribution of the axial and transversal turbulent com- : 


. 
reliably, i in three-dimensional the distribution of ‘turbulence in intensity | 
_ across the cross section. No complete data of this nature are available. Fig. J ir 


--:18**Some Recent Developments in the Study of Turbulence,”’ by G. I. Taylor, Proceedings, 5th ain 
national Cong. for Applied Mechanics, Cambridge, Mass., John Wiley & Sons, Inc., New York, N. Y. 

16**Turbulence Investigations at the National Bureau of Standards,” by Hugh L. L. ‘Dryden, ‘ibid. id., p. 36: 365. 


bject and has presented a 

— 
“source by which the state of turbulence, as revealed in its outward manilesta-- 
a tions, is actuated and maintained. Finally, curve B’—O portrays the gradual - 
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ponents in a n a rectangular conduit. ” As seen, the axial intensity y of turbulence, 


measured by 1 the root mean square | (w’)? , reaches a sharp maximum in n the — 
immediately, proximity of the wall, readily ‘corresponding to the outward 
_ boundary of the | conversion region. From \ there on, the intensity recedes 


"rapidly toward the center of the conduit. 


11. Errect or Watt RovuGHNness - 


basic premise underlying the Prandtl-von Karmdn theory i 
‘the e organic link between the stress structure and the shape of the velocit y 
‘curves.!8 Experiments - reveal that, if the roughness of the w: walls of a . given 


2 


conduit o of radius r, is modified consecutively, the fl flow Q ii adjusted to yield — 

;- each instance the same resistance slope S, and hence the same stress structure 
(Eqs. 3 and 4), then the respective velocity profiles i ins the central “turbulent” 


; Beitrag zum Turbulenz Symposium,” by L. Prandtl, Proceedings, ngs, Sth International Cong. for Ap 
lied Mechanics, Cambridge, Mass., John Wiley & Sons, Inc., New York, N. Y., 1939, p. 345. : 


ve 

z ey “The Mechanics of Turbulent Flow,” -— B. A. eames Princeton Univ. Press, Princeton, N. J., 
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gone will follow iene ‘shaped parallel outlines (see Fig. 16). In fact, if 

y the observed velocities 3s were plotted from a common maximum point, the curves 

_ the different roughnesses would ‘coincide over | ‘the gre greater part of the cross. 


seetion, the eventual differences sine apparent only i in the c com nparativ rely 


ee 


a In regard to the ener energy linen, dik an unvarying resistance slope makes the 
‘total energy | loss W. = y QS proportional to the respective average velocity. 
3 Furthermore, ‘congruent velocity “profiles signify that the velocity gr gradients — 
—— da/dy in the central turbulent zone are the same f for all degrees of roughness, = 
ain therefore in such § a zone the local losses w, by Eq. 34 will be neni in 


- The inescapable | conclusion is that the spending of « energy in the trans- 

mittance region, the latitude of. which 1 may be broadly said to coincide with — 
the central “turbulent” zone, is solely | dependent on the stress structure, and 
« that the spending process in the transmittance zone is implemented by the 
Bs ame agencies that shape the stress ‘patterns (see Section 12). . Deducting the 
7 common loss in the transmittance zone ne |W. |¥ from the respective W. o-quantities: 


hich are proportional to the var bulk velocities), one e obtains the energy 


quantities 


Ww hich for one or another degree of wall roughness are te a spent i in the con- 
version (boundary) zones. Paradoxically, the largest value of We occurs in 

the ¢ case of a smooth surface. ‘The rougher the walls, the s smaller the value of 
and therefore of W., always assuming : an unvarying- resistance slope. 

‘The The relation between | the velocities and the wall roughness i is ‘summarized : 
in terms of dimensionless ratios in Fig. (16.20.21 Velocity profiles: for p pipes of 
different size and different: relative Toughness are plotted versus r/r, in terms ‘? 


so-called “friction 


is the conv turbulent friction f factor 
U2 
4 and and f is the coefficient Pe ies in he pipe formula 
“*The Mechanics of Turbulent Flow,” by B. A. Bakhmeteff, Princeton Univ. Press, Princeton, N. 


20° ‘Strémungsgesetze in rauhen Rohren,”’ en,” by J. Nikuradse, Forschungsheft, Verein Deutscher 


‘*The of ‘Turbulent by B. B.A. B: Bakhmeteff, ‘Princeton Univ. ‘Press., Princeton, N. J. 
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ig. 16 16 to flow in pipes ranging from Li in. to 4 in. in 


¥ the walls of which were roughened artifically by uniform sand grains of different 
- _ The inverse ratios of relative annem ro/e are the quotients of the 


wane of the pipe by the eee size The shape of the roughness grains is as-- 


In the profiles in Fig. 16 congruence 
‘outline, which is to be anticipated, since reducing the velocities to the term ‘ty 

F is equivalent to referring the velocity pattern to a common stress structure. 
Since the curves are dimensionless and imply no restrictions with regard to 
the absolute dimensions of the conduit or the size of the roughness grains, the — 

‘diagram typifies universal velocity patterns for the “ standard” Géttingen sand 
Toughness type. Accordingly the experimental data, ‘endtelying Fig. 16 and 
explicitly tabulated i in the Gottingen Report,” may be used for disclosing, i in 
equally generalized manner, the ‘comparative features of the internal er energy 
exchanges for pipes of tt wd relativ wall In fact, in in 


terms of the “frietion 


“and the total energy loss and the av erage ge spending rat 


aa W.=2f To|To 


= 


Thea analysis m may be ye confined to the cumulative q quantities s only. In expre: ‘essing iB 
borrow replace the pressure gradient in w, (Eq. 11) 41; 


The : subintegral quantity is twice the static moment, with regard to the axis, 
of the area bounded by the velocity curve u/uy i in Fig. 16 16, between the limits 


(r/ro) in Fig. 16 by the u/us. 


After the curves" for the he borrowing and for the phases are. 
computed and plotted, the cumulative spending 1 rate 


‘Strémungsgesetze rauhen Rohren,” by J. Nikuradse, Forschungsheft, Verein Deutscher In- 
genieure, No. 361, 1933, p. 
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s determined by subtracting the ou outline for W; from that of Ws ieee, 43a 
43b contain the quantity 2 7 r, p uss which, with Eq. 42a, becomes to 


and 400 
curves. Th 


; and this value is used as a reference basis in Pree the ¢ c e 

vending, 17), of which Fig. 17(b) is ‘to 
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Fie. 17. —Comvtarive BALANCE IN TURBULENT FLow IN Pres 


terms, ar are most The reversal point i is closest to the boundary 


latitude of the conversion zone also | ‘grows; in the extreme case for which 


in the instance of a smooth wall. . As the roughness increases, the transv erse 
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r,/e = 30. 6, it ie about 10% of r To. As expected, v within he teeneniittenee 
zones, the cumulative loss is shown by a single outline. In fact, W, shows no 


appreciable discrepancies, v up to the closest observable r range (r = 0.98 ro) from 
the wall, which means that conversion is actually concentrated in a very narrow 
layer, the thickness of which is only a small portion of the reachry—fo. In the 


ease of the smooth pipe, the energy destined for conversion constitutes nineteen 


_ of a total of 25. 5, or about 75% Noe For the roughest surface, 


the corresponding energy is. 7 75 units of a total of 13.2, or 66%. On he 
’ hole, however, the dimensionless curves of Fig. 17 are | a rather compact set, 
suggesting that, over a great variety of circumstances, energy exchange p pat holy 


‘in turbulent flow exhibit generically 


THE NATURE OF THE SPENDING Process 


‘it is customary to distinguish betwee een the the ‘ “turbulent” 


p 


lent mixing.” One ait expect that hes agencies would be instrumental 
in “determining the manner in eons” withdrawn from the flow, is 
spent” in the different r regions of a cross section. 
Conversion Zone.— —Paradoxically the function of conversion, which was said 
be ‘synonymous | with the generation of “turbulence (Section 10) reverts 
primarily - to viscous action. _ The “casting of off” of eddies appears to be th the 
ultimate phase in a sequence of cyclic transformations suffered by viscous 
vortex sheets subjected to excessive strain. _ The intensity of the latter is 


yay that, in the vicinity of the wall, the boundary 


\ dy rin 


a> 


which exceeds by hundreds, if not by t hones the ordinary ‘ ‘laminar” strain 
underlying the appraisal of deformation work in Fig. 6. The extreme char- : 
acter of the ensuing deformation i is , suggested by Fig. 18(a); Fig 5. 18(b) shows in 
turn the ultimate phase | of the cyclic transformation, with rg original v viscous | 
- sheet rc rolled up into interspaced eddies, ready to be detached and “cast off” 


into the main flow." ul 


‘The deformation w work involved in the initial straining of the vortex sheet 
can be appraised by me means of Eqs. (> _ Indeed, it is viscous deformation work, 


- %“The Mechanics of Turbulent Flow,”’ by B. A. Bakhmeteff, Princeton Univ. Press, Princeton, N. 
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in the first instance, w hich absorbs: the energy made e available through | trans. 


In the case of laminar ‘motion, as stated in Section 7, all the deformation 
energy is dissipated in loco. _ Inthe generation of turbulence, on the other hand, 


reappears | in converted eddying 
form. Such an apparent difference in 
behavior ean be explained at. present 
% only by hy; pothetical “conjecture. The 
answer lie | the dimensional 
seales of ‘the respectiv e phenomena, 
laminar motion the action of viscos 
ity, 9 when one layer slides past another, 
as is on molecular scale. Any « eddies 
formed in the interface oe layers, which constitute the vortex sheets, 
- would be of ultramicroscopic s size and therefore. (Eq. 37) w ould become | ‘subject 
toa a most intensive dissipative action. Asa result the tiny “molecular” eddies 
extinguished in loco having a chance of breaking away as 


individual units; 


al By contrast, in the boundary zone of Fig. 18, the pendulations, induc ed by 
the instability of the initial sheet, result in concentrated lumps of vorticity. 


of a molar, macroscopic size. : In accordance with Eq. 37, the very scale of the 
eddies militates against instantaneous dissipation in loco, giving the energy 
iene invested into the vorticing units a certain period of life, before they 
are “ground up,” in the process of swarming, and finally extinguished. a 


‘The nature | of the phenomena i inv olved i in ‘the deformation w ork during the 


eonversion may be appraised quantitativ ely. Since, with | Eq. 2, the energy 


loss in Eq. 4 equals “2 A U= T. Pw U the energy available in the conversion 


zone, per unit unit of wall: surface, may be as 


ressed as 


in which a = W |W >, a fractional coefficient on the order of 0.8. | The rate of 
“deformation work per unit volume in the strained boundary den s 

ig. 18(a) appraised on the basis of Eq. 12b, is 


Since deformation w ork at such 1 rate is called to absorb the energy am amount indi- 
cated by Eq. 47, , the volume of the deformable layer (Fig. 18) must ‘measure. 
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Eq. 49a acquires a special significance when reduc ed to the terms of the so- 


‘ealled ‘ ‘friction distance” similitude parameter** ot ot 


mg @ For smooth surfaces, flow i is wholly shaped by the aforesaid similitude 

in § —, . Indeed the thickness of the laminar sublayer 6 corresponds to a 

ent parametrical value of y 8, and the transverse latitude of the entire bound- 
The ary zone, between the wall and t the fully turbulent inner region, to yy. ~ 30. 


‘The in-between reach (8 <y 30) constitutes a “transitory” gone. The 
physical significance of this zone, which lies outward fr from the strictly Taminar 


- sublayer and inward from the turbulent reach where flow is fully “ ‘infested”’ 
with eddies, would be that of a turbulence generating region, “within which 
“eddies are | formed and from which they are intermittently released. _ 
In Eq. 49a, the thickness of the deformable viscous sheet, required for 


absorbing the convertible energy W, in Eq. 47, corresponds to — 
For smooth | pipes, s, with es betwe een 0.0015 and 0.0025, _ a = 0.8, sta 


4 commensurate with the thickness of the boundary zone. About a third « of ya 
‘the | vin Eq. 50b is assignable to the laminar sublayer, in which dissipation occurs i in : 


in accord with Eq. (12d. This: leaves about two thirds, or slightly more, 


of W, to be s spent in the transitory zon zone for actual conversion into turbulence. 
- Again a portion of 1 this remaining energy must be dissipated locally in the course 
of rolling up the eddies, which means that only the smaller part of the flow 
energy W., first available in the conversion zone, finally Teappears i in turbulent 


’ region the contributions 

of the viscous shear in Eq. 45 are practically nil. Equa ally neglig rible are 


dy 


the energy losses , computed from Eq. 126 by using for —— the observed 


local values o of the main flow. ’ It 1 would be natural, on the other hand, to age 


the w,-rates, as they actually : appear in Figs. 17 and 11, in terms of the : same 
indi- “mixing” agencies, which | shape the “turbulence” stress. example, one 
sure may appraise the “mixing action” in the ‘ “Prandtl” manner,” assuming that a 


molar element (Fig. 19(a)), initially to Slament (1) and 
at a velocity u:, is moved by a transverse impulse over a “mixing distance 


‘00 te Mechanics of Turbulent Flow,” by B. A. Bakhmeteff, Princeton Univ. Press, | 
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into o filament | @, , where it becomes “embedded” acquiring the velocity U2 and 


Using the latter for u’ in the expression. for - Reynolds “ ‘ um 


_—assuming further that the transverse exchange intensity v’ is of the same order 


- and indeed is proportional to u’; and finally absorbing the proportionality 
factor by a proper adjustment of the mixing scale, one ends with 


in n which Lis the Prandtl “mixing length’ ’ and C i is the proportionality factor. 


Physically the change in the 1 veloe- 


T rupt transverse theme 
stances, in fact, are analogous to the 

case of impact with lasting constraints, 
as in the familiar instance of “‘inelastic 
bodies.” Energy losses for the latter are 
( to which | the loss suffered by 
element of mass m, the of 


aggregt 


ae exchanged per unit time across the upper and lower hase i is 2 p| \o’|, and 
‘the lost velocities equal ‘> dy) , thus making the energy loss (by ‘substituting 


w (Carnot) = 


reves revealing t that the “Carnot” loss, activity,” 
by the same » relation, which, in Eq. 34, determines the spending rate. In this 
light the quantity Ws (Eq. 34) in the contral turbulent transmittance region, 
may be qualified as 3 constituting a a “mixing” or a ‘ “stress-forming” loss. Note, 
in particular, that the velocities altered through attrition in ‘tesneveres en 


bedment are those of the main flow. Physically, therefore, mixing action ap- ap- 


de Mécanique Rationnelle,”” by Paul Emile Appell, Paris, France, 1911, Vol. II, p. | 508. 
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- pears to imply a direct dissipation of flow energy into molecular. bial arising — 
from the transporting : action of the large- -size eddies. © Such convective transfer _ 
~ of molar masses by the large vorticing units constitutes st specifically the function 
# so-called ‘ ‘turbulent diffusion.’ il _ Accordingly, the local | loss of flow energy 
in the inner turbulent zone of a cross section, determined by Eq. . 55, is a dissi- 
~ pation rate attendant upon “turbulent diffusion” and dependent in magnitude 


Brent atural Watercourse. river sufficiently broad and 


. pen so that the bottom traction force on at 


= 


init 


expressed by the customary du Boys formula 
ctor. A 7510 
Cumulative 


4444 4 
4 6) 
Velocity, in Ft per Sec co —z per Foot of width 


= 


In other words, the stress diagram 7 tis | the product of the multiplier 7S by 
the depth. observed velocity profile = = fy’ ) in Fig. 20(b) serves for 
plotting the cumulative discharge curve | OD i in Fig. -20(c), the ordinates of 
which (a"b’’) for represent the volume of flow 


indicated by the: shaded a area rea ABD’a’ in Fig. 20(b). 


In the light of Eqs. 7 and 9, whereas the total flow energy | ed and 


Spent in the unit of Fig. 20(a) is 


cumulative quantity—_ 


and 
(51) 
— 
(52) 
— 
rder * 
lim 
tl ef w — 
, ani a'b’e’d’ equal in 9 ane _ 
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_ —is directly represented by the cumulative | dischs urge curve OD in Fi ig. 20(c) 


to the scale—2. Likew ise, with E 18 and 5 57, the e energy ate cumulatively 


transmitted the block y across a’b’c’d’ is 


represented (except for a common 8) in Fig. 20(c) by the curve | 
OMC. . The ordinates of this curve constitute the product of the loc al veloci ‘ity 
by the respective depth, or the respective moment of the local velocity vector — 
regard to the surface of the streaming. The energy- curve is 

: _ obtained by s subtracting, point by point, W, from W,, making 


hh Eq. 61 the cnsiadiiy in parentheses represents the difference between a 


partial mean velocity Uy for the block (average ordinate of the shaded f figure), 


and the local velocity for 2 y’ (ordinate : a’b’ in | Fig. 20(b)). tie 


| 


Fie. 


nS The ete velocity profile, presented in Fig. 20(b), is for a lar 
_ern river during a moderately high stage. . With a slope of 0.35 Ry per ok. 


the multipler i in _ 57 to 61 is y S: = 62. 4X 0.66 x 10-4 = 0.0041; ‘and its 


In 21, the curves Fig. 20° e replotted in terms 


nh 


re typical and and shed considerable light 
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the eventual circumstances which the “w vorking” of. a watercourse. 
be fact, the reversal point is located at a depth of 69 ft t to 70 ft, or at about. 
11% to 12% of the total depth (78.3 ft) from the bottom. . Of the total energy 
W, (Eq. 59a), nearly § 92% is borrowed within the transmittance zone and of this 
: about 90% i is conv eyed by transmittance for spending i in the lower 8 ft to 9 ft 
of the flow. The aggregate transmittance loss is about 83% of the total. In 
z fact, the relativ vely small ‘ ‘slip” loss is suggested : a ‘priori by the shape ‘of the 
. velocity curv e. Inthe upper regions this curve exhibits only slight deviations 
from a fe sth nd actually begins to show | a pronounced decline of the local 
velocity only in the lower reaches, below 1 y= = 65 ft. Actually over a depth of 
10% of the total depth the cumulative loss i is only 2% of W.. It is less than 
4% for 0.8 y and 12% for 0.9 3 From there on, the spending i increases 
rapidly. a _ The most intensiv e exchange i is in the closest proximity to the bottom, 


the cgeniinnt in the last foot of the total y y ', = 78.3 ft constituting 51% of Wo. 
An analogous procedure applies tc to circular conduits. 


(14. GENERAL TREATMENT OF BOUNDARY Layens 


cof flow i ina boundary layer. A general account is given in ‘several texts.? .27,28,29 


It i as ssumed that the reader is familiar with the geners al concept and theory = 


te 


‘Premises that n motion” in the transverse | sense is s negligible, | and that the pres-_ 
sures in CTOSS Section do not differ from the pressure py = (x) at the border 
point b. In Fig. 22, and elsewhere in this p paper, the 1 reasoning is applied to a 

_tayer of unit width i in the direction perpendicular to the plane of the figure. _ 


The dynamic equilibrium f for. an element (Fig. - 22(6)) i is is expressed | by ‘equat-— 


-ing the mass-acceleration | quantity p dx dy — to the resultant of the : forces, of 


| _®“Fundamentals of Hydro- and Aeromechanics,” based on lectures of L. Prandtl, by O. G. Tietjens, 
translated by L. Rosenhead, Ist Ed., MeGraw-Hill Book Co., Inc., New York and London, 1934, p. 58. 


Mechanics of Viscous Fluids,”’ Div. G in ‘‘Aerodynamic by L. Prandtl, W. F. Durand, 
q in chief, Springer, Berlin, Germany, 1935, Vol. III, p. 


I * “Fluid Mechanics for Hydraulic Engineers,” by Hunter Rouse, ‘Ist Ed., McGraw-Hill Book ¢ C 
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dps 


whiel h- - * dx dy is the pressure component, and pod dy dx is the friction < 


sultant. Reduced to a unit volume of fluid, the equation oe - 


Forth the usual case of up = = constant, and 


‘indicating that, in the ia of : a pressure gradient, n motion in a the boundary 


layer is ; sustained by the inertia of the | oncoming i. The negative friction 
resultant — (see Eqs. 10) is offset by an re retarding 0 of the motion, 


with the fnertial force” the negative friction 


Obv ously work implied i in overcoming frictional 1 resistances i is 
ae at the expense of the kinetic | energy of the flow, the initial store of 
_ which is gradually reduced. With regard to an elementary cross-sectional 


block (Fig. 22(c)), the the kinetic energy Ex, entering time across ab, 


ceeds the energy leaving the block across line a’ ad by lade, in which the 


quotient measures the energy withdrawn and spent per unit 
time in a boundary layer block of unit length. 4 It has the same significance for 


boundary layer flow as W, (Eq. 4) has for established uniform motion, except 


that 4 ian vary from section to section. Since the quotient is essentially 
dE, 


the inner subscript 6 suggesting that Ws covers the entire transverse extent of 
_ the layer, whereas the outer subscript | indicates specifically the - location of the 

particular unit block under ¢ consideration. On the other hand, [Wale “will 
designate the aggregate energy borrowed and spent in a unit of time in the 
‘body of the boundary layer OBA between O and z. Obviously for 


0, |W a|z is is the difference b betw: een the kinetic energy of dy flow- 


a 


ing across AB equal to the energy p 2 5 dy, originally. carried d by the same 


Vv ‘olume of flow i in n the « on- n-coming u state, across O-O. In in other words, 


= 
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‘tac E, includes energy carried across the cross section from the 
to infinity. Inasmuch : as there is no appreciable change of velocity 
(0 nae. of the boundary layer, one may subs substitute for the right side of Eq. 05: 
( ay] dy. well- known von Kirmén 


= 


_ which serves to determine the skin friction resultant 0 over the length z in terms 


of the momentum lost by the flow between O and z. _ Note that Eqs. 65 and : 


66 are quite general : and therefore are eummaned to boundary —" of all types 


‘whether laminar or turbulent, subject only to ‘the limitation 


fortunately, present know ledge of turbulent boundary layers is anything 
; but complete, and an adequate ana analys' sis of the energy y exchanges i is not possible. 


On the other hand,. for the case of a laminar layer next to a flat plate with | 


0, detailed theoretical calculations are available, and are in excellent 
accord with repeated ieneniibinias, - For example, theoretical and experimental 
aspects | have been presented in detail by Hugh L. Dryden,® whose ‘report, as 
well a as subsequent papers by the same author, 16 contain the r more re- 
cent nenpneone of boundary layer motion, particularly with regard to 


“transition.” 


5. Energy BALANcE IN LAMINAR BounparRy LAYERS WITH | 
o& Basically, the the flow flow pattern is conditioned d by the parameter 


v 
27 


hich R, is the so-called “distance Reynolds of the boun nd- 
er. In particular, the transverse velocity = f(y) at a station 


80 that the ratio ¢ of a local ‘velocity « u distant y from the plate, to the con- 
same 8 stant (border) v velocity 2 to | & is ‘solely a function of 7. In any two 


sen to make = (2) VEG and 

Toss sections, if the ordinates were. chosen to make m = ¥ Ra and 


ol *“Air Flow in the Boundary Layer Near a ay ” by Hugh L. Dryden, Technical Report No. 562, : 
’ ational Advisory Committee for Aeronautics, U. . 8. Government Printing Office, Washington, D.C., ‘ 
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= (#) vi the re respective velocity ratios = — an 


a would be the same. _ Velocity profiles with such similarity of shape are are said to. 
a. be ‘ ‘affine,” ”” per rmitting the representation | of the u = = fly)- -curve fo for any cross_ 


Values of d 
¥ 4 
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tiating 


Affinity « obviously also to the derivative quantities 
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in Fig. 23. 


ENERGY 
he 3 for which are also presented in Fi ase on a 


values of —— 


curve 4 In par particular the dotted line D-D in Fig. 23. 
is drawn at a transversal = 5. 2, by as a conven-— 
tional n measure for the boundary layer thi . At such 
tance, the local value of wu differs from the — re Uo “el the negligible 


amount of 0.5%. Fig. 23 discloses, however, that the derivative values at 
= 5.2 are still quite substantial, and may weigh accordingly in the eventual 
energy balances. A more appropriate conventional v alue for laminar 
layer thickness would be ns = 6.9 making 09 tle 


value of 75 = 6.9 is chosen, SO as to give 6 four times so- ‘dis- 


‘placement thickness” 6a = 1.7 The Pre andtl ns = 5.2 corresponds to 

= 
Bora “affine” ’ profiles, the energy e a: s take definit definite numerical form. _ In J 


fact, substituting ¢ Up for u, ,and 7 ” for y, the aggregate energy Ic loss (E 


Fas 4. 


» 9 

“outline of Fig. 23, the w 1, 044. 
To obtain the local cross- -sectional loss (Eq q. 64), « one proceeds by differen- 

tinting Eq. 70, which makes 


For establishing a interrelation between the energy loss ; and the skin fric- 
tion resultant, note that, numerically, the von Kaérman expression (Eq. 66) 


4 3“*Air Flow in the Boundary Layer Near a Plate,” by Hugh L. Deyden, Technical Report No. 562, 
National Advisory Committee for Aeronautics, U. 8S. "Government Printing Office, Washington, Oe 


4 
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in which w is a dimens 
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and the local w wall stress ses fore se per unit length of plate at the 


16. THe Locau Unir Rat 
Local borrowing in a boundary layer is aes by the same formula 
| 


(Eq. 11) as in uniform established atin although, with = 


: absorbed b by y the friction  tractions is supplied wholly from the kinetic | energy 
3 


In particular, for a a laminar layer, using Eq. 69 in Kq. 11, one one obtains 
us, 


7 Furthermore, ‘dividing sina tiia in Eq. 72 by the transverse dimension of the layer beyond 


— (see paragraph ¢ mnning Kq. 69), one obtains an average w unit ~The 


Eq. 79 into Eq. Ti: “ 


— 
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borrow 

‘Substituting Eq. 75 into Eqs. 70 and 72, one finds 

* 

<2 

multiplier of the parentheses group constituting (in view of Eq. 72): _ 

4 

= 69.4) 

rateof 

| 

= 

— = (ta), 

— 


7 


= 
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is a dimensionless function, represented i in Fig. 24 by oe curve AaO. This 
curve is obtained by multiplying the ordinates of the curves Ao and } Bo : 


Mig. 23. ‘Obviously, € ',, in terms of 26.4 (ws)z, denotes the distribution o of the 

borrowing process in all laminar boundary layers of the type te ~ @. Tt 


[TTA 


Fic. 24 


local borrowing i is zero. near the wall, and it inp declines. into insignificance 


beyond 6. 9. maximum of the borrowing function is reached | 


The Spending Rate. a laminar ayer, the local spending measured 
by the work involved in viscous us deformation (Eq. q. 126), one obtains, with 


a) =f - ( = 26.4 (ws) 2 e’,...(82) 


dy/ dn 
Again, the dimensionless function 


represented Fig. ‘as the curve AS, describes spending process 


generalized universal terms. — Similar to uniform motion, the spending of energy 
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is cc concentrated near the solid boundary. = the intersec- 


point M,. corresponding to = Ws, indicates the region within whichf 
‘spending exceeds the energy amounts by local borrowing. 
_ 8 The Transmittance Rate.—By subtracting curve AS (Fig. 24) from curve ‘thatis 
one obtains the outline ¢ curve or which defines the local trans-f 


_ Expressed 


4 


rs = W — Ws = 
b v8 


wm ulative Energy Rates.— presenting the 

tive energy quantities, it is expedient to add in the direc. 

from the unaffected outward zone the ate, 


3.83 (Ws)s 


in which is a dimensionless 


_ obtained from the data in Fig. g. 23, an 
ically the quantity | Wile indicates the performed on the block B, Fi ig. 
by the e overlyi ing fluid prism nA. "7 The other cumulative quantities take the form: 


L 


(Ws dn = 3.83 (Ws)2 


in hich the g e generalized functions and , are equal t to the 
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The Energy E quation. —By Eq. 17, the local energy 


eb ay" 
{ OH 
(- 
e, 
er 
day tis 


= 
= 


cor 


The: nin matted formula la expressing the relation between the cumulative 


of energ s attendant to fluid friction, the “borrow- 
ing,” or withdrawing of energy from the flow, and its spending do not coincide: 
; spacially. "Spending i is concentrated in the regions s of high vorticity near the 
“ae boundaries, whereas the borrowing takes effect principally in the central 


‘The aforesaid implies an in-between “‘transmittance”’ function, by which 


surplus borrow ed energy is transmitted to the wall region, where local borrowing 


is Is insufficient to cover the requirements for. spending. 7 
= 
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¥ 
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i 3. The aspects s of the three consequent phases are determined by the internal 
4 stress structure and the concomitant velocity profile. 
dy Tat dy* 


applies equally to laminar and turbulent patterns, and serve — —_ for 
the reasoning presented i in the paper. 
For laminar’ r patterns the energy exchange pr process can 
wholly i in analytical form. | Tn turbulent patterns the application tte equa- 
tions must ; rely on detailed data. 
6 ‘The outstanding pines of all turbulent patterns | is that “only. a 
pe and that the larger part (30%, to 85%) is concentrated by transmittance 
for spending in a narrow reach near the solid boundary. This “conversion” 
7 - is ostensibly the seat of the pivotal process in the turbulent mechanism— 
4 namely, that of converting flow energy, concentrated by transmittance, into 
turbulent eddying form. 


The function: conversion ‘in itself reverts to viscous action. The 
part of the « energy, is ostensibly dissipated 
in loco i in th the viscous deformation of the laminar s sub- -layer | and i in the ¢ course 
of forming eddies in the adjacent turbulence generating zone. Thus only a 2 
_ smaller part of the spent flow energy reappears after conversion in the turbulent 
form which animates the familiar turbulent manifestations and activates the 
* convective mass transfer instrumental in shaping the turbulent stress structure. 
a The e final dissipation of of the turbulent energy, continuously replenished 
. from the generating boundary ary zone, is s by way way of viscous ¢ attrition in the course 
of the swarming motion of the eddies across the main flow. 
8. The direct local loss of flow ‘energy in in the central usteniiitenes ian zone can 
_ accounted for dynamically by applying the Carnot impact law to the process _ 
- of mass embedment in the course of transverse mixing. This loss, attributable 
to the convective action of the larger eddies which shape the stress structure, 
is concomitant to turbulent diffusion, 
- 9. In relation to atural watercourses, | the transmittance of energy and its 
spending i in the bounuary zone, shed possible light on the “working”’ of a river 
ke 10. The analysis of of boundary layer motion is. extended to ) consider the energy 3 
_aspeets of the flow patterns. . In the case of a laminar layer v with 1 dp/dx = 0, 


a complete account of the energy exchange mechanism is revealed in generalized 


 Acknowledgments.—The subject matter covered by the paper forms part of 


a + broader investigation on n the “Nature: of Hydraulic F Resistances” conducted 


F or the velocity profile used i in Fig. 20 the writers are gratefully ‘indebted to 
“the Mississippi River Commission of Vicksburg, Miss. _ 
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A PLAN FOR A MOVABLE E DAM 


‘Br Isaac DEYOUNG," AM. Soc. Cc. 


Pilea are two general conditions under which en dams are operated: 
One i in which it is requ red to close a channel i in an orderly predetermined plan 
and Ww ithout i1 immediate ‘urgency; and the other, in which a sudden release of the 


water in a canal caused d by an accident. due to ‘the w recking of a lock gate by 


collision or otherwise, requires a structure that can be operated 3 in . the rushing © 


OF OUNCTY sili 


waters, checking the flow. 
The plan fora movable presented is of the latter type. It is 
an naan dam, and is designed s« o that it can be put in place during the rush | 


held in place by rigid hinged to the 


under the floor of the lock. The dam i is raised through multiple reaved cables 
by two hoisting | engines. The raising operation required only” a short time. 
A single o operator i in a central station manipulates the controls i in an orderly way 


¥ 


as soon as the upper canal is free from all floating craft. wd 


The movable dam cnn in n this paper is a development from the studies — 
made in the e early 1900’s in connection with the design and construction of two | 
movable dams for St. Marys Falls Canal, Michigan, in which the writer had a 
“part. The writer’ ideas: leading to the proposed dam have been 


through - the years intervening since that first experience. 


Forms or Movasie Dams 


‘conjunction with wickets, 1 needles, : and other: (2) submerged 


a Nors.—Written comments are invited for immediate publication; to insure publication the last dis- 
(cussion should be submitted by July 1,1945, 
: on ‘Consultant to City Mer., Sault Ste. Marie, Mich.; formerly Supt. of St. Marys Falls Canal, Michigan — 
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(3) floats; (4) lock gates (sheathed « or skeleton), sliding | 


segmental: gates, and lift gates. 
(1) Overhead bridges of the swing type, from which either horizontal or 
“vet stop planks, wickets, or needles are manipulated, are fairly sure to be 
operated so that the flow of water is checked gradually, thus minimizing g the 

~ shock on all parts of the structure. The operating parts are readily accessible 
and easily repaired. Some of the disadvantages are that they are limited as 


of that a sunken vessel in the may prevent 


a (2) yo dams of the utamean type m: may be mentioned trestle dams, 
bear traps, Taintor gates, ete. The flexible shutter type of dam described in 

| this paper belongs in this class (see Figs. 1 and I 2). _ Some inherent objections 

va 

to the ty = submerged dam are that a dragging anchor m may injure the struc- 


into at jaws of the narrowest part of the canal or lock. 
_ (4) Having two sets of operating gates at each ¢ end of the lock offers some 
“safety because of the double barrier in front of a a moving ng vessel. . A gate located 
in the canal above the lock would accomplish a similar result, the gate being 
operated i in in such @ manner as to be alw ays closed during a danger period when > 


the: lock gates : are open; but the « operations 1 under such conditions 1 would be 


to excessive e delays to navigation. 


~The structure ¢ described in this paper i is composed of 
tome together on the downstream side by plate hinges extending the entire 
= the girders (Figs. 2(a), 3, 4, 5, and 6). On the side the 


tg position, the is folded below the lock floor 
(wee Figs. 3 and 7). — is designed for a lock 80 ft w ide, having 30 ft of water on | 


‘mee (a) Seeae. of operation when required; (b) reasonable facility 
operation ; (©) no obstruction: to navigation; (d) ability to withstand shock; 
simplicity of construction; (f) accessibility for inspection and repairs; 
4 (9) appearance; ; (h) freedom from ene enemy attack : during str stress of war; and (7) cost 


is (a) A movable | dam should be free from the pat ees of being clogged or : 
broken by the —— use of the canal or by the exigencies arising at the time | 


of the accident. 


— 
— 
— 
— 
— 
7 4 
» | 
[These objections, as applied to the flexible shutter dam, are discussed 
subsequently, 
(3) Under the subject of floats, there is the ed by sinking a 
} 
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~ 4 Part Cable Reeved 
Over Equalizer Sheaves 
Each Endof Dam 


Cable Clamp — 


4 Patt Cable’ 


6 Part Lifting Cable 
Side 
of Wall Recess 
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In the proposed plan the dam is folded i in small compass in the bottom of : 
the lock and protected by transverse breast walls, prev enting da amage ii 7 
dragging equipment, such asanchors. 
| 


SY 


+) 
‘SECTION D-D 
= Fig. 4.—Sections THRovGH Hinces 
As debris may be moved by the current created i in the canal when a lock gate 
4 is w recked, some provision should be made to prevent the collection of such 
. debris on. parts of the dam, so as to to interfere with its successful il operation. — hh 


_ the plan, the anchor slots, eight in number, are covered by the top anchor having | 


= 


= 


‘ode iron attached, whose outstanding legs cover the slots completely from end 
4 o end of anchors, preventing the debris. from entering these slots (see Fig. 6). 
pny funnel-shaped | pits betwee een the anchor slot walls catch this debris. = 

It is not considered likely that a vessel would become stuck, because then 


it would itself become a barrier, building | up pressure sufficient to move it 
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It would be desirable if the dam could -™ closed before the current in the | 
canal ha has reached a dangerous velocity, but this ideal is impossible i in any form 
of dam. . The flexible shutter dam can be raised in a few minutes, and as soon 
: as all floating craft in the upper canal has passed through. 
i (b) The operating m mechanism should be simple, | not likely to get out of 
“order, and it should be easily accessible for quick operation . The flexible 

- shutter da dam is operated by n means of two hoisting engines, one on each side. ad 
lock. The operation may be performed from one ‘station through appro- 
priate electrical connections similar to the Selsyn method of operation. . (The 
- Selayn system of of signaling and control consists of two s small ‘specially wound 

-induetion motors electrically interconnected in such a w yay that the movement — 
of one armature causes the other armature to move by a corresponding amount 
and i in the same direction. One motor is operated at the sending point as a 
generator is called the transmitter and other is as a motor and 


noe width, the operation may be performed eure through a train of nl 
(ec) The structure re should be of such a form as to offer no obstruction to 
navigation. These requirements have been met in the flexible shutter dam in 
that all parts are outside of the limits of the navigable cross section of the lock. 
a 7 @ Since the proposed structure is held rigidly in place by substantial anchor 
beams, which revolve slowly about th their centers, there will be no sudden change 


of quien @ of an: any part « of the dam to receive any s shock. _ The flow will be 
- stopped gradually and there will be a smoothness of operation. — ‘ It is expected 
that the apg not be raised until the upper part of the canal has been 


any floating craft that have been in 


proper is ‘composed of plate girders. The anchors for. holding the 
of H-sections, channels, angles, and eyebars, _ The hinges are of steel plat ates, 
continuous on the downstream side the full width of the lock. , These hinges 2 are 
Tiveted | to the f flanges of the ts = with the hinge pins in place; the Ll 


side, holding the dam, are riveted to the w vebs of the plate girders. All the 
‘necessary that these because 1 of the parts is slow 
and infrequent. The upstream flanges” of the plate girders consist of angles 


having their meeting legs machined, a lose and joint. The 


The flexible shutter dunn is within the lock betw een 
the upper guard gate and the upper operating gate at the elevation of oft the canal 
bottom of the upper ‘approach. When the lock is unwatered it can be com 
pletely inspected, painted, and repaired, if required. 
@ As all parts will be placed out of view, the proposed dam has no ob- 
 jectionable features from an appearance point of view. 
_ (h) Because the structure is submerged, it i is not a likely subject for e1 for enemy 


= attack. The power plant can also be placed out of view, if 
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As the flexible shutter dami is is light ix in other types 
of movable | dams, the cost will be considerably less. Cost of maintenance is 
small because of the smaller | number of parts and the ease of access and — 
when the lock is unwatered. The weights (in pounds) of parts of the 
proposed dam (lock 80 ft wide and 30 ft deep) are estimated semana as. 


Dam proper, including girders, connections, 00 of anchors 310, 000 
Other half of anchors (supported o on lock floor 53 ,000 


7,000 
"Embedded anchorages 40 ,000 


425 ,000 
The w eight of two hoisting engines has not been included in the foregoing 

estimate. Because of abnormal wartime conditions, ‘no estimate of cost is 


presented. Under normal conditions the a1 average cost of steel and iron in 1 the 
structure should not exceed 15¢ per lb. 


_ Probably no engineering structure has been g given more study in connection 

_ with the operation of locks than a movabledam. Should a lock gate be wrecked 
7 by collision or otherwise w hile restraining the waters in a canal, , sudden release 
of the water would cause a wave that would probably result in carrying g away all 

=_— moored at the canal walls, and such vessels would be e subject to collision 

below the locks and be wrecked. The magnitude of the catastrophe wot ould 
_ depend on the head depth and size of the lock and canal. Havitig a movable 

dam will not avert the wrecking of a lock gate, but such dam will serve to to stop 

_ the flow of water, permitting r repairs to be made and restoring the lock for service. 
pe The importance and necessity of having a movable dam i is evident. The 


“task of stopping a flow of water in a . ship canal without such a structure w vould 
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AMERICAN SOCIETY OF. Civil. ENGINEERS 
Founded November 5, 1852 


REPORTS 


4 
ORGANIZATION, FINANCING, AND 


PROGRESS REPORT OF THE COMMITTEE OF F THE 
SANITARY ENGINEERING — line 


INTRODUCTION 


vou « 


—This progress report is intended to stimulate discussion. 
The« objective of the committee is to formulate the fundamental concepts of the 


subject assigned to it. pupa provide a basis for the preparation of state 


renin acts or for the amendment of existing enabling acts regulating the 


organization, the financing, and the administration of sanitary districts. . Iti is 
hoped by the committee that there will ill be an opportunity during 1946 to receive 
and consider discussions of their progress report by engineers, administrative 

- officers, and lawyers, and to present a final_report in January, 1946, that will 
be ca to the Local Sections of the Society, 
2. General Statement. —A sanitary district, is a local al governmental unit 
created under a state enabling act to: finance, construct, and operate works for 
the collection and treatment of si sewage, , including industrial wastes amenable 
7 to joint treatment, and to provide adequate final disposal of the effluent and 
of the ‘by-products of the treatment processes. The principal Parts of the 
works to be built and administered by a sanitary di , district include: © . aborted 
(a) House « connections, mainly from the street sewer to the curb; on 

(b) Sewers (lateral or collecting, trunk, and intercepting) 


The three parts of the title—organization, financing, and administration— 


are interrelated. All spring initially from state enabling acts which give 
wide or only limited powers to the administering bodies that they. create. 

Sanitary districts differ in size and in their objectives. Some undertake 
only the construction of intercepting sew vers, m main pumping stations, and treat- 

ment works and may serve a considerable number of municipalities. A Some are 

_ coextensive with an existing city. A few include lateral or collecting : sewers 


Nore. ritten comments are invited for immediate to insure — ‘the last 
should be submitted by July 1, 1948. 


3 
| 
— 
— 
— 
| — 
4 
d — 
a 
| 
| 
| 4 
A H 
4 = 
= 


17 SANITARY "DISTRICTS 
in their construction, whereas others take over existing collecting sewers for 
0 operation : and | extension, or leave all of this to existing municipalities within 
boundaries of the sanitary district. ‘Thus, local conditions and practices” 
- areimportant and will affect many procedures as to administration and nd finance. | 
Such physical and social differences explain i in part the different procedures” 
that have been adopted in the United States. 
_ Inthe past, local governmental units have been | giv: en the power to accom- 
Te one or all of these objectives and these units have been created as taxing — 
units with the power, in general, to levy ad valorem taxes and special assess- 
4 ‘ments. This grouping would include s sanitary districts with taxing power and, 
in some cases, with power to collect charges for services s only, the taxing power 
= denied. In a few instances these powers have been combined. The | ‘ 
scope of legislative neliion is limited or prescribed, in some instances, by con 
_Stitutional limitations upon the power to issue bonds and upon provisions 
- relating to the levy and collection of taxes and special assessments, 
8. Present Status. -—There i is a considerable difference of thought and prac- 


tice regarding these ‘matters by courts, commissions, governmental authorities, 


q 
-relativel rely 1 new sound precedents are to In some cases, 
Ww ell-established practices cannot or should not be disturbed. _ The many — 


different enabling acts and methods of financing are ample evidence of this 


situation. here ‘competent: professional service has been used, the design, 


constructian, and operation of of sewage dis disposal projects i in sanitary: districts 
have been : ably accomplished. 


The fundamental considerations leading to the best methods of organizing, 


i financing, and administering sanitary districts are the subject o of this gs 
_ The present status of sanitary districts, as a matter of historical record, i 
ion in several publications, and reference is therefore made to the Bibl 


PART 1. 


as organization of a sanitary district, as 
expressed i in in the following 


 @4H How shall the district ct be originated and 
— (b) How shall the boundaries of the district be ew 
(c) How shall the. administrative board be formed? 
@ How shall the | district be financed? 
aa [On ‘What shall be the powers of the district? 


The answers to these questions constitute the fundamental steps v shich 


“more or less control the p preparation of an enabling act. In general, the « com- 


district should be taken Ww within the locality that ae comprise ‘the sanitary 


Origination and Creation—Under existing enabling acts, 
district is generally "originated and created in one of the following ways: _ 
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a) ) By a petition from a certain number (often 100 to 500) ‘qualified ~~ ers” 
(b) By the a action of one or more existing local governing bodies, su such as as 
city or village councils, county s supervisors or tow nship. or 
(c) By the action of the state through its legislature or some designated 
department of the state, such as the Department of Public Health. | 
—— committee rere the origination a and creation 1 of sanitary districts by | 
d by a a general vote by - qualified p persons within the desig-- 
nated boundaries. The petition should be presented to a duly elected local ] 
official or governing body who shall consider it, hold public hearings, approve 
or disapprove the action, and conduct the vote. The emphasis in the first 
place should | be on local action. 
oe In cases where the need for a sanitary district is great and the locality does 
not make sufficient progress, tl] the state should notify the the locality of this fact and 
—after a a reasonable period of time, perhaps one year—should create the dis- 
trict under specific authority granted by legislative action. 
7 6. Extent of Sanitary Districts—The boundary of a sanitary district is 


important. Districts having the following found: 


(a) Part of a municipality; 
An entire single municipality, y, generally a ie 
A ‘single » municipality | and some adjoining unincorporated area; = 
A group | of municipalities, generally adjoining, with | some unincorpo-_ 
rated area, the whole. constituting a metropolitan area; 
(e) A single county or an area within a county; 


f). ‘An area within an arbitrary line, ¢ comprising a region, and — 


(q) A river valley or drainage area. - 


The committee Poses the determination of. the extent by a competent 


the locality. 7 The ¢ expenses ; of this survey and peony preliminary costs prefer- 
ably should be appropriated by the existing local greene bodies affected. 
ated for part or all of the costs by the state. . These. costs m may later be repaid 
by the sanité ary district if and when it is created, organized, and financed. 

a  . it becomes ‘necessary for the state te come the district, the legislature 
- should appropriate moneys for the preliminary engineering and administrative 


costs. 


| jae Administrative Boards.—Enabling acts provide several methods of nities 


up the administrative board, among which are the following: — 


(by | an existing local governing body. or person such a as 
Jocal judge, a mayor, a munici sal council; or by the governor of 
© By designation i in the enabling : act 1 as, for instance, the e appointment of 
certain existing local officers elected. — 
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SANITARY DISTRIC 


Some enabling acts provide for a combination of the foregoing methods. 
7 he committee favors, in general, the election locally of the members of the 
: administrative board. © _ Many appointed boards have been eminently successf ul, 
a however, and the e appointment of board members by some locally elected 
official may be acceptable method. Where municipalities are 
7 included, the chairnian or president of the sanitary district might be appointed 
by the governor from candidates nc not residing i in the « district. _ Experience also” 
indicates that, except in very large districts, the board members should receive 
- only nominal compensation. In In very large or ‘special districts th the chairman 
of the board should receive an adequate salary for full-time service and, i in 
some cases, each member should be so en A board of three to 
re members i is considered to be w workable. 
«8. District Financing. —WwW hen a district is pany it should automatically 
_ collect a charge on the area W ithin its boundaries to 0 cover the reimbursement of 
the preliminary ‘expenses a and the e payment « of the ‘administrative, engineering, 
and other expenses of the district prior | to the receipt of moneys through the 
Pow ers of a District. —The powers given to a sanitary ‘district the 


~ enabling act are not unusual; they are details of secondary importance to the 


_ fundamental considerations. The district should have the pow er to: sd 


(a) Make and for the use of its facilities; 


(d) Preseribe the disposal of objectionable ‘industrial 
not amenable to tr eatment with domestic sew ‘age, and regulate the 
om ( (e) “Sel or otherwi wise dispose « of effluent, ‘sludge, or other er by-pro jodacts; a 
Require the use of its facilities whenever they are accessible; | 
a (g) Make and apply rates of charge for sewage disposal service and dete b 


(h) Make use of suitable patented processes necessary to the work of the urb 


Purchase material, supplies, and equipment without ‘competitive or 
a formal bids up to a stated total cost for each transaction, the amount a 3 


= rie to be larger for large « districts | than for small districts. owe At 


Parr II. Finance 


ae money roe cover all of its proper r costs, including co1 construction, operation. 
- maintenance, and administration. _ These can be considered as the total annual 
- cost, which includes debt service and all the current expenses of operation. 
This part of the subject reaches into the fundamental ec equities of the allocation 
of the over-all costs of providing sewerage or sewage ‘disposal facilities. No 
single method in the United States has 1 receiv ed predominant support by courts 
commissions. variety ‘of methods and combinations of methods call 
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SANITARY DISTRICTS 


be cited. . This diversity is caused, in in part, by the historical aspects of methods 


of financing sew age disposal projects as follows: 


 @ One of the ‘earliest methods was to raise all necessary funds by taxes 


on ‘on property. This was accomplished by the issuance of general obligation | 
_ bonds supported } by property taxes and by an annual levy to support operating © 


_ (b) A second method i in ‘point « of time has been to raise construction costs 
by taxes on property and operating costs by charges for sewage disposal se service. 
‘This method h has been used toa considerable ex extent in Ohio aml will be found to 
apply to many present sewage disposal projects. 
more recent method has been to raise all ‘funds by charges 
to) persons or property receiving service. The distinguishing characteristic of 
this method is that the entire cost of the sewage disposal service is allocated to 
4 persons or property (or both) | on w hose account it is provided. The effort is 
ey made to allocate the total annual cost on the basis of the cost of providing the © 
_ (d) For many years, the cost of constructing sewers has been met by a2 
frontage assessment, such as a charge per front foot of benefited property. 
ay Methods of financing to be adopted should be considered in the light o of 


local conditions, such as established methods and habits, public relations, con- 
a venience of administration, enforcibility of collections, and the security of the 
oa rev enue from year to year. The extent of the work to be built and operated. 


i will also influence the method of financing. : Thus, practic e and experience a are 

‘The foregoing methods of financing have different results as regards the 

of ‘If it has been or is fair to rs all the costs a 


“funds by a  chazge against direct users; that is is actual connections 
with the sewage disposal facilities. In some districts, such as the Brockton 
(Mass.) and the Buffalo (N.Y sanitary districts, and the Washington Sub- 
Sanitary District in Maryland, the methods of financing seem to 
"approach equity in 1 that an attempt is made to allocate the costs to | persons and. 
‘Property reasonably in. in accordance with the cc costs of prov iding the ‘serv vice. 
re _ As between different methods, the sums involved should be considered. 
A total average annual cost of $5.00 per capita amounts to about 1. 4¢ per 
‘eapita per day. - This is a small unit and in some cases a moderate departure 
from the fairest charge in favor of a well-established aed convenient to 
administer seems to be tolerable pr actice, 
‘Methods of financing in present use might be tested to determine how 
closely they result ‘in an equitable distribution of the charges n necessary to 


defray the total annual costs. It seems unlikely that an altogether equitable _ 


method of financing can can be set 1 up that is not too o complicated, and considera-— 


tion must be given to what is convenient and expedient. Iti is believed, how-_ 


-ever, that substantial equ equity can be accomplished through methods of financing . © 
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including cl ecllocting oon sewers as well as intercepting » sewers sewage treatment 


works, is taken from at an ‘unpublished oe by Frank. A. Barbour, M. Am. Soe. 
te * General taxes are based on the valuation — property or on the 
ability of the individual to pay, regardless of the par ticular benefit derived 
by him from the expenditure of the money so obtained. Thus the cost of 
education, of policing of the community, the support of the poor and other — 
; public necessities, the benefit of which to the individual cannot be definitely 
_ estimated, i is reasonably provided for. An assessment, on the other hand, 
is levied in proportion to the particular benefit accruing to individual parcels 
of land by the construction of a certain public utility. Unlike a tax it_ 
“ does not take something without any assurance of compensation but rather | 
is intended to express, with as much accuracy as possible, the value added to 
the property by the works to meet the cost for which it is levied. An 
_ assessment implies permanency in the improvement and in the value added 
to the property. Obviously it must be apportioned on some logical, uni- 
{ form and unchangeable basis such as the physical dimensions of the property 
> and the rate per unit once established cannot be altered from time to time. 
The value added to the property may be merely potential and does not > 
imply actual use of the public utility. Such use will constitute a further — 
benefit which in amount may be estimated by rental charges based on | 
units involving the elements of quantity and time and which are not 
necessarily fixed in amount as in the case of assessments. 
a “Having thus called attention to the essential characteristics of general 
taxation, assessments and rentals, it remains to be said that various com- 
binations of these methods of raising money for sewer purposes have been 
used in different municipalities and, as will be shown later, with statutory 
eee construction of a sewerage system results| in a a general hygienic 
betterment to the community from which must follow an indirect return 
5 to all citizens regardless of actual use or opportunity to use the system. 
Further, in all sewerage works a portion of the expense for outfalls and 
- purification works is due to topographical conditions or sanitary obligations 
4 to neighboring communities, which constitute a disability to be borne in 
4 some considerable part by the Town as a whole. It therefore follows that 
a portion of the money necessary for the installation of a sewerage system 
is logically raised through the general tax rate and, as will be seen by 
reference to the statutes, the law requires that at least one- quarter (1/4) 
3 not more than two-thirds (2/3) of the total cost of the system must be 
so obtained. This is equivalent to a statement that at least one-third (1/3) 
; ‘a be collected on some other basis than through the general tax. — The 
az is to determine what proportion of the whole cost between one- 
third (1/3) and three- -quarters (3/4) shall be thus raised and how can the 
amount to be so raised be most equitably apportioned. © Obviously the 
answer will be found in that method which will most fairly express tl aa 
‘relative benefits derived from the system by the individuals taxed. __ — 
“The construction of a sewer in a street constitutes a substantial 
oar increment to the value of the abutting property—irrespective of how this” 
” property may be at present developed and aside from any actual use which 
may be made of the sewer. This betterment pertains to the land and can 
_ be best expressed on the basis of physical units of frontage or area area W ithin 
a certain distance of the street ora combination of thetwo. | 
“The greatest benefit from the construction of sewers will be conferred 
on those who actually connect with and use the system to a of their 
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daily wastes. In amount this benefit i is logically by the quantity 
of sewage discharged from the individual premises and most equitably 
expressed by rental charges based on the meter records of the water depart- 
er: or, where private water supplies are in use, on estimates or gaugings 
of the amount of liquid contributed to the. sewer. alt _therefore 


construction of a sewer in any =F and —s and — the benefit 


Two comments on the 0 foregoing ab abstract are of interest: 


(a) An assessment, rate per unit, once established, cannot be altered. _ Iti is 


a8 a directive to be applied toa specific project. and not as a general permanent 


vi (b) The relation between the benefit from the contention of sewers to 
5 those actually connected and to the municipality as a whole should not in all 
| cases be subject to the interpretation of the term “greatest benefit” as including 
all or practically all of the benefit. In so far as the sewers contribute to general — 
sanitation and public health, the public benefit is a factor. A Oo 


il. 1. Kinds of Service Provided by a 1 Sanitary District. —The kinds of service — 
provided by a a district depend u upon local conditions, st such as whether separate 
or combined sewers are used or whether there is a large amount of drainage 
through infiltration. — Differences in the extent of the service are also pertinent 
especially as to v whether or not collecting sewers are included. As basis for 
discussing the problem of a fair distribution of the costs of. sewage disposal 
service or of providing for the use of the facilities, the following items seem to 


Be 
(a) The cost of. structures and equipment - needed to provide capacity y for: 
‘The initial population, future g growth, storm w water in combined sewer sy ystems, 
illegal storm-water connections, ground- water infiltration, and industrial 


wastes; 
— (0) The cost of operation, maintenance, and administration needed to to 


permit the use of the sewage . disposal facilities; and 

c) The cost or value of the service given to protect the public h a 

for the ¢ general welfare of the locality. 


The problem i is to find a fair an and reasonable relation between the costs on 
the one hand, and the service and use on the other. Some of the foregoing 
cost items provide service mainly to area and others mainly to users - Thus 
rf the costs of providing ser service for storm water, ‘ground water, and future growth 

“May, in some cases, apply to, and | depend on, area more than on direct use. 
For particular district, a careful computation of the costs of providing these 
- several services and uses may permit the formulation of an equitable schedule 


of charges. Thus the total annual cost of facilities to serve only the present 
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y method of raising the required income must include a part through the iia 
> ae 1 tax 1 part by t on land and a part by rental. I mY 
a _ general tax levy, a part by assessment on land and a part by rental. In i 
. what proportion the total shall be thus divided remains to be determined.” — 
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SANITARY DISTRICTS 


population as compared to the total annual cost of the usual facilities whieh 


7 ‘include capacity for future growth will” indicate the co cost of providing for 
- uture growth 1. This might properly be charged to area, orat least to presently 
unoccupied area. The cost of handling s storm water and ground water n might 

Bor be estimated and this cost regarded as a measure of service to area. = The: 


remaining operating costs would be a charge on users. _ Pie Gis 


ie The service given to protect the public health and for the general wel in: 

“of the the locality i is broad and ‘general, not subject t to actual computation. 7 his 

" is a real service, however, ‘much of it being to users but some of it to p property. = 


Perhaps a reasonable percentage « of the total annual cost can be allocated to 

oa _ The application differs in accordance with local conditions. _ Thus, a sani- 

tary district comprising a relatively new sparsely settled area will justify and 
permit a a wider application to area than an older well-built-up area. — 7 
12. Rate Schedules—There are many different kinds of rate schedules ins 
250 schedules are a admirably ‘summarized in the Third | Progress 
Report: of the Committee on Sewer t Rental Laws : and Procedure of the Sanitar ary) 


schedules are _ adopted to provide for the collection of the needed. 


— in accordance with Bens general method of financing so as to avoid 
unreasonable discrimination. | 

Even in water "rate schedules, Ww vhich have a much broader basis of experic nee, 
there is a marked difference a as, for instance, ,in the number of brackets, the slide. 
a spread between small and large users, and the amount. allocated tot the 


public through fire protection ‘service. 


at 5) 13. Iilustrative Schedule of Charges. —F or For illustrative purposes only, a com- 


‘putation of a a schedule c of charges, rates or rentals, is presented, as follows Ss: 


Distribution of Units.” 


Residen bis 


Tor il ronting on a sewer line (85% « of Item D: 
Parcels smaller than 7 7, 500 0sq ft (08% 


Not on a sewer line of Item 310, 000 


314,000 314,000 

Fronting on a sewer line of Item 7).... 220, 000 
f 9 Not fronting on a sewer line (307% of Item 7). ,000 
Total area units (Items land 7 2,404,000 404,000 
io’, Numerals in parentheses, thus: 8: (6c), refer to corresponding items in. in the Bibliography (see Appe Appendix 1) a } 


Item 
“15 
Ttem 
— 
rg 
23 
— 
1 (4,800 acres) 2,090,000 
26 
99 
39 
33 
5 
a” 


Minimum number of consumers... oth 12,500 


Minimum Number of Intermediate Users: 
Total, Items 12, 13, and 14....... 25,000 
7 ‘By “Area Units” F ‘ronting: ona Sew er Lit ine— 
Small pareels in in ‘residential zones 690, 0,000 
Large parcels in zones (90,000 units" 


Business and industrial z zones (220, 000 units” 


¢ 
Total, Item Items 16, 17, and 206,500 


By ‘ “Area Units” Not Fronting o on a Sewer Line— roe ae! “ad 
Residential zones (310,000 units at $0. 03)... errs 9 300 4 


00 24 Minimum users. (Item 12), 1 2,500 at $6. 75, 000 
8b Intermediate users (Item 13), 9,700 at $12. 00... 
4 B26 Large users (Item 14), 2,800 at $42.00... . 
Special industrial 
Contracts with adjacent ¢ areas... 
otal revenue by use distribution (Item 24 
30 ‘Total revenue (Items 2 23 and 29)....... 
= 31 Less 10. 4% margin for for delinquent account avery 
32 Total estimated net 000 


34 in use e 
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SANITARY DISTRICTS. 


This computation on is related to ‘several units of a an hypothetical “district given 
in Appendix II. The illustration applies more particularly to a relatively new 
oo settled area for or which the s sanitary district is is to and administer, 
the service of collecting sewage as well as treatment. The total annual cost 
5 of $500,000 is assumed to be allocated to “area” and to “use” under the the assump: 
tion that | 40% is s charged to area (Item 33) and 60% to use (Item 34). An 
ape unit”? used in computing the charge to area is 100 sq ft. The charge to 
“ease” is based on the water consumption assuming a completely “metered: 
= system. ‘Some general similarity to methods of financing in use since 
about 1925 in the W ‘ashington | Suburban Sanitary District will be noted. — al 
ok ‘The foregoing illustration is not fully related to established and sometimes 
convenient practice. Local conditions should be considered | within reasonable 
limits of tolerance; 
«14. Financing Rights under Enabling Acts.— —The enabling act, subject to 
the provisions of the state constitution and related court decisions, establishes 
the rights of a sanitary district as regards methods of financing. It either 
: prescribes a method of financing or instructs the board to use an equitable 
method. — The committee favors a broad general grant of pow er rather thana 
Specific: direction. Thus, the of a district might be by ‘moneys col- 


sewage disposal facilities in the district, allocated by the governing belly « ona 
fair basis and by a fair method. 
16 5. Governmental and Proprietary Aspects——The question has been raised 
as to whether sewage disposal and, in general, the operations of a sanitary 
‘district are governmental or proprietary functions. Throughout the forty- 
eight sta states, court decisions and legal opinions can be cited on both 1 sides and i in 
some cases the e opinions se seem to be conflicting. The ‘committee is not able to 
state an opinion on this important aspect. = | feels, however, that the adminis: 
pring of a sanitary district (because an important objective i is to protect the 
_ public health) is more a governmental than a proprietary function, 
Some cases seem to indicate that the method of financing depends « on 


SS the | governmental or the proprietary aspect applies. This may be 


(a) Under the proprietary powers, the district would be financed — “users 
and not by property. - Under this assumption, it is held that : sewage disposal 
is a business comprising the furnishing of a service (or possibly | a commodity 
in reverse) for which users only st should ‘pay. — Under this s application, ‘sewage 
_ disposal would be regarded as a service similar i in n nature to garbage removal 


(db) Under the governmental —r ers, the district may be financed by any 


equitable allocation of the costs sewerage or sewage disposal 


i Further guidance is desired « on this aspect of the subject assigned. to the: 
committee. Reference large of court decisions and legal opinions 
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SANITARY DISTRICTS 


likely that such local citizens w ill be » willing to serve. Ww ae 
oe In general, the administrative board should not only perform its current 


look forward to increasing standards of sanitation, to improved methods o 
sewage t treatment, and to the future needs of the district. 


SUMMARY 
The request for wide discussion during 1945 denotes ‘that many of the 
foregoing paragraphs should b be considered as tentative. The committee hopes 
that, by the e expression of the « opinions and comments of others, ‘they will be 
able later to — a useful final report. : 


Respectfully : submitted, 
ELson T. 


GEORGE E. BARNES M 
SAMUEL A Chairman 


4 Committee on Organization, Financing, and 


January 18, 1945 
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<< 16. General Statement—The administration of sanitary districts is largely *. 
matter of detail and does not involve such fundamental considerations as 
a enter into a discussion of organization and financing. Experience has amply ae 
demonstrated that the ability, character, and intelligence of the members of ——_ 
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Statistical 


"Area, Distribution in Acres— 
257% 


sew: wage (milion daily). . 


4 
5 


Residences (100, ,000 at 4.5 persons 22,200 


Business ess and commercial properties... ,800 


Tote al occupied properties. 25, 000 


-Vaeant parcels of (tem 17m minus 14)... 10, 
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23 
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25 
a 
Units ror A TypicaL MUNICIPALITY WITH SEPARATE SEWERS 
Quantity 
tem 
35 
Business and industrial, 9%................... 720 
Total (12.5 sq miles).......... 8,000 34 
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‘Deseription Wa ‘Thousand dolla 
@ fate « at. 


‘House connections to the stree 


Intercepting sewers.................. 000 


Sew: age treatment plant. 1,000 1,000 


6,000 
Annual Cost of ‘Operation— 


Administration 


ste ations, and sewage 


Intercepting sewers, pumping | s 
we 
treatment plant 


— Total operating cost (Items 24, 25, and 26).. 


Annual Average Debt Service— oe 


a 
@ 2, 
Collecting s sewers and house connections......... 9 
Intercepting : sew ‘ers, ‘pumping stations, and sew: age Ba bat 


Total debt 


Total annual al oper ating cost once 27 plus Item n 30) s 


«Unit Cost Data, Operation, and Debt Service.— 
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other streets parks (Item 
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Total 31/Item 


Gov ERNMENTAL | oF MUNICIPALITIES 


distr rets, ete ) are agencies created usually under ‘the authority 


of a state law for the purpose of constructing and maintaining sewer systems, 7 
. or sewage treatment or disposal Ww orks, for the benefit: of a particular limited — 
area. In many cases they are or » ceniont to > furnish sev sewerage in areas which are 


not incl 
wu cluded w ithin a1 an incorporated municipality. — Be In some cases such districts 


_ May be orgs nized as a means of joint sewage disposal by several municipalities 
with or without adjacent. unincorporated territory. In ‘other cases sewerage 


districts (usually of the assessment type) may | be organized by a munici- 
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to furnish sewers to an unsewered area within al boundaries. 5. Then 


many variations in type and p purpose. 
general, such districts have many of the municipal 


corporations, although their activities are circumscribed within a special field 
of service. In some jurisdictions they are designated to be municipal corpora- 
=a either by statute or by judicial decision. Usually, however er, they are 
termed ‘‘quasi- municipal” corporations. Some of their functions a] appear to be 


cog ‘governmental in character, in a m manner similar to that of 


There appears to be a wide difference in judicial decisions in the several 


states a as to whether the functions of municipalities (and ‘therefore of sanitary 
districts as qo -municipal corporations) i in relation. to sewers age | or sewage 


disposal are ‘gov ernmental”’ or “proprietary” (also called “ministerial,” or 


” 

“corporate,” or “municipal’’) in character. In governmental functions 

are those which : are performed by the municipal c or quasi- municipal corporation 

for the welfare of the s' state, and include such activities as police and fire pro- 

tection, public health, and public school education. — Proprietary functions, by 

Bey those which are performed primarily for the benefit of local 

inhabitants, and as a rule include such public works (frequently revenue pro- 


: - ducing) a as W ater works, gas and electric plants, public markets, tourist camps, 


wharves, public transportation, and m many others. Historically, most of such 
functions were originally performed by private individuals or corporations in 
the anticipation of profit. hether the municipality. or district derives a 
_ ~profit from such activities | appears i in some jurisdictions to be. a factor in 


= = it as proprietary; in other jurisdictions it _ not appear to have 


ON ot only are the ¢ court decisions « at variance between ween states, but i in some 
decisions there a appears s to be confusion of thought and logic ¥ within the decisions 
themselves. 3 The reason for this variance is apparently due to 1 the fact that 


many of decisions from suits to recov er damages from municipalities 


‘the time w ea ‘the king can do no w rong,” that the sovereign eure can- 
not be sued by a private citizen, provided the function or act was “govert- 
mental” ” in character. As s long as governments w were relatively simple in their 
operations, no gr great confusion in thought occurred; but under today’s condi- 
_ tions, with local, state, and national governments: branching o1 out into many 
fields of activity far outside of the early concepts of governmental functions, 

and furnishing services and materials formerly purvey ed by individuals and 
‘corporations as a business, it has become less possible to discern | clearly the 
distinctions between ‘‘governmental” and ‘ “proprietary” functions. 


_ Some of the confusion must arise from a too 1 narrowly cireumscribed view 
of the purpose and effect of sewerage and sewage treatment, and from a lack 


of knowledge of the public health aspects of these matters. — Courts hav re tended 


to look upon sewerage, in particular, as a matter undertaken and performed for 
purely local benefit only, and of little or no state-wide interest. ane oy 


Nothing could be further from the facts under modern conditions of rapid 
transportation by train, automobile, and airplane, in a nation whose people in 
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normal times travel freely over considerable distances. Sanitary defects in 
“even comparatively : small communities t thus present opportunities for the 
‘infection » with communicable disease of a relativ ely large number of } persons, 
who may in turn transmit such diseases by secondary infections over wide 
areas « &B is only necessary to call attention to the epidemic of amoebic 
dy sentery W hich was spread - to nearly all states in the nation in 1933, as the 
‘result of comparatively 1 minor defects i in sanitation in two hotels i in Chicago, Ill. 
‘This was merely a a ‘spectacular illustration on a large scale that defects in pec oa 
tion in any community are at least of state- -wide, and at times of national, 
concern. The epidemiologi cal records of any state will furnish n many testes 
tions on smaller scales that communicable diseases resulting fr from unsanitary 
‘conditions in one local community are carried to other communities. — 
Lad For this and other reasons, the sewerage of communities is a matter 7 


_ public health concerning the state as a whole. 4 It is not merely a local benefit 

to: a particular” area. It is directly and indirectly a matter of public health 

significance. In fact, phase of sanitation which provides for the 


‘disposal o of human e: eaten is a basic factor in public health. It i is s generally 
courts that those activities w deal with 


* urthermore, some of the confusion on this matter apparently arises a 
r 


the fact that most of the decisions result, from litigation involving suits fo 
damage of some type. © . In many such cases the courts seem to have endeavored 


to escape from m the concept o of the governmental ; nature of sewerage and sewage 
‘disposal, as a phase of public health "practice, >, by declaring sewer erage and 


sewage disposal to be proprietary and local in character, in. order that the 


injured party might recover damages from the governmental agency ‘sued. 


As a ‘commentary on this situation, M. Borchard states (7 


such as sew ers, pr eg etc., we find a greater eosin to hold the city 
liable for negligence. The explanation of this disposition can hardly be 
bs found in the usually ascribed reason that such undertakings are not govern- 

mental, but ministerial in it is found rather in the fact that when 
a public enterprise creates a direct nuisance to private property, the 
become “subordinate ‘or immaterial. "Certainly there is nothing» ‘less 
gov ernmental in providing a city with sewers than in providing it with — 

schools and fire protection, and the courts often find themselves involved in 
Bs the problem of determining whether sewers and drains do not partake = 


means to safeguard the public health, usually regarded as a governmental - 


Furthermore, according to Mr. Borchard (7a): 
“Immunity or liability in these matters upon the view of the 
_ particular courts as to whether the duty is governmental or corporate, , and 
in this respect there is the usual w ide divergence, depending often upon the 
particular line on which the courts in jurisdiction got started. There 
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oo On the other hand, certain duties and functions are well settled as 
_ being corporate and not governmental, including the construction and | 
_ maintenance of municipal water and light plants, * * * the construction 
and repair of sewers, * * * and generally the management of property 
owned by the municipality * * *. Furthermore, any business conducted 
by the municipality for profit i inv olves the exercise of cor porate rather than 
_ governmental f functions. * * * And even if a duty is primarily a govern- 


mental one, it is converted into a corporate one w hen it is being conduc ted 


Again, quoting Mr. “MeQuillin (8a): 


“Tt is well settled that a municipal corporation is not responsible for the 

negligent acts of its employees who are endeavoring to carry out the regu- 
‘ lations of the city to promote the public health and to care for the sick and 
destitute. In the collection and disposition of garbage, undoubtedly the 


city acts for the public health and discharges a governmental function.” | 


Rather « obviously, if Mr. MeQuillin ac admits 1 that i in in the collection at and dis- 
posal of garbage the city § ‘undoubtedly discharges a a governmental function,” 
his argument for the proprietary. nature of sewerage fails, for the reason that 


sewerage e and sew age disposal are much more closely” related to the public 
health than i is garbage collection and disposal. 
Fur ther illustration may be had by reference to the views sof J. A. Tobey (9): 
= “Although the proper and safe disposal of sewage is recognized as an 
2. gryems health measure, ‘it is a well- established rule of law in this country 


F by. the disposal of its sew age, since this i is a cor porate ee not a government 


- ae exception to this rule is in the case of the discharge of sew wage into 


Furthermore, quoting Mr. Tobey (9a): 
_ “By the weight of legal authority, nevertheless, the construction and 


ave. 
institution of-a municipal sewer system is a governmental function, but its — 


operation and upkeep is a proprietary function” 


it rell a city i is liable for nuisances sewage dis disposal, 


7 a although not necessarily on the grounds stated by ‘Mr. Tobey, but to state be that, 
z= construction of sew ers is governmental i in character, w hereas the oper: ‘ation 
of a sewerage system is proprietary in character, appears ‘illogical. The mere 
existence of a sew yer system in the the. streets w ould be of no value Ww hatev er to 
the public health; o only w when the sy system is connected to | premises, and in actual 
‘operation, do actual benefits to the public health occur. 
Another fallacy is implied in Mr. Tobey’ statement—that there is 
liability for nuisance when sewage is discharged into ) tidal waters. The exten 


sive fouling of the beaches along Santa Monica Bay by the inadequately treated 
sewage of Los Angeles, Calif.; the os smells along the - easterly and southerly a 
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sewage at the shore line; and some of the conditions i in other harbors, 
to be mentioned as sammaiea to. dispose of this fallacy effectivel 
_ The committee has not attempted to make an exhaustive sear 


decisions bearing u upon either the matter of gov ernmental” as contr 

y “proprietary” aspects of sew verage sewage disposal, or the status of sanitary 

d districts as ‘“quasi- -municipal” or ‘municipal’ corporations. ¥ It is hoped that ° 
“ engineers from the several states will furnish data as to the leading decisions in 

d their states, statements as to the apparent trend of judicial opinion. 


| Suffice it to sti ite th at the decisions Vv — and only a a few will be quoted 


“appear to hav e been aaa, will serve as the first illustration. “Th both cases, 
ie 
a child had drowned by falling into an open tank at a sew age plant, which was 
‘not protected by fencing, and the parent sued the city for damages. In the 
case of Barker vs. City of Santa Fe (136 Pacifie 2d 480), April 14, 1943, the 
‘New Mexico court ruled that the construction and re repair of sewerage facilities — 
was not a gover nmental function, and that damages could be recovered; but 
in 1 Gatcher vs. City of Farmersville (151 S.W. 2d 565), , May 7, 1941, the ‘Tense 
court ruled that sew - disposal | was a governmental function precluding the 


recovery of damages. 
In the New York case of Hughes vs. Auburn (55 N ‘N. .E. . 389) in 1899, it was 
tuled that in the construction and maintenance of a sewer or drainage system a 
municipal corporation exercises part of the governmental powers of the state. 


In the a Carolina case of Dillon Catfish Drainage District vs. Bank of 
Dillon. (141 S.E. 274), January 9, 1928, it was held that the civil code provided — 
for the creation of comer distric ts to promote agriculture, public health, 


>convenience, etc., which is a ‘ vell- recognized gover ernmental function. 
ito = Inthe Georgia case of City Council of Augusta vs. Cleveland (98 S.E. 34: B)y 


February 13, 1919, the co court said (page 346): 
ay e es the opinion that iii ofa city to maintain its sewera rage 
drainage system in a good working order and sanitary condition is a govern- 
mental function. That such maintenance is connected with, and has 
reference to, the preservation of the public health is so well known and so 
generally recognized that courts will take judici al cognizance thereof.”’ ai 

th the Kansas case of Foster vs. Capital City Gas and Electric Company 


(265 Pacific 81), March th 10, 1928, it was ruled that the construction and main- 
tenance of a sewer system is ordinarily gov ernmental function. 
Pi In a Texas case, City of Gle ladesw ater vs . Eva ans (1168 b 3.W. 2 2d, 486), March 
10, 1938, it was ; stated that a municipality, in building a storm sewer, was as 
engaged in a governmental function. In the North Dakota case of Hamilton 
vs. City of Bismarck (300 N.W. 631), November 5, 1941, it was as also held that, 

in ee ting a sewer system, the city ' was acting in inits governmental capacity. 

% However, in a Pennsylvania case involving | sewer 1 enti als (Petition of City 
of Philadelphia, 16 Atlantic 2d 32, October 28, 1940), it wa as ruled that sewers, 
like water sy stems, are owned and operated by municipalities: in their pro- 


prietary capacity and 1 not in n their gover nmental ce apac ity. * 
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. City of Cleveland E. 2d 57 1), 
12 decided 
is in the nature of a governmental function. 


647), March. 12, 
is a gov 
Courts have used the matter of r revenue as a criterion for: the determination 


of whether a municipal function was or was not proprietary. Mr. McQuillin 


“A municipal government has a double iii first, the private, 

_ proprietary function, and second, the legal duties, from which the munici- 


: Pp pality obtains no revenue or other special benefit in its corporate capacity.” 
In the Louisiana case of Floyes, for Use and Benefit of Floyes, vs. City of 
snare (194 Southern 102), March 4, 1940, it was stated | that the pivotal fact 
in determining the exercise of a governmental function was whether it was 
done purely as a public duty, or for pecuniary profit. Directly contrary “_ 
this was the eine decision in Watkins vs. City of Toccoa (189 S.E. 270 
‘September 18, 1936, where it was held that the manufacture and sale of s: oath ; 
toilets by the city was | & gov ernmental function, and the fact that the city 


derived | a ‘profit thereon did not make it liable for negligence. 
2 Them manner in | which revenue for the m maintenance and operation of sewerage 
- facilities is ; obtained, perhaps, has no real bearing upon the determination of 


_ whether sewerage is a governmental or a proprietary function of a city ora 
7 district. The cost of sewerage may be defrayed either wholly or in part from 
p taxation or from usage charges, without affecting the real nature of the function 
as a phase of public health protection of both the state and the local community. 
‘Whether or not a profit may be realized from such taxes or charges appears 
somewhat ¢ to be academic, for the following reasons: (a) From the nature of 
‘the s service, a true profit « cannot be obtained in actuality; and (6) in practice, 
it is extremely doubtful if electors would permit the 2 imposition of excessive 


charges o or r taxes would result in an apparent profit. vey 
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SOCIETY OF CIV 

Founded November 5, 1852 


AND D THE UNIT HYDROGRAPH 


‘oo co. ‘ 
— | +. 


0. Crark, JUN. Am. Soc. C. The discussions represent much 


4 thought and labor at a time 5 hen most of the writers are actively engaged in 
na entirely different t matters. The expressions which corroborate the paper are 
ty gratifying 1g and greatly y appreciated. _ 7 ‘The divergent views outline some fertile 
wid fields for further investigation of problems which are still largely dependent 
~ | on judgment and should stir stimulating and constructive thinking. Some of 
af “the queries and counterproposals presented lie beyond the ability of the writer 


in the closing discussion. . All are sincerely appreciated. 
. Sw eet pertinently emphasizes | that the labor involved i in unit hydro- 
M4 - graph and flood- routing techniques is excessive for routine river forecasting. 


F answer; some have been answ ered by other discussers; and a few are review ed 


ty. 
‘Redue tion of this labor by the development o of simplified empirical coefficients: 
a is necessary and highly desirable. How ever, ‘the simplification 1 might be based 
upon a large number of hypothetical floods computed by appropriate methods 


as well as upon records of past floods. As Mr. Sweet states regarding the | 
James River, nonuniformity of storm distribution may produce variable stage : 
relationships between the e gages along the river. If records of only a few large 
floods are e available, it is easy to assume that relationships are fixed, only to be 
embarrassed when nature proves otherwise >. Howe ever, unit hydrograph oh study 
of hypothetical adverse distributions of storm rainfall will reveal many of the 


variable possibilities wh Ww hich may serve as guides for good forecasts rather than 


as alibis for poor ones. — ‘The James River, w ith its two halves on either side 
of the Blue Ridge Mountains i in Vi irginia, “aliens opportunities for illuminating 
= to any one w ho believes that nonuniformity of runoff distribution on an 


is not significant, or that g -to-g gage stage ‘relationships | are simple. or 


i Nore.—This paper by C. O. Clark was published in November, 1943, . Proceedings. Diseussion on 
this paper has appeared in Proceedings, as follows: February, 1944, by James 8. Sweet, and Otto H. Meyer; 
March, 1944, by L. K. Sherman, Gordon R. Williams, and Ray K. —— Jr.; May, 1944, by E. F. Brater, | 
LC. Crawford, Robert E. Kennedy, and Victor H. Cochrane; June, 1944, by Franklin F. —_— Sep- — 
tember, 1944, by Howard M. Turner; and December, 1944, by Don Johnstone. 


_ “Hydr. Engr., U. S. Engr. Office, Winchester, Va. 
Received by the Secretary January 12, 1945. 
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Meyer clarifies several points in the ‘paper, and appropriately warns" 
of situations in which K is not a constant. Unit hydrograph theories ai and 
time elements of flow are convenient: oversimplifications of a most 
complex form of open channel flow. are no means universally 
applicable. The paper treats only the elements common to unit hydrograph 
theory and to the theories of flood routing based upon constant values of K 
“Lag routing,” to which both Major Meyer and Mr. ‘Snyder refer, is a very. 
practical a1 and expeditious method of flood ‘routing. Fig. 3 and its supporting 
discussions ns demonstrate the essential ¢ agreement between this method and the 
much more laborious application of the Muskingum technique. The con- 


- struction of the time-area concentration curve of Fig. 8 and its routing modifi- 


—eation conform essentially to “lag routing” procedure. 
Lag routing is an empirical approximation of some solutions of the stor: rage 
equation; it i is not itself a solution of the storage equation for any fixed rela- 
_ tionship between storage capacity and discharge. = ‘ti is an expedient method 
‘for use in a natural state of river development, as in flood forecasting or in 
computation ¢ of past floods. _ Howe ever, the writer w would be unable to wel 
4 mine, with it, the effect of improvements which. change the storage capacity 
along a river, such as extensive levee systems—or particularly those ‘improve-_ 
7 ments which change the relative influence of inflow upon storage, such as a_ 


- Thus, for solving routing g problems wh hich involve > a change in storage- 
discharge the vriter believes that the storage equation will be 
r reliable t than empirical “Tag routing. As Major Meyer v w rites, there 
are several methods other than the Muskingum method for the solution oi 


" Mr. Sher man—the ci creator of the unit hydrograph—re eviews the underlying 

theory of this excellent hydrologic tool with masterliness. Possibly through 
i misinterpretation o: of the bar pattern of i Fig. 8 8 and through t the assumption | that 
curve a is a net rainfall pattern of twelve unit periods | rather than a r repre- 

- sentation of shape of drainage area, Mr. Sherman concludes that the resulting 
derive ed | hydrograph w was not a unit -hydrograph. In Fig. 8, a bar p pattern 
7 representing net rainfall would be of infinitesimal ‘duration at zero time. Th 


Fig. 9, 1: rainfall w ‘as principally confined to that on April 25 and that early on 


April26, 


An apology i is made to Mr. | Sherman and to others who found ¢ difficulty in 
checking the computations in ° Table 1, particularly the conversion of Col. 8 to 
Col. 9. The conversion at a rate of 69, 000 (cu ft per sec)- half-days per in. 
a a is a rounded value a appropriate to a . drainage a area of 1 ,300 s sq miles s, rather than 
to one of 1,335 ‘sq miles. The former area is that for w hich the first derivation 
on this’ tor dam site above the gage, “was prepared. 
Intent upon presenting an unretouched, original, routine solution, , the writer 
failed to notice and properly « explain 1 the small difference in drainage | area. 
The routing calculation illustrated is for 15 hours, as stated in Table 1, r: = 


for 9 hours, stated i in text. Subsequent investigations indicated 
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CLARK ON STORAGE 


there is much room for improvement in the analysis. _ The entire origin: ul 
analysis was prepared in in 4 hours and was not then intended | as a , part of a 
_ published paper. The presentation ¢ of these data was intentional. _ The real 
= of a hydrologic tool is much better judged by the rough-hewn original 
product: with its faults, than by the prettier refinements w en investiga- 
tors: possess and ultimately adopt for use. 
~” ‘Mr. Sherman comments on the size of watersheds to which unitgraphs may 
be ‘applied. 7 He correctly states that size is not a barrier to the adequate 
of nonuniformity of rainfall with respect to time. On the other 
hand, nonuniformity \ with respect to areal distribution of the rainfall creates — 
yariations such as those shown in Fig. 18. E ach of these hydrographs resu ulted 
from are rainfall brief enough to be considered a unit rain. — The relative size of | 
hy drographs : at the three stations (gaging, r respectively, 23%, 55%, | and 100% 


of the Appomattox River above _ Petersburg) and rainfall records indicate 


heavier rainfall on one end or the other in each case. Of the hydrogr: aphs al at 
Petersburg, none would be ‘suitable for derivation ofa a unitgraph by any method 
w hich assumed that the hydrograph resulted from a uniform distribution of 
= In the proposed method, which makes no such assumption, an equally — 


single unit hhydrograph, because, although a would represent fow 
from uniformly y distributed runoff, the record does not. . Subdivision of the 
area an and t the use of three component unitgraphs like those in Table 2 produce 


better er results. The methodology a: as extended by y Mr. Turner, incorporating 
= area and runoff depth i in the time-area curve, gives much better results. 


= he possibilities of storage acting to increase outflow drew varied comment 
Besse Messrs. Meyer, Williams, Linsley, Brater, ‘Snyder, Turner, : and Johnstone, 
ie from ‘ ‘controvert[ing | accepted concepts and experience” and “‘erro-— 
neous” to “familiar to all persons who have worked with reservoirs storing 
appreciable volumes of water under the backwater profile.” The writer 
“agrees heartily with those who su suggest ‘ ‘These conclusions should be examined 
critically. The possibility should not be passed ox over lightly, however, 
failure to recognize the possibility. leads, directly to « overoptimism about the 
- effectiveness of certain types of flood control dams, and to underestimation 
. of the importance of avoiding fixed-pool operations at “power and navigation 
dams ‘and unnecessarily rapid gate closures. 
‘The writer’s comments about the -inflow-controlled storage w were intended 
to be qualitatively explanatory (rather than ‘quantitatively assertive) the 
statement (see heading, “Valley Storage’’); “Valley storage does not always” 
decrease flood peaks excellent discussions reviewing the case 
against | such a view rpoint appear t to require some presentation - of the data 
The writer’s use of a channel in Fig. 1 for the qualitative explana-_ 
tion was ill chosen from a quantitative viewpoint. The uniform channel was” 
J seleeted because of the simple, familiar form of the two backwater nol 
however, a uniform rectangular channel would not have an -value 
in excess of 0.5. T herefore, ‘Professor Brater’ 8 quantitative ‘criticism is well 
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founded. A quantitative demonstration might consider that the ‘rectangular 
channel possesses a non-water- -carrying flood plain se several times as wide as the © 


channel at the upper end, ¢ converging to the channel limits at the lower end. — 
"By ay appropriate. selection of the flood plain width at the upper end of such a 
_ channel (which would be quite similar to the channels of many natural streams 
“carrying relatively , small flows in their extensive overbank areas), the volume 
of stored w vater under the backwater ¢ curve could easily be made to exceed 
720 acre- iy: which is the volume Professor Brater’s analysis indicates a as the 


limiting volume that could be passed without surge. Moy 2 


. _ Channels—the variable cross sections of which effectively converge in the 


direction of flow—might be expected to possess storage capacities influenced © 


more by inflow than by outflow . Such channels might | include: | eed 

© Channels converging in variable depth (such | as low-slope channels 

cai _ approaching a a steeper. channel and thus flowing from a zone of g great 
A ‘Tange in | stage to one of small 1 range in stage and at the e same t time 
possessing flood plains of about equal width, or such as reservoirs 
ae fixed levels at the dam. with pronounced fluctuation at 
their upper ends, or such as reaches ending in tidal estuaries) ; and 


7 (6) Channels conv erging in width but of relatively constant range in 


“was shown in Fig. 6, lies on the Piedmont Plateau. Two gaging stations are 
7 located above the Petersburg gage: Mattoax, ‘Va., with a ‘drainage area of 729 
~ && sq miles, is on the 2.0- -day isochrone of Fig. 6, and Farmville, Va., with a 
"drainage area of 306 sq miles, is on the 4.3-day of Fig. 6. 
F % The typical flood rise at Mattoax in feet is much larger than that at either 


ed i Farmville or Petersburg although flood plain widths are comparable for the 
as three ‘stations. The variable cross section of flood flow can be » considered 
wl | diverging i in the direction of | flow from m Farmville to Mattoax and converging — 
on from Mattoax to Petersburg. a Hydrographs of flow at these three stations 
on 


from rainfall of short, are shown i in Fig. (Note that the 


Matton provable ii indication of of ‘storage 


of the Petersburg hydrograph as compared v wits the Mattoax 
hh the second class’ might be short reaches of many rivers, as most flood 7 7 
plains alternately widen and narrow. Many rivers have long reaches through 

= flood waves flow without appreciable chations 40 It seems a little more ~ 


, ©“The Floods of Mareh 1936,” W ‘ater-Supply Paper No. 800, U. S. Geological Survey, pp. 89 and 90 
two hydrogr aw 
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ne : In the first class may be found certain South Atlantic streams which cross the — 

ce Piedmont Plateau flowing toward the Fall Line. The Fall Line is a cascade a 
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ed tially the flow component from above F’armville.) ‘To occupy such a position q 

the in the hydrographs shown, it appears that the flood wave might have under- sy oe 
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rational to interpret ‘such unreduced passage as the result « of alternate 1 slight 
increases and decreases in peak rather than as a passage of an unvaried peak, 
AU record of aw ave passage having slightly more increase than decrease was 
shown in Fig. , in which the hydrograph. of flow at Danville ap appears to have 
built up slightly ‘and shor tened to become the later portion | of the hydrogr: aph 
at South Boston. ‘The discharge from the intervening drainage area Ww ould” 
E have been insignificant on the twenty-first and twenty-second day on 
Most ; rationalizations about the possibility of stor: age operating to build 


2 
4 
af —72' x6! Pavement 
on 0.1% Slope 
Fia. 19. FROM EXPERIMENTS By C. F. ARD AND M. T. Aua 
shows the experimental hydrographs w vhich C.F. Taaard, Assoc. M. Am. Soc. 
Cz. , and M. T. Augustine. obtained by | sprinkling a a , pavement. a a In this case, ik 
the: change from one flow distribution profile to another did he place rapidly 
~ enough to cause a surge. How fast does the change take place in rivers? yi 
How fast does storage volume react to a change in flow? 
“dd ‘Several discussers mention the assumed instantaneous reaction of storage ™ 
to a change in flow (which is characteristic of most flood-rc routing techniques an 
_ and favorable to reduction of peak flow in te method < discussed wl hen z is” tio 
Ness than 0.5) as the error which leads to th® author’s conclusion. The formula dh 
7 for reaction time presented by Messrs. Ww illiams and Brater is quite favorably | th 
regarded by many, whereas Mr. ‘Snyder a appears to ) regard any value less than 
as inadequate. Several intermediate Viewpoints have some adherents. tin 
=, The time elements indicated by Eq. first discussed by Mr. Williams, & 
be om very long values for a shallow depth at the beginning. of a flood to crest o 
a depth values which are very short in relation to the time elements of mean ay 
a at ‘Preliminary Report on Analysis of Runoff Resulting from Simulated Rainfall on a Paved Plot,’ ~~ 
i pe C. F. F. Izzard and M. T. Augustine, Transactions, Am. Geophysical Union, Vol. 24, 1943, Pt. I, p. 5 500. ‘ .. 
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-flood-wave e travel. _ Thus, in the reach of river to which Mr. Williams —e 
the formula, the channel depth was more than 50 ft at the time of crevassing. © 
For this depth, the formula discussed v would indicate less than 40 min. (N either 
the methods of observing stage, nor those of defining the time of crevassing of 
a levee are sufficiently reliable to distinguish betwe een the time of transit as 
originally | ‘stated and that given by the formula. * This value \ would be in n-— 
significant | in flood routing ¢ over a vensle involving flood-wave travel time of 


15 hours. 


_ The writer does not claim th: at the upgrading 0 of flood crests by storage 
action has been proved, but merely | offers the explanation and the examples of 
the paper in support of the statement, “There is some valley storage | which 

does not. decrease flood peaks.” The upgrading of hydrographs does appear 
possible, however, and should be investigated. A closely analogous situation 
exists in electrical _ circuits, and the upgrading (in height, not discharge) of 


hydraulic 1 Ww: aves in 1 converging channels is a recognized phenomenon i in tidal 
-hydraulies. 
~ Mr. Linsley clarifies several ideas in the paper a his discus sion. ' { 
neds he presents for determining time of concentration in watersheds . 
without records shows better cor relation with data in the paper than the 
author’s formula. | To aid in further study by others, additional data are 
_ presented in the same manner in Table 6. The majority of v alues were 
‘developed from records available ‘prior to 1938. variation in the co con-— 
stants: suggests judicious local application. Attention is again: directed to the 
Wy two sentences ‘under the heading ,* Determination of the Instantaneous: 


Hydrograph for Streams Without Flow Records,” and to Mr. Snyder’s closing 


twoparagraphs, 


Messrs. Linsley and Johnstone properly question the writer’s use of uniform 
velocity in time contours along and 


A 


may be of tw OW ‘making a right, as as the 
likewise neglected principle of using greater storage values, _K, on flows from 
more remote : areas instead of a single average value is a logical refinement of 

- Opposite effect. In general, the use of a varied velocity of channel flow and 
a single value of K -tends to increase the fault originally recognized by 1 the 
in item 6(b) of the “Summary. 

essrs. Brater Johnstone \ voice an interest in a more detailed descrip- 

tion, of the time-area concentration curve. As Mr. Snyder states, the curve 
is an ol¢ old device*-* for calculating flow > It has fallen into disuse because its 
2: assumes es the reductive influence of storage to be negligible. A storage 
correction factor was included in the paper to revitalize this useful ok The 
time-a concentration curve is derived by marking ti time along t the 
"watercourses, the time to th 

Ne w York, N. Y., 1922, p. 94. t 


a ‘Calculation of Flood Discharge by the Use of a Time-Contour Plan,’’ by Cecil N. Ross, Trans- 
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. _ "Sewerage and Sewage Disposal,’’ by L. Metcalf and H. P. Eddy, McGraw-Hill Book Co., Inc 
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a concentration.’ ” For ‘the latter term, a standard definition such as that quoted 
by Professor Johnstone i is intended LL ‘Fig. 4 illustrates why the terminus of this 


- time is not. at the end of the hydrograph (theoretically infinite in length) but 
Ww ellup onthe recessional leg, = 


TABLE 6. —ADDITIONAL CoMPARISON ‘18a AND 18b 


(sq_ 

i 7 iles (hours) miles) | 


4 North Branch, Potomac River, at Bloomington, Md. 100 
Savage River at Bloomington, Md 
Georges Creek at Franklin, Md. ee 5 
Wills Creek near Cumberland, Ma 
_ Patterson Creek near Headsville, W. Va 
South Branch, Potomac River, near Petersburg, W. Va. 
South Fork, South Branch, Potomac River, near 


Cacapon River at Y ellow ‘Spring, Ww. Va 


0.089 
0.018 
0.071 
0.044 
0.052 
0.068 
0.16 
0.04 
O17 
0.026 
0.11 


343 


= 


Back Creek near Jones Springs, W. Va............. 
Conococheague Creek at Fairview, Md 
North River near Burketown, Va. 
Middle River near Grottoes, Va.. 
South River at Waynesboro, Va 
North Fork, Shenandoah River, near Strasburg, Va. 
Passage Creek at Buckton, Va 
4 Cedar Creek near Winchester, Va. 
_ Antietam Creek near Sharpsburg, Md.. 
Monocacy River near Frederick, Md 
_ Goose Creek near Leesburg, Va 
Seneca Creek at Dawsonville, Md 
Jackson River at Falling Spring, Va. 
James River at Lick Run, Va ass 
Cowpasture River near Clifton Forge, Va...........| lk 
Calfpasture River at Goshen, Va 
North River at Rockbridge Baths, Va 
North River near Lexington, Va 
North River near Buena Vista, Va. 
Hardware River near Scottsville, Va a. 
Rivanna River at Palmyra, Va 
Willis River at Flanagan Mills, Va. . 
Appomattox River at Farmville, Va. 
Appomattox River at Mattoax, Va 
Appomattox River near Petersburg, Va............. 
-Meherrin River near Lawrenceville, Va 
River at Emporia, Va... 
Roanoke River at Roanoke, Va... 
Blackwater River near Union Hall, Va 
Pig River near Toshes, Va 
Snow Creek at Sago, Va 
Goose Creek near Huddleston, Va. 
Otter River near Evington, Va. 
Falling River near Brookneal, 
Dan River near Francisco, N. C 
Mayo River near Price, N. C 
Smith River at Bassetts, Va 
Smith River at Martinsville, Va...._ 
_ Smith River at Spray, N. C 
Sandy River near Danville, Va 
Banister River at Halifax, Va 


—) 
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une 
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Time of concentration, not storage facto. 


[r. Craw ford | presents ‘interesting analyses of channel-storage, correlation 


stages : at the discharge point and ‘upstream from it, ‘and shows the superi-_ 

_ ority of correlations which give weight to upstream values of inflow and stage. 

He calls attention to the potential possibilities of determining from discharge 
records other storage relations which would be too | costly by cross- -sectional 
surveys. The writer has seen area-capacity curves" se small ‘dam sites 90 
a developed which would require very detailed field surveys for equal a accuracy. 
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Much of this fallacy may originate in ofa finite 
base length. is storage theories and unit hydrograph theories are to be cor- 
related, the Presi of finite base length i in a unit hhydrograph must give way 
to a concept of an infinite length, a characteristic possessed by the writer’s 
instantaneous hy drographs and their r derivatives. } Nevertheless, the practical 
differences, as shown by Fig. 15, are within the range o of : enenlly acceptable 
_ hydrograph computations. This is probably the reason that the practical 
simplicity of finite constant base length has so long outw reighed the theoretical 
ies Mr. Cochrane has very thoroughly examined the procedure for deriving 
unitgraphs" utilizing storage concepts, and presents so some explanations w hich 
| should be helpful to those who were mystified by the Ww riter’s presentation. 
; Ina simple manner, he clarifies the twofold effect of storage: ,: First, the greater 
: effect of storage in creating time lags between flows from different zones, or the 


time-area diagram; and, second, the further modification of Thy drograph shape 


* ~The order i in w hich these steps have been utilized is the reverse of that i in 
— ‘Fig. 8, but the result is the same regardless | of order. | The order in Fig. 16 
‘gives: a clear ‘picture of the mechanical effects, and would be necessary if one 
wished to use variable values of K for different : zones. — At the sacrifice of this. 
possibility, the order of Fig. 8 involves less work. However, a unit hydrograph © 
for an entire area may produce very disappointing results if applied to a storm 
characterized by great nonuniformity of runoff distribution. _ Mr. Cochrane 
a a very practical application of storage principles which considers 
both nonuniformity of runoff distribution and some variation in storage effect. . 
__ Another method of accounting for nonuniformity, suggested by Mr. Turner, 
is to extend the e time-area conc concentration curve to express not ‘only shape but 
“also different runoff. distribution—that is, multiply the zone areas: by the 
of runoff in each zone before mening on the routing procedure of 


and is. 


z Both procedures abandon the unit hydrograph and require complete re- 
calculation of each flood with attendant increased labor. c Other less logical, 
_ but extremely useful techniques to account rapidly for nonuniformity include 

P adjustments of chronological sequences of runoff used w ith the unit hydrograph. 

J Thus, concentration of volume on more remote parts of a simple watershed can 
be simulated in some degree by using an n effective runoff volume for a later 
time than actually wa vas the case. Extreme intensity. of rainfall within a unit: 
_ Period can be simulated byi increasing t the computational amount in that period - 

: and decreasing the values on either side of the pe period. The ‘most flexible 
_ technique yet utilized by the writer is the development of unit hydrographs for 

very small zones, which are routed downstream to gaging stations and tabu-_ 
lated in the manner of Table 2. This procedure is followed for as small areas’ 
- desired to express nonuniformity to the necessary degree. These sub- 


‘Unitgraphs ca can always be added together for use on | areas | as large as justified 
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CLARK ON STORAGE, Discussions 
‘a the uniformity of any storm and can be used separately if necessary. = Al 
though this procedure we was discussed under the heading, “Unit Hydrographs 
Large ‘Drainage Areas,” watershed is large from the viewpoint of 
Refined techniques are necessary and desirable to more suitably 
correct: procedures and to confirm academic viewpoints as_ to the effect of 
certain qualitativ e variables. _ How ever, a computation is only | a guide to 
; judgment. After a hydrogr aph i is computed with a unitgraph or other appro- | 
priate tool, judgment demands some modification o: of the answer if the premises 
- _ at thé methods and the facts of the occurrence are not in accord. At times a 
wide line, a small scale, or a little freehand “artistic license’ may have as much 
merit as as hours of computing. : In their proper place, these too are proper tools. 
a ‘Mr. § Snyder presents sowed of the practical al aspects of hydrograph ealcula- 
tion and the mathematical basis for a very simple and useful method of routing 
4 over long reaches of river where storage-discharge relationships are such as to 
produce essentially constant time elements. He calls attention to elements of — 
logical and mathematical fallacy which may appear small to some investigators, | 
but which must be faced squarely and solved before the hydrograph computa- 
7 tions can be considered accurate enough | to prove. or - disprove some of the 
Mr. § Snyder ‘states, the Muskingum method of flood routing | loses its 


a usefulness when applied to too long reaches. = However, this may also | occur” 


because the storage-discharge relationship for a long reach of river is much ; 
less closely approximated by a straight line than is the relationship for short 
reaches. The storage item in flood-routing formulas such as Eq. 5 deals with — 
. the volume stored in a reach at any instant of time, whereas the storage capacity 
of channels s having nearly ¢ constant, elements of flow time (n = 1.0 in Eq. 3) 
refers more > closely to the volume under the high- -water profile. Although this 
latter volume is the total storage ¥ volume utilized i in the reach of river during 
flood passage, it is more than is used at any instant during that time. The 
7 ‘difference becomes a material one in long reaches. Nevertheless, if storage 
solutions | are confined to the length of reach for which E q. 5 would be appli- 
cable, Eq. 6b is a correct solution. — 7 The Muskingum method is a correct 
solution of Eq. 66 only to the e extent t that a a finite value of T is a satisfactory — 
‘approximation of dt. Mr. Aumes expresses ses the op opinion that values of 7. used 
the preparation of f Figs. 2 and 3 were tooshort. at 
a . inadequacies ¢ of methods of mathematical solution receive too little 
recognition. _ Most engineers find it necessary to adopt arithmetic approxima- 
‘tions, such as Eqs. 1 and 6c, for the solution of calculus 1s equations, such as” 
Eqs. 2a and 6b.“ he practises are so common that they escape much com- 
_ ment. _ The simplicity of the ‘Muskingum m method encourages these procedures, 
since it depends first on the substitution of a finite value of time T for the 
infinitesimal time dt. it. Then, too, many convenient, quick "solutions evolve 
from various lengths « of time Tw Mah sabe an one or more of the constants either 7 


a Oor0. Values of T as large as K or larger conveniently obscure some of the 


44 “Graphical Integration of the Flood- Wave E quations,’’ by Harold A. Thomas, Am. 
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effects of the inadequacy of of the theory recognized by t the writer ‘in his dis- 
‘cussion of Fig. 2, and have been accepted by some as a basic requirement. 
The differences are as given by Mr. Snyder. In an effort to prove whether 
such a basic requirement, added anything to the over-all accuracy of hydro- 
graph solutions, the writer has for several years assumed that the shortest 
practical time element, ¥, was the most appropriate. on Tn | general, practical 
ee have been no more — and i in some cases souenoetly better, 


To be the may need some provision (1) to cover 
“the flow conditions in which the lag of the crest discharge i is gre: sater than the 
lag between the centers of mass of the inflow and outflow, and (2) to cover . 
the conditions in which the recession rate of outflow exceeds the recession 

rate of inflow. Flow conditions requiring this modification are evident in Fig. 


18, in which the hydrographs of the Appomattox River at Mattoax are shown n 
to have a materially smaller rate of recession than does that portion of the hydro- 
graph at Petersburg | which encompasses the flow volume from above ‘Mattoax. 
. Inspection on of Fig. 17 and Eq. 25 ‘suggests the applicability of a negative 
value of T.. Thew titer finds the mathematics of a solution evasive, however, 
and doubts that this is a practical ; solution. : Nevertheless , it is explanatory of 
a relationship betw een K and time of concentration (T.), which still baffles 
him. However, it is s assumed for the ] - purpose of the next paragraph that nega- 
tivevaluesof T,areapplicable, 


a he f failure of of the w riter to solve this possible adaptation of lag routing is 
the source of his confusion in statements under the heading, ‘“T he Hydrogr aph 
as an Index of Storage,” about which Mr. ‘Snyder r rightfully comments. — The 
values of K applicable to smaller watersheds | appear to be as large as 0.5 of 
“the time of concentration (7); yet those for larger watersheds have a much 
‘smaller ratio of T, and occasionally a smaller dimensional quantity. Obviously 
“this j is possible if 7, could have a negative value. For the three gages on the 
Appomattox River—Farmville, Mattoax, and Petersburg, each site being 
ised downstream from the preceding site—the values of K which the writer ‘regards 4 
«BB as most appropriate are, respectively, 9 hours, 12 hours, and 9 hours, whereas" 
ttle the most acceptable times of concentration, T. ' are 32 hours, 114 how, and 
hours. (These values for Petersburg differ from values i in Table 1, which 
were the very first computations ‘made for that stream. The differences be- 


tween the values of to 144 and 168 is wand large, 


application, 24 hours or less, would but 4 which i is close to 
limiting error fixed by rainfall data. that Gifierence 


flow conditions, a factor usually ignored, but by no means absent. ap: beeroxqou 
The J Appomattox River presents the most perplexing flow conditions to _ 
Which the writer has yet applied the unit hydrograph theory and the a 
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of concentration for comparable drainage area that he has ever encoun- _ 
tered, as well as a unique shape. Although the illustrated hydrographs for 
4 the Smith River s at Bassetts and the James River at Lick Run are certainly i 
the more common type, exceptions, like the illustrated hydrographs of — 
_ mattox River and Meherrin River, are necessary to prove the rule. Judging 
from records of rounded and tree-shaped watersheds typical of the glaciated 4 
ah alluvial streams of the Midwest, it would appear illogical t that a stream 
could rise for a longer period than it fell, or that streams could have 
: hydrographs like those presented for the Appomattox x and Meherrin r rivers, 


Mr. ‘Sherman’s doubt that any unit hydrographs presented i in the ‘paper 


Mr. . Turner provides a a well-balanced onintanut the paper, confirming 
statements and expressing appropriate reserve and doubt where his experience 
7 and understanding of conventional concepts and procedures do not justify 
7 presented conclusions. His ct comparison of results obtained from the outlined 
~unitgraph procedure and from his similar procedure™ encourages confidence 
principles applied. extension of the method to inelude ing 


refinement whose merit would depend on the extent of the use for dun. “He 
expresses several views in sharp conflict with those of the paper. The pro- 
-priety of defining a unitgraph to include nonsurface runoff is basicly criticized. 
Granting that some investigators have presented hydrographs attributed to 
_ surface runoff, the writer questions w shether such a claim can be substantiated 
or disproved, since there is no feasible method of distinction, no acceptable 
_ specific definition of either kind of runoff, or any impervious boundary 
5 _ betw een the flow channels utilized by e each. _ Practise i in the application of the 
- terms seems to justify distinctions satisfactory to the user and determinable 
only to the extent of their compatibility with a _ preferred theory. ‘The w rriter 
‘indicated that separate unitgraphs | could be used for surface and subsurface 


flow, if warranted, and that there was nothing fixed about the 70:30 ratio. 


4 that definable quality ‘which the use of the unit hydrograph theory hee 
In his demonstration that there “* * * * should * * * be no no possibility of 
A multiple solutions” of conventional unit hydrograph derivations from m rainfall 
and stream flow data only, Pr Professor Johnstone adequately presents the reasons — 
why an investigator who uses all the data available can secure: : (a) A distress 
ingly large number of different solutions if he assumes he is able to determine 
runoff from rainfall; or (b) no ‘solution if he admits inability to determine runofi 
- from rainfall. = In the storage concepts presented in the paper, there are » addi- 
tional considerations which eliminate the necessity for dependence on rainfall 
and runoff determinations. . Until runoff determinations become more reliable, 
however, the 1 real accuracy of ber method must remain partly : 
me ‘‘The Flood Hydrograph,” by Howard M. Turner and Allen J. Burdoin, Journal, Boston Soc. . 
Civ. Engrs., V Vol. XXV Til, N No. 3, 1941, p. 232. vr 
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RY WAVES 1 N NATURAL ; 


M. A. CHURCHILL. 


A. Assoc. M E.4«—The Society i is indlebted 
ser to the author for the seemutation and interpretation of a wealth of observed, 
He not theoretical, data on translatory waves in natural channels. In an attempt 


TO- | to make the p paper of still g greater value, the writer presents herewith some sup- 


to. TABLE. 11.— SPEED or TRAVEL 4 AS OBSERVED, In 1943, on 


ary bas (a) Wave TRAVEL, (b) Wave TRAVEL, 
; Fauurne Face; Ristne Face; Nov 12 

time at | Time - time at ime ime - 
face — Mile ‘midpoint | in | served |served| in served 
tio. of fall =| hours] rate? of fall j|hours} rate? | hours} rate? 
Nanees Ferry 38.2 4.67) 3 83| 4. 9.00| 1.80 

nfall Highway 70 Bridge 5A m. 12:10 

mine * Average flow was 7,900 cu ft per sec. In miles per send _ ¢ Theoretical wave travel rate, 3.88 miles — 
wall perhr. 4 Theoretical wave travelrate,3.95 miles perhr, 

nfall porting data « on observed waves in natural channels and in reservoirs, together 
able, with. corresponding data on observed rates of water travel. 
> in _,, Nore.—This paper by J. H. Wilkinson, M. Am. Soc. C. E., was published in June, 1944, Proceedings. _ 
1as appeared in Proceedings, as follows: January, 1945, by Harold A. Weeeel, 


ia Hydr. Engr., TVA, Knoxville, Tenn. 


_ “@ Received by the Secretary December 28, 1944. 
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CH URCHILL TR ANSLATORY 1 WAVES 


a Wave Travel in Lower er Holston River.—In connection with a determination of a 


- time required for w: ater to travel 46.6 miles from Cherokee Dam on the Fg 

Holston River to a point near 


Below Cherokee Dam | the head of Fort Loudoun a 
Reservoir above Knoxville, 
a Tenn., the velocity of @ nega- 


of. translation was 


crease 7 ,740 cu ft p per sec. 
negative wave was ob- 
served at four stations on the 
‘Holston 1 River,below Cherokee 
as ‘shown by Fig. 10(a) 
2 (stages) and ‘Fig. 100) 

charges). _ Stages were ob- 

= at all four stations, but 

"discharges are known only at 
oa the station immediately below. 


also observed. The negative 

wave was produced by reduc- 

tion for a period of 24hrinthe 
discharge of Cherokee Dam on 

November 7, 1943, from 7,530 

cu ft per sec to 3,840 cu ft per 
12 sec, and ‘then an abrupt in- 


3 


le 5. 


Chloride Concentration, in Parts per Million 


“y lio Cherokee Dam and at ‘Straw 
Plains. (Stages below Cher- 
’ 8 okee and at Straw Plains 
7” are from continuous recorder 
| charts. . Stages at Ni ances 
é a Ferry and Highway | 7 
5S Bridge’ are from staff gage 
readings at intervals of 2h 
less. Discharges at Nances | 
9$ erry at Highwa ay 70 
3 
on this negative wave. a 
computed wave velocities, when the computed velocities are based on only 


two, relatively short, rated ‘iver reaches, is “more luck than sense.’ 
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4 author neenmeaily based his comparisons of observed and computed rates on 

; similar data and so it is rather remarkable that his results on Clinch River 

show such good agreement with the theory. ae pane 


he 
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Fic. CoNncENTRATIONS AT OBSERVATION STATIONS ON THE LOWER 
_-Housron ‘NOVEMBER | 12 AND 13, 1943 
s dong? sal) To 


4 


> 


W ater Travel in ‘Lower Holston River he ‘Tate of water travel, as d dis- 
tinguished from wave travel, through the same reach of river was dived 
_ between November 12 and 14, 1943, during which time the release from 
Dam was nearly. the same as both befor re and after the negative wave 


“TET ALI 
(©) STRAW PLAINS, MILE 17.0 | | 
(d) HIGHWAY 70 BRIDGE, MILE 5.4 | 
4 


CH {URCHILL ON _TRANSLATORY W/ Discussions 


generated on on November s 1943. (The release from Cherokee Dam on No- 
- vember 12: increased gradually from 7, 770 cu ft per sec at 1:00 a.m. to 8,070 
eu ft per sec oc at midnight. ) To deteendine the rate of water travel, samples: 


observed. These samples were analyzed chemically for chloride concentra- 
tion, with results as shown in Fig. 11. The presence in the river of an industrial 
waste of varying chloride concentration made possible the use of the chloride 
concentrations for identifying a , particular m mass of water at the suc successive 
down nstream ‘stations. ions. Table 11(c) gives pertinent data on the rate of water. 
— Ratio of Wave Velocit, y to Water Velocity, Lower Holston River—A com- 
parison of the rate of wave travel with the rate of water travel in the lower 
Holston River indicates that. the Tate of wave travel is slightly more ‘than 
twice the rate of w ater travel, for the particular rates of discharge observed. | 
This ratio tends to ‘support the author’s contention (see heading, ™ Translatory. 
Ww aves on the Clinch River: Velocity | ofa Translatory Wave’’) that “The w “| 


velocity is about half the wave velocity at ordinary stages * * *.” In fact, 


> data presented in this discussion on the relation of wave to water velocities 
* probably better proof of this statement (as it applies to the general case) 
then the data presented by the author, since the author used water velocities 
{ at only four points in the 64.2-mile reach of the Clinch River between the 
ave and W ater Travel, Ww atauga and South H ‘olston Rivers. —Observations 
have been made of the time required for water and for or Waves ( of translation to 


_ travel through a reach of 40 miles on the Watauga oni South Holston rivers. 


The waves were produced by variations in release from a power dam at the 
head of the reach. . Discharges varied from 200 to 1,300 cu ft per sec. The 
average time for a series of waves was determined to be 16 hr. Water travel 
the same reach required 91 hr. 


Comparison of wave- -travel and water-travel rate es s shows that, for ‘the 


particular | rates of discharge e existing, the waves traveled ed 5. 7 7 times e as fast as. 


the water, 
‘This ratio of wave to water velocity is rather high for a —r river 
channel, particularly considering the rather wide range of flow between crests 
and troughs of the various waves. As the author states, however, the ratio 
will be higher for low discharges, and during these observations the av onl 
Ol rate of flow was about on one third of the mean annual rate. _ The existence of 2 
relatively large number of long, deep pools in the lower Watauga and § South 
= Fork Holston rivers could explain the high ratio. In any case, the wide differ- 
ence in the relation of wave to water velocities found here, as compared with 
the relation | found on | the lower Holston River, indicates the possibilities of 
‘error in trying tou use wave NS as a a basis for computing water velocities 
and Water Travel, Fort Loudoun Reservoir. —Observations of wave” 
travel through Fort Loudoun Reservoir (42.8 miles) showed an over-all w rave 
velocity of 21.4 miles per hr. _ The theoretical Tate of wave travel | was co com 


puted by ‘Eq. 2b. The mean depth of the reservoir was found to vary e mail 


were collected at 2-hr intervals from the same four points at which the wave 
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February, 1945 CHURCHILL ON TRANSLATORY WAVES ad 
tially as a straight line “#" om 14 ft at the upper gage to 48 ft at the dam, when 


_ the pool elevation is 811.5. The wave velocity, except for the influence of the 


water velocity, ¢ can be expressed as vg Ya = ; or: 


0.176 (8) 
J, 000015007. 14 + 0.00015 0. 1504 L 


al Jo 
\ 206,000 
7, 460 se sec = 2.08 


‘and t = 0.176 


14+ 0. 0001504 A 


pore in the upper end of the reservoir. 


; zero. _ The theoretical wave e velocity i is thus computed to ‘0 be 2081 6 miles per hr 


a as as compared with the observed velocity of 21.4 miles p per hr. aha 
_ Measurements of actual water travel by chemical identification a 

indicate that a period of 190 hr is required through the 42. 8 miles of reservoir. 

| This rate of travel is 0.22 mile per hr. . Had the moving water occupied aol 
‘full depth of the reservoir throughout the entire length of the pool, the the water 
velocity would have been even slower than the observed value. QQ 
the The w vave velocity was 95 times S$ as fast as the water velocity, stall 
This ratio explains \ why the author’s conclusion is true that the water velocity 
has little influence on the wave velocity in a reservoir. “yey 

_ Basic Concepts of the Relation of Wave to Water Velocity. = an attempt to. 
get a clear understanding of exactly how a translatory wave could run ahead 


of the water which ey formed it, the w _ formulated t the following 


_ Assume a smooth, , straight natural channel in which the flow i is gradually 


Femnos from a steady flow of Q, toa higher flow SS ~ Under su such conditions, 
4 profile of the channel and water ‘surface w ould “appear as: shown in ae 12. 
Since the hydraulic radius is” is gres ater the 
“higher flow Q2 , the velocity V2 is greater than 
assuming, that is less than bankfull 
discharge. Since this is the case, water moving - 
ath mean velocity V2 tends to to overtake water 
moving at mean velocity Va . This action pro- CHANNEL oF aa 


: duces a result t similar 1 to an n express | train cr: ash- 

ing into the | rear end of a slow freight. In that ‘ease the rear cars of the 

- freight are telescoped byt the ‘express. In the river, the slow er moving Ww ater is 
“ 

a. ‘telescoped” by that moving at velocity Vo. ¥ ‘The quantity of w rater per second 
80 “teleseoped” is expressed mathematically by (V2 Aa, | in which A, is 
the cross-sectional area . corresponding toV;. This quantity of w ater is actually 


Squeezed out of its original ee in the river channel and is forced into that 
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WAVES 


portion of the » channel. Tepresented by Az — A, 
thereby added each se second in front of that at the velocity Vs 


7 ‘Siti: the rate of advance of the wave front is — , relative to V:, 


The actual rate of advance of the face. of is 


‘the v water i provided the ont has an initial flow i in it, ‘since ‘the 
second term is always: positive. e. By simple algebra, this expression for the 
ave velocity ea ean | be. ‘shown to be equal to the familiar expression, Eq. 1 1b, 
“except for the constant used to convert feet per second to miles per how. | 
W ave-Height Reduction. —The author explains: what happens to a transla- 
tory wave produced by a release of short duration, when its height i is reduced 
(heading, “Characteristics s of Translatory Waves”), without explaining the 
Teasons: for these happenings. ‘The author states that “* * * the ‘spreading | 
fillets at. the end- of-rise and at the beginning-of-fall gradually : shorten this 

[erest] period as the progresses downstream.” It seems. to the wv writer 
4 that, in addition to any so- ealled “ channel-storage” effects, the basic reasons: 
why the crest of a wave of short duration is gradually | lowered are as follows: 


- As the wave moves downstream, it pushes the water already i in the channel 


of the channel represented by As — Ai. However, this pushed water is not. 

spread uniformly over the area At — Ay but rather it tends to resist movement. 
and stay near the | lower edge of this, area, , where the channel width is the 
narrow est. This nonuniformity tends to. weal the beginning-of-rise point out 
farther and farther ahead of the crest of the Tising wave. . As the nose runs 
farther ahead, the rising face of the y wave ‘is thereby flattened; only a little 
water is added to the rising face of the wave, at and near the crest. a — | 
_ Although little water is added on to the forw ard end of the crest, a con- 

quantity i is removed from the rear. Since water the crest stage 
‘moves at a velocity higher than that at at the _ er stage to the rear, a qu: antity 
of water equal to (V2 — J V3) A; is removed from the rear face of the wave e each 
| second. pe _ This ecnstant drain on the rear face of the wave eats it away and ‘thus 
shortens the crest length until finally the original crest is all gone. Then the 
remaining crest is still further flattened until finally no wave exists. ae 
An additional reducing influence i is the fact that, due to the increased slope 
of the r rising face, the water r composing it tends to run faster than the falling 


face where the slope is flatter. Thus the wave tends to stretch out in length 


at the base, and consequently i is lowered at the crest. 1 


ahead of it, as explained previously herein, and forces this water into that : area 
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“AMPLIFIED SLOPE DEFLECTION. 


al Zz JOHN E. GOLDBERG, 


E. GoLpBE Assoc. M. _ Am. Soc, E amplification o of the 


slope- deflection method is . presented i in this concise and excellent paper. “ The 
fundamental geometry of ‘slope deflection is combined with the author’ s con- 
cept of the traverse angle as applied ‘to members under flexure. _ The writer 4 
is gratified to perceive that, as this and other interesting papers indicate, the 


UL 


value of f slope deflection as a pow —_ method of analysis is steadily becoming 


_ Several of Mr. Stew art’s hideaii deserve critical consideration. One 
of these is the sentence (see he: uding, ‘ ‘Demonstration Problems’ | “In con- 
ventional slope deflection the number of simultaneous equations required | to 
‘solve a given problem is s definite.’ _ This standardization is a distinct advantage 


of the slope-deflection method foun which both the experienced | and the i inex- 


perienced ‘engineer will benefit since any engineer can see at once the proper 
direct slope-deflection ‘approach toa given problem. On the other hand, the 


slope- deflection | method is so so extremely flexible that it allows much leeway, 


when desired, in the solution of special problems. 


_,, Consider, for example, the problem of Fi ig. 4. _ Some years ago the writer 
“developed : an extremely simple analysis for problems of this ty This 


method has the advantages of being readily applicable to any number of 


stories and of being extremely rapid since the basic equations ‘converge very 
‘quickly even in extreme cases. . Applying the method to Mr. Stewart’s al . 


-réblem, two equations would be set up immediately 
(6 X 24+) 1) = 0.5 (20 


— 


On (6X2+ 1) = - 0.5 5 (20 X 10) + 0.5 5 (20 20) + + O4.....(198) 


—This paper by Ralph W. Stewart, M. Am. Soc. C. E., was published in September, 1944, Pro-- 
ca ree on this paper has appeared in Proceedings, as follows: November, 1944, by Leon | 


bs "ie. Consolidated Vultee Aircraft Corp., San Diego, Calif. go, Calif. 
Received by the Secretary November 30, 1944. 
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DEFLECTION Discussions 

The solution « of Eqs. 19 Ox = 9.524 and @p = 23. 81) agrees exactly with the 
author’s solution. Thus, the method is shown to be ‘simple and direct. Of 

. course, this i is merely proof of the fact that the most useful me method nd genera ie 


is the one with which the engineer is most familiar. 
Mr. Stew art goes on to mention the same frame with the lower oil 
pinned at one end and states that six simultaneous equations are required by 
the conventional slope-deflection method. This i is absolutely true; but, when 
these equations are set up (and setting them up is in itself very simple), thre 
of the equations may be eliminated in a very direct manner. 7 Note also that 
a conventional strain- ~energy analysis can be made in terms of three redund: ants, 
Mr. Stewart’s point concerning ; sign conventions (see heading. “Definitions 
of Symbols and Other Data: Signs,’ ’ and text following Eqs. 7) is well taken. 
In his original | presentation!? of slope ‘deflection, G. A. Maney, M. Am. Sor. 
CE, adv: anced the sign convention which the author advocates. In sub- 
sequent work,131415.16 both Professor Maney and the writer continued to use 


convention, finding it simple and conv enient. 


The author’ statement, ‘discussing Eqs. 7, that * * the final ‘moment 


ima 
the end of a beam for which the moments are not ‘dominated by joint 


is less than the fixed- end moment” is likely to be misinterprete! 
iy the reader. _ The : statement w vould be clarified greatly ‘if the author hat 


written that “the moment is less, algebraically, th than the fixed- end moment.” 


CAMILLO Weiss,” M. Am. Soc. C. E 7¢_The four “ properties of an 1 elastic 


curve traverse presented in this paper are principles derived in part from 
“moment areas” and are ‘important enough to be compared the well- 


known p principles s of the moment-area method. 


The eo concept of the A- A-angle i is helpful in visualizing t the effect of variation: 


2 cross sections or end r restraints. _ ~The author is justified i in stating that in 


cases the number of simultaneous be reduced by recog: 
nizing the significance of the A- angles. 


The method would seem to deserve a more designation 


“amplified slope deflection.’ “Although 1 related to slope deflection, the un- 


-= ms of the author’s method consist e exclusively of slopes and the method 

is based on its own distinctive principles, 

‘The author should be commended clear and able presentation | of 

thm principles and | for | his skilful derivation and interesting of 


: 12 Engineering Studies No. 1, by GA. Maney, Univ. of Minnesota, Minneapolis, 1915. —_ — te 


‘5 a, 13‘*Vertical-Load Analysis of Rigid Building Frames Made Practicable,’’ by John E. Goldberg, 
Engineering News-Record, November 12, 1931, pp. 770-772, 


14 Simplified Methods for the Analysis of Multiple Joint Rigid Frames,” by George A. Maney and 
= E. Goldberg, Northwestern University Bulletin, Vol. XXXIII, No. 7, 1932. 


15**Wind Stresses by Slope Deflection and Converging Approximations,” John 


_ Transactions, Am. Soc. C. E., Vol. 99 (1934), pp. 962-985. 


es. ‘‘Natural Period of Vibration of Pas Frames,” by John E. Goldberg, Journal ¢ of the pe 


Concrete Institute, September, 1939, p. 81.00 


Designer, Bethlehem Steel Co., Fabricated Steel Constr., Eng. Dept., Bethlehem, 
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DISCUSSIONS 


tire | ANALYSIS OF STATICALLY INDETERMINATE 
STRUCTURES USING F REDUCED EQUATIONS 


tions Diseussion 


4 


: 


@ 
JosER SORKIN 
joint JosEr Sorkin, 5M. Am. Soc. C. The thesis presented in this paper is 
pretel offered (see “Synopsis’ ”) as one that “* useful in the analysis 


or had of certain types of statically in indeterminate structures * 7” 7 paper 


nent.” belongs to the category | of a number of others published since 1930, all at- 
_ tempting to simplify the application of the classical elastic theory and nearly 
all, unfortunately, ending i in utter failure. The: only ‘notable exception is the 
“Analysis of Continuous Frames by Distributing Fixed-End Moments 6 by 
Hardy Cross, M. Am. Soc. C. E., in which a radically new approach | substitutes: 
simple for the s of simultaneous equations. This statement 
implies no reflection upon the sincerity of purpose of the other w Titers. 
recog, a The simple fact is that most statically indeterminate structures, as en-— 
countered in actual practice, do not lend themselves to solution by specific 
afl simple mathematical expressions. Vv ‘arious solutions may be devised for any 
* age assumed academic ¢ cases, but it is impossible to assume enough ¢ cases to meet all, 
or even most, of the conditions encountered in actual design That is is. 


primarily: the r Treason \ why efforts | to develop a “cure-all formula’ are futile. 


shat in 


“methods and that it does not apply to all structures. Its primary purpose 

“appears to be simplification of approximations s of members preliminary to final 

design. design of a statically indeterminate structure by any method 
yoldbers _Tequires that the size and the properties of the members 1 must I bes assumed before 

onal "proceeding w with the stress analy: sis. . Iti is not clear w hy the same e principle does 

ni 

rower - not t apply to the method proposed by the author, regardless of his statement to — 
mated ; ‘the contr: ary. Irrespective of the — used for determining the reac reactions, 


; the shape of the member originally assumed must finally be checked for stresses. 


American 


a Notr.—This paper by Lee H. Johnson, Jr., was published in November, 1944, Proceedings. 
Engr. of Design, Howard, Needles, Tammen & Bergendoff, ‘Kansas Mo. 
Received by the Secretary January 2, 1945. 
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Any part of the structure feud to be either deficient or extravagant must be 

: modified correspondingly with suitable readjustment of the reactions and the 
‘stresses. — In prognosticating the proportions of such a . structure, nothing can 

ever take the place of a designer's judgment based on past experience with 
_ similar structures. . The question appears to | be whether the methods proposed 
are in fact ‘simpler | than other methods available. 

Ino far as two-hinged arches are concerned, certain rudimentary ‘ate have 

in existence for many years rs. Thus, for e preliminary d a 


parabolic arch, it is customary to assume OF 


5PL.. 


in which k L = distance from support. ~ Incidentally, Eq. 32 is the a 


: Eq. 16a, since j = 2k. This formula can be employed with only a a slight err error 

for the preliminary. proportioning of a circular arch of usual proportions. Os 

That no greater precision is necessary | is indicated by the use of Table 6 as 


{ Satieail in Table 7. Values for the maximum influence ordinate, as obtained 
‘TAB ABLE 7 .—ComParison OF INFLUENCE ORDINATES IN TABLE 6 wiTH Compt- 


| Arch _ ew | | 


Reduced” Actual Eq. 32 Reduced” | Eq. 32 


7 by the use of E Eq. 32, are w vell within ‘the required accuracy for preliminary 
design, compare favorably with the values obtained by the use of the author's 
reduced equations, and are more simple to compute. Preliminary determi- 


nations within | 3% are fairly close. Variation of details may account for 
At its best, the method appears to be directly applicable only to solid con- 
crete ribs of such shapes as the author (see paragraph ereyeg Eg. 16b) has 
‘selected ‘ ‘primarily because of their mathematical simplicity.” alt tis is clear ‘that 
for steel ribs composed of shapes and plates, where changes i in sections are 
abrupt, the mathematics would become so involved as to make the method too 
cumbersome for practical use. The author claims that the method is als. 
d applicable to arch trusses and to three-span continuous trusses; it is yet to be 
demonstrated, how ever, just how such application would be useful i in any but 
“very simple cases, where any relation to a practical design problem would = 
ior: It may be p iailinsd to state that the actual time required for the analysis 
of any usual structure is a very small part of the total time required to prepare 


__ 17*The Theory and Practice of Modern Framed Structures,” by J. B. Johnson, C. W. Bryan, and F. E- 
Turneaure, 9th Ed., John Wiley & Sons, Inc., New York, N. ¥., 1910-1916. 
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designs, uns, and specifications Hence e, the betw een one or two, or 
designer as a his own approach to a recog- 
nizes “short cuts” as his experience broadens. The quickest est design method for 
any structure is the method that the designer knows and understands best. I 
Of late entirely too much emphasis seems to have been placed on methods 
and not ‘enough on the physical meaning of the mechanics of 
Thus, as an example, the significance of Eq. 8a, for a horizontal reaction of - 
two- -hinged arch, is simply that the horizontal reaction is equal t to the horizontal 
movement of an arch end due to superimposed loads (if the end were free to 
move) divided by the simils ar movement caused by a unit load applied at that — 
end. “ Such a a definition is easily understood and the problem then is cone of 


‘meaning of Eq. 8a can easily te by. moments, ordinates, 
differential elements—the actual significance of the meaning being lost in the 
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"ENERGY IN FLOW EXPANSIONS: 
Wwe th 

| of 


oO 


pone a statistical date, 1 Kalinske 

paper which is based on a st aggering number of photographs. se 

Ww Thenever w ater decelerates in a conduit under conditions o of steady flow, 

there is a conversion of kinetic to potential energy. This conversion is not 

made without loss in excess of that induced by the surface of the conduit. 

“i The losses _ which accompany | this conversion of energy have been studied 

_ intensively, and the methods of study can be classified | broadly into two 


(1). Quantitative experiments on a particular: installation or model from , 


oo _ which over-all losses incident to the conversion can be deduced ; and 


ee 


(2) Mathematical analy ses of general and often abstract motions of an 


deductions and can contribute much directly to the fundamental understanding . 
* the mechanism of energy ¢ conversion. 7 Both approaches to the study of the 
conversion of kinetic to potential e energy have yielded significant contributions. 
_ The engineering approach is usually to seek a particular solution of 4 | 
problem—leaving to others the related problems of correlating, generalizing, 
and rationalizing the data into a general pattern. — The mathematical approach — 
usually begins with an idealized problem, leads to a general hypothesis, and 
leaves to others the task of applying these generalizations to. particular cases. 


Too often these two approaches to the problem— springing f from opposite s sides 


of a chasm of unknow ns—do not meet to form a well-integrated structure of 


af theory and experiment. investigation reported by Professor Kalinske was 
‘directed toward studying the internal mechanism of conversion in 


Received by the Secretary January 6, 1945. 
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duits with straight axial expansions and toward forging anoth 


integrated rational structure. 

na The purpose of this discussion is (a) to classify the elementary phenomena 
of conversion of kinetic to potential energy in hydraulic engineering, and (6) 
to present selected bibliographic data which may suggest the course of further 
laboratory studies and illumine the weakness or ee the mathematical 


‘Table 5 classifies the fundamental problems of the hydraulic conversion of 7 
kinetic to potential energy and the signific: ant geometrical variables for each 
eae. - Common to all the problems listed are such physical variables as rate 
of discharge, roughness of conduit, state of motion, head loss, and velocity 
distribution. ‘There are but two classes of flow. WwW ithin the reach of 
the conduit under consideration the flow i As either fully confined by the walls” 
of the conduit (Table 5(b)), or not fully - confined | (Tables 5(c) and 5(d)), in 
which case there is a free water surface in the reach. The interplay of un- 
balanced forces and accelerations of liquid elements is maintained i in the first 
‘case by changing pressures along fixed boundaries, and in the second case by 
an aperiodic pulsating free water surface or by ‘@ complex sone © of air- entrained 


A conduit w which fully confines the flow within a reach under consideration 


may induce ‘deceleration of the liquid by the flow ¢ or by 


changing its direction. There are four primary « 


Mas 


* A straight axial expansion which is represented ir in its its simplest by 
the lower leg of a venturi meter, by a straight, conical draft tube, and by similar 


-— 3 ks axial-to- radial expansion of which the spreading draft tube is an 


example 


8 An rove in which the flow can be decelerated by y changing the direction 
a flow only, or by a combination of changes in direction and changes in average 


velocity. Elbows may, therefore, be subdivided further into— 


‘Prismatic elbows in which the cross section and shape of the conduit 


remainunchanged; 


ales ‘Report of the Committee on Hydrodynamics,’”’ by Hugh L. Dry rden, F rancis D. Murnaghan, and ~ 
_ H. Bateman, Bulletin No. 84, National Research Council, W oer ~ C., 1932, pp. 469-481 (bibli- 


of the Committee on Hydrody namics,” by Hugh Francis D. Murnaghan, and 
H. Bateman, Bulletin No. 84, National Research Council, Washington, D. C., 1932, pp. 481-492 (bib- > 
__ _18“A Study of the Hydrodynamics of Spreading Draft Tubes,” by Andreas Luksch, thesis presented 7 
to the State Univ. of Iowa, Iowa City, in June, 1935, in partial fulfilment of the requirements for the 
= of Doctor of Philosophy (bibliography, 52 entries). 
“Hydraulic Tests of Small Diffusers,” 
Bulletin o. 13, Studies in Engineering, Gale. a Iowa, Iowa City, 1938. 


18**Flow Characteristics in Elbow Draft-Tubes,” by C. Mockmore, ‘Transactions, Am. Soc. C. E., 
103 (1938), pp. 402-464. 


at. Bibliography of Important on Draft 1913,” National of Standards, 
U.S. Dept. of Commerce, Appendix I to report on investigation of draft ‘tubes conducted for the Tennessee 
‘alley Authority, June 30, 1935 (bibliography, 78 entries), 


8“ Plow of Water Around Bends in Pipes,” by David L. Yarnell and Floyd A. Nagler, ‘Tran 
Am. Soe. C. E., Vol. 100 (1935), pp. 1018-1048. 
"Flow of Water Around 180-Degree Bends,” by David L. Yarnell and Sherman M. Woodward 
- Bulletin No. 626, U. S. D. A., 1936 (bibliography, 19 entries). 
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vous 


ad Bisymmetrical ell elbows i in n which the oTOSS section or shape is symmetrical 


~ with reference to two » planes—one in the plane « of the elbow axis and 
the other normal to that plane and bisecting the straight | approach 


Expanding elbows by elbow draft tubes. 


that ‘the water is discharged into a lock chamber 
during | the filling stage, the manifold may also be classified as an example of a 
in a conduit with a (flow not fully confined 
by the walls of the conduit) has two major subdivisions :*- 24,25,26,27 0 
4, if the velocity in a reach is supercritical (velocity head greater than 
half the hydraulic mean depth), the water either passes through an hydraulic 
jump, or continues in some form of jet. The hydraulic jump has been studied 
intensiv ely in rectangular conduits, circular conduits, trapezoidal conduits, 
and conduits of other sections. — Surface waves and plunging jets have been 
studied in open . channel flow and submerged w weirs s chiefly. 
Ii. If the approach v elocity in an n expanding conduit is subcritical (velocity 
‘ton less than half the hydraulic mean depth), the hy draulie "problem of the 
transition becomes primarily a matter of backwater curves. This topic has. 
received extensive experimental and mathematical study. Practical problems 
. in this | category ¢ can be handled quite satisfactorily by combinations of labora-_ 


tory stu study and conduit analysis. be 


Essentials in the study of energy conv ersion in flowing are Newton’s. 


aws of motion, a complete physical conce pt or hy pothesis, exact definitions of 


terms, ration: al induction and deduction, and objective comparisons and | ex 


periments. The engineer has usually “appre oached the study with the objective 


of determining « over-all losses of energy y rather than w ith the purpose o of studying 
the internal mechanism of energy conversion. me Many guide posts to engineering 
judgment have been established on the basis of so-called one-dimensional 
ts or asia — solutions of problems i in hydraulics. » Five illustrations of ‘ ‘one-dimensional” _ 


a _ approaches to to problems involving conversion of kinetic to potential energy in 
ait 


f flowing water are sketched in n the following paragraphs: ‘ali 


A. Newton’s second law of expresses the relationship the 


**Lock Manifold Experiments,” wm Edward Soucek and E. W. Zelnick, Proc eedings, . Am. Soc. 
October, 1944, pp. 1255-1274 


of Open Channels,” B. A. Bakhmeteff, McGraw-Hill Book Co., Inc., ‘New 


. _4“*The Hydraulic Design of Flume and Siphon Transitions,” by Julian Hinds, Transactions, Am. Soc. - 
% “Hydraulic Jump in Trapezoidal Channels,” by | Cc. J. . Posey and 8. 8. Hsing, Engineering News- 
Record, December 22, 1938, pp. 797-798. 


"928-929 Determining the Energy Lost in the Hydraulic Jump,” by J. C. Stevens, ‘ibid., June 4, 1925, » PP. 


9, 


. " “Hydraulics of Steady Flow in Open Channels,” by Sherman M. W oodward and Chesley J. Posey, 
ohn Wiley Cone, Inc., New York, N. 
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Discussions 


ON FLOW EXPANSIONS 


Precisely consistent relationships can be expressed between the work done on 


a ‘particle and the change i in its kinetic er energy, or between the impulse of a 
4 force acting 0 on a particle and the change i in nits | linear momentum. _ Consistently 
applied, the principles of force and. acceleration, work and energy, and impulse 


and momentum lead to identical conclusions. Evidence that this fact is. 


J sometimes misunderstood is the inference that Bernoulli’s theorem (essentially 


ow ork and energy) ‘‘applies” in some cases s and that the momentum principle | 
_ “applies” i in others. — When all factors are considered i in any given problem, 
identical conclusions thew the logical application. of the e principle of work and 
energy or the principle of impulse and momentum. Students of hyd draulics 
go sometimes do not realize fully that apparent ‘ . ‘differences” are but logical 
which follow from differences in viewing the ‘complete physical 


picture oF concept 
| . Be The principle of work and energy a) applied to the flow of a liquid in a 


stream tube leads directly to Bernoulli’ s theorem. x ‘If real water flows in a 
- real conduit, the walls of the conduit exert a drag on ‘the flowing liquid w hich 
opposes the motion. In mathematical derivations and algebraic manipula-_ 
tions, if the physical pictures are incomplete or inconsistently simple, “other 
factors” are sometimes introduced presumably on rather 
es as hy potheses whose validity must eventually be tested by experiments. 
4 C. ‘The principle of f impulse and momentum applied to declerating flow of 
water in a level « open conduit leads to an expression for depths of flow adjacent 
et the hydraulic: jump. ‘Simple assumptions lead to an expression for h head 
loss in terms of alternate or sequent depths of flow v Ww hich is in substantial 
agreement with laboratory. and field observations. Suffice it to state that, if 
the principle of impulse and momentum is logically applied to accelerating . 
“flows: of water in a level open channel, ‘the valid conclusion must be reached — 
that the water surface promptly | becomes level as a consequence of the basic _ 
definition of a fluid. On any horizontal plane ‘the fluid cannot indefinitely 
-Tesist shearing stresses which must be ‘set up by unbalanced pressures parallel 


ul D. The principle of impulse and momentum applied to Pen eae in 
a 90° expansion | leads to differences i in pressure above anc and below the e expansion — 
which depend on the assumed p pressure distribution on the face of the expanding | 
section. a Reasonable assumptions of this pressure distribution | lead to. 
_ puted head losses which are in substantial agreement with laboratory obser 


_ E. The principle of impulse and momentum applied to “nonresilient im- 
pact” of two particles moving in the same direction with different initial 
velocities leads to an expression for head loss which is identical with the ex- 
pression for head loss in a 90° expansion ina closed ed conduit. — In the ¢ over- -all 
7 analysis of head loss in a sudden expansion, 1, nothing ne need be nema a as to the 

detailed mechanism of conversion. . Particles at the boundary of the jet 
approaching the expansion come in contact with ‘slower moving particles te 


mixing—with reduced The term ‘ ‘inelastic impact,” 
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4 confined liquid is ‘‘elastic” in the sense that it has a well-defined bulk modulus Ss 

of elasticity. Individual globules of an unconfined liquid are not resilient, 
“however, i in the sense that mechanical en energy can be stored in an ‘enced 
globule by virtue of its elastic deformation. Henee, the term nonresilient 
impact ¢ seems more descriptive. ; This nonresilient | impact of liquid 1 masses may 

‘7 ' be viewed as a physical basis for anslysing the conversion of kinetic to potential 
: energy, and it may throw added light on the author’s second conclusion that 


( maximum total turbulence energy is a small part of the total energy 


change taking place. 


Professor Kalinske has outlined briefly certain concepts relating to tur- 

| ulence an and he has referred to selected literature. _ The engineer will doubtless _ 
- feel a void between these concepts and their application to even an elementary 

"engineering situation. _ If statistical methods of studying the phenomena of 
turbulence and the nonrecoverable ‘components of kinetic energy in a stream 

are to be applied to a given conduit, one must know the primary pattern of 


“flow through the conduit, the intensity of turbulence at ph gy point in in ate 


expense of refinement in techniques. If the ‘data are are 
obtained from model tests with the intention of translating results from one 
conduit to o another, it follows that the two conduits must be > geometrically — 


r similar, the intensity of turbulence at corresponding points in ‘the two ) conduits 


_ The practical significance of the statistical approach to fluid turbulence i in 
the study of the conversion of of kinetic to potential energy has not yet reached : 
far beyond the laboratory and the « computer’ S desk. _ The engineer certainly 


itely cannot look to the analyst of turbulence problems for the design of a diffuser, 

allel a draft tube, or the hydraulic system for a lock. 4 However, if the basic theory 

of fluid turbulence prompts the experimenter and the to look 3 more 

yw in critically into the effect of baffles and the internal similarity of flow patterns, | 

nsion the mathematical approaches—as yet of limited direct. application—may 

ding “eventually make their effects felt in more important problems of hydraulic 
bser- _ The study described in this s paper was initiated by the wr writer, in 1 1936 at 

oe “the invitation of J. C. Stevens, President, Am. Soc. C. E., chairman of the 

t Society’ s Committee on Hydraulic Research , and progress reports have been 

nitial abstracted periodically in Civil Engineering. Between 1936 and 1939 a bib- 

liographic study was “completed by. Andreas Luksch, graduate assistant in 

er-all” hydraulies at the State University of Iowa, and selected tests were conducted 

the by Mr. Luksch; by Edward R. Van Driest, Assoc. M. Am. Soc. C. E., and 

jet Arthur Luecker, Jun. Am. Soe. E., , research assistants in mechanics and 

les hydraulics; and by Hesi- -Hou Chang, graduate student. In 1939 and subse- 

pact, ently the work on the project was directed by Professor Kalinske. 
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for. prior or development of the photographic hich played| 
4 such a significant réle in the studies reported by Professor Kalinske is due to 
_ Edgar E E. Ambrosius , John C. , Reed, and Henry F. Irving** who pursued a_ 
pre ogram m of studies in hydraulics laboratory of the University 
Z Illinois at U Urbana whic hich was under the direction of M. L. Enger, _M. Am. 


“Study of the Flow « of Ww ater Through a Glass Pipe, by Edgar E. ‘Ambrosius, John C. Reed, and 
Henry F. Irving, preprinted papers and program, Aeronautic and Hydr. Divs., A.S.M.E., Summer Meeting, 
June 19, 20, sana 21, 1934, George Reproduction C Co., San a Calif., 
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AMERICAN ENGINEERS 


STABILITY AND STIFFNESS OF CELLULAR 


ite DEAN P. Testa « wg 
Dean P. ‘Am. Soc. C. little has been. 
written or published on the rational design: procedure and the behavior ex- 
“hibited by ‘this type of cellular « cofferdam construction. Therefore, “assump- 


~The p paper constitutes : a contribution to 

_ the engineering profession on this subject . As the author states, tests should 

be made, but results should be used with discretion, s since they have often 

following comments discussion are confined to cofferdams in which 


The inherent economic advz antages in using - circular cells with ‘connecting 

ares are due primarily to simplicity of construction and re-use on subsequent 

- stages and projects. Furthermore, cells can be filled independently without 

_ interruption, and possible interlock failure w ill be localized in one cell. As a a 

general rule, if a cofferdam does not exceed 50 ft in height and is carefully — 

‘designed for iieaseul interlock, allowing for saturation of fill within the. ai 

resulting construction will be satisfactory and safe. 

Designing on ar an “equivalent width’ basis a length of unity) was 
found to be sufficiently accurate after a comparison was made with the —— 

‘obtained from using a typical cellasa unit. This expedient obviously simplified | 


this 


—This paper by Karl Terzaghi was in » 1944, Proceedings. Discussion on 


this paper has appeared in Proceedings, as follows: December, 1944, by G. G. Greulich, and Raymond re 


Pennoye 


4 Structural Dept. TVA, Knoxville, Tenn. 
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2260 TSAGARIS ON CELLULAR COFFERDAMS 


_ ini of the great height of the Kentucky cofferdam (100 ft ii rock) 
and of the magnitude of this undertaking as a whole, very careful consideration 
was given to the following factors: (a) Saturation of fill material within 1 the 
cells; (0) interlock tension; and (c) internal resistance to shear. Iny view of 
this concern, the w riter cannot agree wi ith the statements | made on the subject 


2 internal shear, Ww which read (see second paragraph following Eq. 8): ; eee 


before this failure [by overturning] occurs, the octierdam i is likely to fil from 
another ca cause which has never received any attention” 4 and (see second pa gle 
graph preceding Eq. 28a) ““* * a potential source of failure by shear along 
"vertical planes (Eq. 19b) w hich has been ignored i in the past * * 2 ou 
(a) Saturation of Fill Material Within the Cells —On all co cofferdams built 
by TVA, a mixture of sand and gravel with little silt content was used to obtain 
as much free drainage as possible. Filling was done by the hydraulic pri a 
except f for a few cells built in the ary, such as tie-ins. - "When the first ae ; 


dams Ww ere constructed, the need for dependable information on saturation — 
was realized. a For ‘or purposes 0 of f study, se several rows of pipes were provided within | ; 


7 


the cells, and observations were re taken, after unw atering, x, on the level of water 
resulting from bot both a rising and receding river. Ree wow 
ae general, results « disclosed a rather flat gradient across - the cells, _ 
ginning from a point somewhat below the outside river elevation. No | 
plete drainage was evidenced in spite of the customary weep holes on the dry 
a side. In any case, much depends on the care exercised in the field and the 


extent of the provisions made to ameliorate saturation. wide 
saturation line w Ww ith an approximate slope of on starting from the 


“4 on 1 were where ‘special ‘care was to cells” w drained. 
Hence, all calculations were made accordingly on subsequent jobs, using 110 lb _ 


below the saturation line. — As a result, all cofferdams required a berm on the 
dry: side to ) satisfy the aforementioned Tequirements. 7 This. does not necessarily 
‘mean that cells cannot be made wide enough to avoid berms , provided they 
check in every | other respect and are economical. A horizontal line, at an 
elevation so chosen as to represent the average expected condition of saturation, 
_ should serve just as well and at the same time simplify computations. 
— (b) Interlock Tension. —This is one of the most important features in 
cellular cofferdam design and should be kept w w ithin safe limits. ae The p practice 
of computing pressures on a straight-line -line variation all the way to rock is 
fallacious since cells almost invariably e: exhibit a characteristic bulge — 
one quarter to one third of the way up from the base. ‘Driving apparently — 


forces the steel sheet piling into rock sufficiently to prevent any appreciable 


per cu ft for material above the saturation line, and 65 lb per cu ft for ene 


movement at this point. hen berms are. used a similar condition exists, 
except that: the overburden o or berms ser serve to restrain the the 


& tension is computed | should be considered 1 maximum at three quarters of 


the way way down from the top in cases without berms. if used, how- 
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-TSAGARIS ON CELLULAR COFFERDAMS 


pore opening of interlocks have resulted n near the bottom, 

but repairs ‘can be made easily, especially i in cases where no berms are used. 
- All’ possible conditions should be investigated to determine maximum interlock : 
. tension, particularly during filling o operations w vhen the river elevation is low, 

because this condition most often controls, 
‘The a author does not mention the amount of "tension contributed by the a0 
“connecting ares to the existing main cell tension. This been found all 


on the main cell through the tees or w yes. The neglect cannot be pion 
and computation as a free force does not tonsider the passive re resistance of the 
fill within the connecting ares, which resists a large percentage of the tension | 
before it can n find its way into the main cell. 


by TV A engineers with cofferdam design, and investign- 
tions in this connection have been thorough. No rigorous mathematical 


actual pressure the i of the piling, and ‘the: resistance by 
interlocks are so variable. | _ Very ably, Professor ‘Terzaghi explains the im- 
plications involved i in computing actual shear and evaluating the resistance to 


this shear along a a selected critical plane. 
nie Relatively small horizontal deflections are capable of mobilizing high pres- 


sures s within the fill, especially ag: against the inside surface of the piling, with 
large resistances. This tendency will continue as the lateral 
increases 3 provided the cofferdam is still safe against interlock. failure, sliding, 
and overturning. Furthermore, if slippage occurs between sand and gravel 


"particles, the new conliguestion of material may conceivably offer a greater 
to than results of computations indicate. Readjustment 


“nots so highly will come into play to the horizontal 
“movement due to slippage stops. This behavior transcends: the realm | of 


hypothesis; otherwise the paradox ¢ of nonfailure for cases with factors of safety 
less than 1 unity cannot be explained. . Since all cofferdams built by TVA were _ 


temporary, a factor of safety of 1.25 for shear, sliding, and interlock tension 
was considered quite adequate. 


2 resistance, followed the general by the 

- that it was Tecognized that active pressure w ithin the material was a 
“by the effect of external horizontal forces tending to | compact. the fill confined 
the cells. A friction coefficient of 0. 0.55 was used. To the foregoing, the 

resistance aj against interlock slippage was added, using a coefficient of friction 

“of 0.3, which gave the total effective resistance to shear. ‘Under the most — 
7 unfavorable conditions, the smallest factor of safety was 1.2. It is felt, how- 

ever, that the cocficienta of friction for both internal shea ‘and slipping of 
‘- interlock were quite conservative. 
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ON CELLULAR COFFERDAMS 


a Recommendations.- —To summarize, the writer would like to offer 


er- leaf t ty with connecting ares should be used; 
‘Selected free dr draining ‘fill: material, such as sand and gravel, will be 


found best; 
3. Computations should be based on satur: ated fill below line with 
‘approximate slope c of two horizontal to one vertical (for simplicity : a horizonts l 
line may be adopted at an elevation which represents average conditions) ; 
= * Shear resistance should be considered and computed on a logical and 


effects. of the characteristic bulge j in.sheet pode he 
7 6. Grouting of the foundation rock will result in a nr and st afer t coffer- 
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DISCUSSIONS 
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BY E. GRINT ER, , GEOR 


E. Grivrer,® M. Am Soc. the column is probably the 


most critical structural member, it is well for the structural engineer to r review, 
at intervals, the background of procedures for column design. Design of 
columns has long been associated with the use of a “column formula” and, 


therefore, the author’s excellent. paper again stimulates the sion to re- 


4 even w the were all per in one laboratory. This 
‘is true because of the virtual impossibility of producing a centered load or a_ 
perfec tly if it must be from Parts. 


structures w here soil movements, temperature movements, and 
“a he choice of a column formula, therefore, is something of a psy ee 
or philosophical mi: atter ather thea a strictly Scientific decision. A straight: 
line, a broken line, a parabola, or a reversed curve can each be . chosen to 
represent test results, probably with equal satisfaction and accuracy, within: 
the. common Ww vorking range of slenderness ratios. If this” is true, then the 
choice of a formula should be influenced by another consideration: What © 
column formula will produce the best results in terms of design when placed in _ 


the hands of the: ordinary designer who i is not always” a highly competent 


9 oar | 


There are two mistakes that can be made in the use of a column formula. 
~The formula itself may be misinterpreted ; and this mistake w ill become more 


“probable as the formula becomes more complex. For this reason either the: 
Rreiiger formula or the Aarflot formula mentioned by the author could become 7 


a source of errors in design. - ‘The second mistake is in extending the use of the 


2 Nore.- —This paper by William R. Osgood was published in December, 1944, Proceedings, 
Vice-Pres. and Dean, Graduate School, Illinois Inst. of Technology, Chicago, Ill. A 


“Received by th Secretary J 4, 1945. > 
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GRINTER ON COLUMN FORMULAS 


rmula beyond its specified range. This can be omssomuiall by siti { 
necessary limitations on the formula or by selecting a formula of the kind 
saint by the author for which extrapolation bey: ‘ond test results will 
reasonably safe results, 
Perhaps it is ‘unfortunate, but, nevertheless, much design work is s actually 
_ performed by persons s of limited competency i in ‘structural design. — For them 
the straight- line formula or the Rankine formula is the reasonable limit of | 
Be complexity. — Since no greater accuracy can be obtained by use of a | 
more complex ead. the advantage of adding terms i is not very clear. The 
best form of printed specification and one which would meet all. ‘objections , 
would be a complete statement of the design formula to cover each range of | 
slenderness ratio. The only ‘mistake that has been common | 
writing has been ‘to cover only the normal range and simply to neglect entirely 9 


the | upper and lower ranges with resulting confusion in the ae of mary £ 

"persons who must, and do, perform design work. 


GEORGE Wi INTER, 7 M. AM. Soc. C. E. Until a decade or two ago prac- 


simply to approximate, as ‘possible, the 
results of a greater or smaller number of column tests. Little attention was | 
= to the physical significance of these formulas, to. correct representation | 

of end fixities, manufacturing peculiarities, and sometimes even to physical | 

characteristics of the material. ‘Since then the | development has been along 

more rational lines . Methods were developed | to account for column ‘strength 
te considering, i in 1 a rigorous : analytical way, properties of material, end fixity, 

initial shape of the member, manner of | aad application, ete. \ Two types of 


= (a) T The one designated by the author as double-modulus theory which 
takes caaanh of the deviation from the straight line of the stress-strain curves, 
particularly at higher stresses. This method proved rather cumbersome in 


f ‘its is application » and, as as the author states, is of only limited practical consequence. 
i (0) The ‘other takes account of the fact that | the strength of a column de- 
pends to a great extent on its initial shape and the manner of load application. 
_ Initial deviations from straightness, eccentricities of loading, and other factors 
which might be termed imperfections were found to decrease the strength of 
5 columns considerably, particularly i in the range of low and medium slendernes: 
“ratios. _The secant formula is of this latter ne as is the method proposed b by 


gr great advantage of methods (b) consists in fact that any term 
- entering such a column formula has a definite physical significance and the 


af entire formula i is s based on ri igorous elastic theory rather than — represent 


“Rational Design of Steel Columns,”’ by D. H. Young, Transactions, Am. Soc. C. E., ‘ai 101 
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ane ing | designer such a formula can be adapted to alr almost any | set of particular: 
circumstances, such as eccentricities of loading, imperfections of shaping, | etc., 
which is entirely impossible with the purely empirical type | of f formula. a! Too 
much reliance on the latter a actually 1 may lead to unsafe design, inasmuch as a 
column, with la larger amounts of eccentricity or crookedness than y were inherent 
in the test specimens, actually fails at loads often considerably below those 
_ given by the empirical formula. After many years of extensive experimental 
research the Society’s Special Committee on Steel Column Research® recom- — 
‘mended a formula of this rational type—the secant formula—as the basis for | 
design specifications. It found this formula not only to be rational and 


definite i in the physical significance of its. ; mathematical terms, but also estab- 


tirely tests. formula is » to any any gives s of 
many eccentricity of imperfection of initial ‘shape, factor of and 


Engineering Association (A.R.E. A. ) are presently 
prac- ‘this formula and, taking advantage of its versatility, are adapted not only to . 
atical @ “standard t structural steel but also to a number of special steels used in bridges. 
>, the r* The author’ Ss proposed formulas a1 are > exclusively | of the first, purely empirical 
n was type and, therefore, in the writer’s opinion, represent a definite step backward 
tation in the direction of formal expressions merely adapted to a given set of tests 
ysical - devoid of fundamental physical meaning . Therefore they are not ap-_ 
along plicable to any structures except those w hich duplicate | the test setup in all 
ength significant details. They are purely mathematical expressions representing 


fixity S-shaped curves with zero ) slope at the lower end and zero asymptote for large 
pes ol wee of Wt; that i is, they ‘merely have t the formal mathematical characteristics: 

of all column without being capable of physical interpretation. 
= 4) In the writer’s opinion the desirable development i in this field lies’ in the 
which direction of rational column formulas of the type mentioned. For standard. 
_— design conditions the way a column curve was originally | determined may be 
of little consequence. Under nonstandard circumstances—such as columns 
of unusual alloys, special loading conditions: (eccentricities), : , and particular 
vale fabricating methods—which n may it involve either larger ¢ ors smaller imperfections 
son of shape than a are commonly met. with, formulas of the rational type ¢ give con- = 
atl ‘Siderable | guidance to the designer even w ithout, or with very limited, test 
sth evidence. The authoi’s formulas, to be applicable, woul el 
jernes The au 1or’s formu as, to be app icable, would require e abor: ate 


tests j in each of these special cases. It is precisely for this reason (to take care 
sed by 
of unusual situations) that the A.R.E.A., in an appendix to its specifications, — 
. =. the complete secant formula y with the physical significance of all its 7 
y term terms explicitly stated. From his own design experience with ‘structures of 
nd the -honstandard material and nonstandard but at least ‘approximately known 
resent ‘Imperfections, the writer can testify to the reliability of this method which, 


minat- ‘in some cases he dealt: v with, was very ‘succesfully verified by a small number 


‘The secant formula, and others of the same type, permit the engineer to. 


"predict the maximum fiber stress under a given load for a column of any 
: 


ve Ol 
5 
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| 
D1 (1936), 


INTER ON COLUMN FORMULAS» 


material provided following characteristics are known with reason: 
accuracy: Properties of of material, end conditions, eccentricities of loading, and | 

‘initial deviations from : straightness (crookedness). 4 By ‘equating this ms aximum 
fiber stress to some stress considered as dangerous (such as the yield point for 


= and by applying ; a suitable factor of safety, a practical column formuls, 
is immediately derived which covers the entire range of slenderness ratios. . The 
_ on that, despite these advantages, the explicit secant formula has not 
become too popular with designers lies in the fact that it is s rather ‘unwieldly. 
‘mathematically. — It requires the use of i trigonometric tables and eannot be 
‘solved e: explicitly for the w vorking stress | for a given _ slenderness r ratio; - that i is, it 
necessitates in each case the drawing of a complete column curve. . This prac. 


P~ difficulty is easily overcome by developing general column formulas of 


the standard familiar type (Rankine and Johnson formulas), which are made 
to fit very closely the values of the exact secant formula. 
The general form of the secant formula is 


(P 


in which (supplementing 1 the notation of the paper) y =w otking stress ‘ 
i i limiti that is, d l f t fibe stress, . 
in compression; oy = limiting (tha ‘is, angerous) value of extreme fiber stress 


4 generally the yield point of the material; 7 = factor of safety; anda = “degree 


- in which ¢ = eccentricity of loading; d = maximum deviation from straight-_ 
= (crookedness) of column; and c = distance from center of gravity i 

outer fiber of cross section, 

ai - Accepted values for these constants are: » = 1.7 to 1 .8 for bridges and 1.5 : 


; to 1.6 for buildings; a“ 0.25 for standard manufacturing conditions and 


- centered connections; and k=0. 75 for riveted and welded connections and P 
A very close approximation of the values of the secant in: the 

and medium range of 1/i is obtained by the expression | 


Eq. 24 is of the Johnson: parabolic type familiar to designing engineers. Most 
column specifications, such as those of the American Institute of Steel 
struction (A.LS.C.) and the American in Railw ay Engineering “Association: 

“(AR R.E.A.), are written in of the Johnson formula for the low and 
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‘wala range of //i. The values obtained from Eq. 2: 24 are extremely close to — 
‘of ithe secant formula, Eq 22, for all values 1s of the working stress down to 


and Eq. 22, 
it for = 0.5 ( -) = 


mula 


in which is the stress at Yi - =0. The of i 


is, it 


For values of L/i larger: than those given by Eq. 27, the Johnson curve deviates 
increasingly from the secant curve and, therefore, should no longer be used for 


| design purposes. — It will be found that for structurally useful alloys the value 
_ _ ranges from about 80 to 160, depending on mechanical properties. 


. as the end point for the Johnson formulas was first sug- — 


(The use of 0.5 


gested to the w aviter by Prof. J. H. Cissel, M. Am. Soc. C. E., of the University : 
of | Michigan, Ann Arbor, _in a discussion concerning general questions of column 


stress 


stress. 


slenderness ratios s larger than |=} — avery close of 
sa ‘the secant dune mula ‘is obtained by t the expression 


‘aight-— 
ity 
s and 

as and J 


4 500,000 


9 


Rg. 28 is adjusted tog give stresses extremely ¢ close to those of the 
(24) ‘secant curve up to values of l/i = 400 and even higher. This. was 
the value ), rom Ba. 28 8 equal to (Euler 
+f Sage Vie = 500. At such large slendernesses the values of the secant formula ar are 
‘ss almost identical with those of the Euler equation which results in the close 
‘‘pproximation of Eq. 28 with the secant at 
ciation Iti is seen that Eqs. 24 and 28 are given in the same general terms as is the 
‘rant formula (o, , E, n, k, a) and in column formulas of the type 


EO 

i 
= 
made 7 | 
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familiar to the designer, that i is, . the Johnson formula for small and medium 
slenderness ratios and the Rankine formula for large slenderness ratios. Thus, 
_ for a given ‘material and given conditions, design formulas, Eqs. 24 and 28, 
: are computed on once and for all by , substituting : appropri: ite values for the con- 
3 stants. _ The resulting formulas are then used in the same simple manner as 
any standard column formula in current codes. In fact, they are of the same 
4 form as those used in in the A.I.8.C. specifications. In | this: manner the use of 
the arithmetically cumbersome s secant formula is entirely eliminated, whereas 


all the advantages are retained. ‘The | formulas are adaptable to any type of 


material of known mechanical properties, which is likely to become advan- 

tageous v with the increasing use of various ferrous and nonferrous alloys for 
structural purposes. The approximation is extremely close for all values of 
@ from 0.1 to about 0.6; that is, in the practically. important range. mo 
“a The use of these leita and the excellent approximation they furnish to 

_ the values of of the secant formula are ill illustrated 1 in the 1e following t two examples: 
Example 1.—Consider a column of structural grade steel, under the fol 
lowing design conditions: oy WA = 33,000 lb per sq in. ;E = 29,500,000 Ib per sq 
in.; 7 = 1.76; a = 0.25; and k = 0. 875. These ‘data are precisely the same 
as those lor pin-ended 1 members in the A.R.E.A. Specifications for 
Steel Railway Bridges, 1941 (see Appendix A). Substitution of the foregoing 
values i in Eqs. 24 to 29 gives the following working formulas: For 1/i = 0 to 


P\ 


15,000 — 0.323 


or 1/i larger than 152— 


is almost identical with that of the A.R.E. A. Bridge Specificati 
(Section III) fo for slenderness ratios great as 140: 


P \ 1/ l 2 
(7). = 15,00 -5 (iy. 


For a slenderness 1 ratio the A. A.R. E. A. specifications suggest the direct 


use of the secant formula, «Ati is thought | that the use of the more convenient 
‘Eq. 30b, instead of the : secant formula, is preferable in a practical way ae j 
~The exact | secant ¢ curve ‘(solid | curve), as Ww ell as the approximating w orking 


“curves” (dotted curves), is shown in Fig. 1(a). The close approximation 


Example 2. —Consider a column of high-strength, -corrosion-resisting steel 


- — the following design gn conditions: Oy = 65,000 lb per sq in.; E= = 29,500, 000 


lb per sq in.; 7 = 1.55; a = 0.35; and k = 0.75. § Substitution of these values 
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(a) Structural Steel, in Railway 
E=29.5x10°, Ty =33000, 1.76, - 
a=0.25 and Ke= al 
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29 gives ‘the following working formulas: For 
venient 31 — ( 0. 


t curves for high- strer 


Seeant formula | are in 10). 


a 

iim 

0 to 

(32a) 

ngtoEqs.32aand32bandtothe 


Discussions 


Examples 1 and 2 2 demonstrate the ‘the writer’ 8 proposil, 
mple form of the resulting working formulas, and the good degree of 
_ approximation as compar ed with the values of the secant formula. — . 
- Lb For each of his nine formulas the author has indicated how they ean be 
made to fit a particular series of tests. The same can be accomplished by the 
secant formula, a as was actually « done in the Society’ s column investigation? 
_ Whereas fitting an empirical f formula to a series of tests consists in the adjust- 
7 ‘ee of arbitrary coefficients, the rational w ay to coordinate the : secant formula 
with test results is to determine experimentally the degree of imperfection : 
entering the formula - To establish the magnitude of this quantity it is only 
“necessary to know e+ -d; that is, eccentricity plus crookedness (compare Eq. 23). 
‘The effect of these two factors is known 1 to be almost identical; that is, an 
a eccentricity. of given amount has the same effect as a deviation from straight- 
ness of the same amount. 4 For this reason it is superfluous to determine e and 
ad separately, since only their sum m enters the secant formula (and Eqs. 24 and 
28). quantity (e+ d) is easily determined in standar ways” by the 
_well- known ‘Southwell method. By this method it is even possible to deter- 
mine the degree of napatteiiian without loading the element to destruction, 
= which i is of considerable advantage i in the case of complicated special structur es % 
Once this s quantity and the physical characteristics of the material are known, 
the secant formula and the hee simplified expressions, Eqs. 24 and 
‘It is quite possible that one of the s several formulas can be made 


ot to fit the secant equation equally as well as as Eqs. (24 and 28. The writer has 
eo ‘these particular formulas because they are of a form familiar to de 
a signers and widely used in design codes. — Experience has shown that there is | 
‘ no particular disadvantage i in using two different formulas for different ranges 
l/i. The main ‘point th the writer w ishes to convey is that, in his opinion, 
ae analytically rigorous” formulas which lend themselves to complete phy sical 


-‘_ aes interpretation are vastly preferable to empirical expressions of a purely formal, 


a A discussion by the author ¢ of the full significance of his arbitrary constant 7 
referred to as the secant yield would contribute to a 


neer in what work, it as 

by the author in the ‘Synopsis’ and in the ‘ “Summary,” from postulates 
ich, from the view ‘point of such an. engineer, are untrue or ‘unimportant. 7 

moe. The‘ ‘Synopsis”’ states that a prime » requisite of any column formula is is that 


“thes average stress ‘ “approach the Euler value as the ratio of slenderness becomes . 
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ifnite.” ” i s the Euler value that ratio i Is zero, and as the column formulas 
in practical use do ‘approach zero as the r: ation, of slenderness becomes infinite a 
statement seems of little importance. However, it is here implied, and 
elsewhere stated, that a a desirable property y of a column formula i is that it shall 
_ be susceptible to extrapolation. _ The practicing engineer is not concerned with — 
_ extrapolation | (as to values of slenderness ratio), as his entire work is w ithin 
s range of slenderness ratios defined by specification. 7 7 
eal ~The ‘ ‘Synopsis”’ further states that a desirable requisite 1 is that “the first - 
Foose of the average stress with respect to the ratio of slenderness be 
continuous.’ This condition certainly does not present itself as a necessity 
the practicing engineer. Why is it stated as axiomatic? 
It is stated in the that “there i is no point in using a formula that 
is s applicable only * chute * up to a given value of the ratio of slenderness * sali — 
AS a matter of fact, almost all of the ‘critical: work done by the structural a 
designer i is in a range up | to a given v: value of the ratio ¢ of ' slenderness, an and a for- 
ula that is applicable w within that r range is ‘exactly what is needed. 
Iti iss stated in the ‘ “Summary” 


“There is to to be said a on the grounds of its sinaplicity, 
5 ja any formula can be plotted to as large a scale as desired and values 
read from the graph, or it may be tabetaten wi w - increments in the sauna 


Asa matter of fact, a practicing engineer will be confronted frequently with a _ 


‘need for an estimate of column capacity, at a time and place where neither 

‘ graphs nor tables are a available; the writer is certain that he has made such — 
_ estimates, in the course of his experience, thousands of times, by the simple 

: application of mental arithmetic, either at a saving of effort as compared with 
looking up a graph or a table, or when neither was at hand. It is very largely 
for “this reason that the ‘prosticing engineer is so so fond of ‘ “straight- line” ¢ 


“ 


ae The attitude of 1 the structural engineer toward column formulas has aa 
ia been oversimplified; but the correct approach for him to take is 
entirely different from the path marked out in this paper. pee a 
The theory ‘leading to the ‘ ‘secant formula” is is as valid a piece of mechanical — 
reasoning as that leading to the customary beam theory, the “theorem of three z 
moments,’ and other too tools of the structural engineer. ‘urthermore, the he secant 
formula has been shown to conform well to tests, w here the tests have been - 
controlled with respect to the variables entering into the formula. | The point 
which the structural ‘engineer 1 must remember, how wever, i is ‘that, in order to 
Bice the secant formula, there must be established not only the properties of the 
: “material and the radius of gyration of the column but also an effective length» 
(between points of contraflexure) and a constant representing the possible 
external moment applied through frame action, eccentricity of end ort 
and crookedness resulting from fabrication, , shipment, or érection. — If practical : 
‘column formulas need to | be rendered more precise ¢ and, therefore, more com- 


plicated, the need i is to reec 
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over the entire e column fi eld. They may be safely estimated for one of 


is to attempt what has been tried in the 
estimate them o once a and for all over the entire field. 
“is ‘It would appear ar to the writer that the most fruitful field of i inquiry now 
open to the structural engineer is to determine, for each different type of column 
(compression members in framed trusses, porncnnty in mill buildings, columns i in 
story buildings, derrick masts, etc.), the ‘most: reasonable allowances for free 
length and for applied moments that should be inserted in the general secant, 
formula, to s afely and economically meet the probable loading” conditions i in 
that particular category of columns. 
‘Having established these several allowances for moment, there can 
a down, fc n or steel of any one yield- -point strength, : a secant formula 
ae expressing the average unit stress at that load which causes a yield- -point stress 
. on the e most severely stressed fiber. _ For a series of other steels, having other 
‘yield points, all these formulas ean again be written, and the difference of 
yield point will make them all different, in detail only, from the corresponding © 
members of the first series, 
yP- It is reasonable to believe, from the work already reported by a committee 
oft the Society, that each of these formulas would quite reasonably represent the 
maximum loads that would be obtained from a series ; of tests, introducing in 
~ each particular case the material and the moments postulated i in set setting up tl al 
ae, The structural engineer would then affect each of these formulas by a factor 
of safety, representing the number of times (usually between 1.50 and 2), by, 
which the anticipated total load could be assumed to be safely multiplied before 
a the extreme fi fiber stress would reach the yield point i in question. _ These then 
- be a set of safe-load formulas, all based ‘upon the same reasonable 
premises, and each safe for practical use. — _ However, one serious practical 
would the secant formula does not lend itself to 


formulas, easier to handle, and not unduly hon series. 


In nthe ¢ case 2 of members i in trusses, is common 


= 15, 000 | 25 


= 17,000 — 0.485 


—for columns i in buildings, up toa a slenderness ratio of 120. Beyond these 


ratios it has been found, these specification-writing bodies, 


— 
ths 
— 
an 
api 
> 
jus 
3 
4 
for 
fee 
by 
— 
| he 
— 
| 
j 
— 
4 
Ls 
—— 
— 


February, 1945 «JONES. ‘ON COLUMN FORMULAS 


oon than those just stated—A.R. E.A. choosing the secant formula itself, 
and A.1.8.C. choosing the Rankine-Gordon type, as a reasonable equiv alent. 
In actual bridge a and building practice, inasmuch as columns of greater — 


= ratios than those to which the second- degree parabola can be _ 


applied either are extremely rare, or are relegated by the specifications to non- _ " 
critical elements of the structure, it would seem to the writer to | be fully 
; justifiable to adopt, for simplicity, for the second portion or end of each curve, _ 
a curve of the Euler form—namely, ‘some easily remembered constant 
in the numerator and the familiar in the denominator (ine aut 


“_ or an inverted parabola, which also can be kept: simple in form an 


\t/ 
oe writer would not wish it to be. inferred ‘that he considers any of the 
- foregoing remarks to be applicable outside the field of usual columns of struc- — 
- tural steel, in bridges, buildings, and like » structures; ir in that field, however, he 
feels that the discussion of column formulas ‘should be entirely uninfluenced 


by the hypotheses and developments « of the paper. er. 
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STRESSES: IN THE LININGS OF SHIELD- -DRIVEN 


TUNNELS 


Discussion 


M. A. SIGVALD 


J OHANNESSON 


D. ‘Assoc. M. Am. Soc. C. E. 7a_The primary ry value 


of this paper is that it presents : a rational method of determining the soil Te 
actions for shield-driven tunnels and ‘dispels much | of the “fog” ‘hitherto : sure 
rounding an analytical treatment of the subject. The high cost of modern, 
_ shield- driven, vehicular and rapid transit tunnels makes it of prime importance 
to attain the utmost economy of design consistent with reasonably | ‘assumed 
loading conditions. | The author’s proposal for reducing the weight of metal 
= the linings of such tunnels deserves the careful thought } of tunnel builders. 4 
An elastic t tunnel ring subjected to a system of active loads will | deform at 
points in a definite and, consistent with these deformations, : soil 
reactions of definite ‘magnitude and direction will appear in the surrounding 


material, dependent only on the elastic properties of the soil and the flexibility 


the tunnel . Th he customary practice of assuming that the 


Tradial or tangent lal forces, yie e momen rust, and shear at that joint, 
adial tial Id th h 
—-: the system is in static » equilibrium. _ When referred to the developed 
- periphery of the ring, the A-constants ‘represent the ordinates to the influence | 
= 
for moment, thrust, and shear at the invert joint due to a unit load 
acting radially or tangentially at sixteen points on the circumference of 3 
For example, the constants in ur and Am, when so plotted, yield the in- 


- fluence diagrams shown i in Fig. 17. To locate the maxima and minima of the 
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influence curves, first derivatives of the expressions in Eqs. 8 are equated to 
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Joints (See Fig. 
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zero, yielding the angular values of these points. _ The influence ordinates — 
corresponding to these values m may then be found by resubstitution in Eqs. 8. 


In T: able: 1 the A-constants have been determined for intervals of 22° 30". 

‘Where force concentrations are applied anywhere between the sixteen tabu- 


+ 
_ 


6 7 8 8 0 ll 
Points (See rig. 4) 
17.—Momunt INFLUENCE Line For ForcEs AND 


lated points, a linear interpolation of the A-values might lead to inaccurate 
results. 7 Influence diagrams | plotted to a large > scale give a truer indication of 
point -to-point variation, permit 1 the intermediate values to be determined more *_ i | 
accurately, and furnish a means for determining graphically the reactions for = 
any system 0 of loads i in ‘static equilibrium. . Since any joint may be assumed as : m5 
4 point of origin similar to the invert joint, the same influence diagrams may ee 
be used for any joint by | shifting the joint and force notation—in @ manner ; 
‘similar to the analytical procedure described by the author i in Section 11. or 
8 are basic for any ring-shaped structure, greatly simplifying the 
mathematical work of deriving ‘equations for moment , thrust, and shear in 


tunnel rings due to : any given loading, as an will s serve to illustrate." is 


forces P concentrated at 6 and 2 An. equation must be 


“Stress Coefficients for Langs Horizontal Pipes,” by Jones M. Paris, Engineering News-Record, 
November 10, 1921, p. 768. 
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th Fig. 18, ee horizontal forces P have been replaced by their radial and 
is, P sin 6 and P cos 0, respectively. Substituting in [J s0i 


iw 


Ea. 6c for. Ae and A. mt from Eqs. 8 and replacing a by and — 0, respec- 
tively, the separate effects of ‘radial and tangential forces may be w written as 
follows: Moment produced by radial forces—_ 


sin 6 (sin @ + 1 — cos 6) — sin8 — 8) sin 


cos 6] “[- — 6) sin? + sind — — sin cos $6). (30a) 
and moments produced by 1 ang force 


= +22 cose x 0 (1 — 8) — — 6) (1 — e086) 


Pr 
- — 6) cos 0 


Adding the s separate effects: 


Pr [— (w — 6) (sin? + 6) + — 8) cos 


sin 0 — sin cos 6] = + 7-6 — sin 16) (1 — cos 

fc The writer would like to call attention to the effect caused by a change in 
the value of the soil compression constant K. For this purpose the moments — 
Sn produced in the tunnel lining s shown in 

Fig. for successive values of K 


from 3 to k= 15 were computed 


| 


| 


| 


The results are summarised in Table 15, 
and moment diagrams corresponding to  — 
_ various values 0 of K have been plotted # 
against the developed ring circumference 
Fig. 19. Values for K = 12 12 in Fig. 19 
Fig. shows the progressive r redue- — 
tion in ‘crown moment with an increase 
in the value of K. if this moment is ; plotted against the ‘soil compression 
F constant as in Fig. 20, it is seen that the variation i in this reduction i is nonlinear. 
Mr. Drucker, although | treating the problem of soil reactions from a different 1 Te 


Assuming, as a basis, a ‘sandy soil a K-value of 1 12 to which the 
= given by the author applies, it is seen from Table 15 that, for de 
‘creases in K to 9, 6, and 3, the crown : moment i increases 9. 5%, 23.7%, and 48.5%) 
espectively. "Similarly, the increase in crown tensile fiber stress is ‘12. 


_. §**Determination of Lateral Passive Soil Pressure and Its Effect on Tunnel Stresses, ds by M. oa rs stay 
Drucker, Journal of tl Institute, May, 1943, p. 499. 
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32.3%, and 65% for the of K. When it is considered that a a 
- having a K-value as low as 3 might be somewhat plastic, and would exert | 
‘TABLE 16. —EFFECTS OF THE VARIATION OF THE ‘Som _ COMPRESSION: 


ON THE TUNNEL LinING 


=" 


“Som, Companenon Constants, K 


Properties (I = 136 In.‘; S = 23.4 In.4; anp A = 27.6 Se. In.) 


(b) Som Reactions, Kips, AND MomeENts, IN Incu-Kips 

-14.0 
6.071 
| 5.164 


. 


4 Joint G 
Joint H 


Joint J 


(c) SerrLemENr 3, IN — 


o> | 766/234 = —3.27|64.6/23.4 = —2.76|56.8/23.4 = —2.43 512/234 = 46.9/23.4 = —2.00 
-|87.4/27.6 = —1.36|37.5/27.6 = —1.36| 37.6/27.6 = —1.36|37.7/27.6 = —1.37|37.8/27.6 = —1.37 


— 


86.4/23.4 = = +3. 69 72.0/23.4 = +3.08)63.7/23.4 = +2.72|58.2/23. 4= +2. 49 54. 2/23.4 = 32 

34.7/27.6 = —1.26/35.1/27.6 = —1.27/35.3/27.6 = —1.28]35.5/27.6 = —1.29|35.6/27.6 = —1.29 

‘ee oo 19.1/27.6 = +0.69/ 19. 1/27.6 = +0.69} 19.1/27.6 = +0.69]19.1/27.6 = +-0.69) 19.1/27.6 = +0. 69 


‘Total | 

for K = 3 are perhaps too high. _ Unless conclusive evidence is at hand as to © 
the n nature of the soil, it would appear prudent to assume a conservative value — 

From the foregoing, it is e lent that i increasing o— soil compression con- 


- Mant ha has the same effect as increasing the flexibility of the lining. _ The simple 
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Haunch 
1a. 19. —Momenrs IN THE TUNNEL LintINnG RELATED TO THE Sor ComPresston Constant 136 N. 
Denot reaction ‘segm 1 to 
te Denoting by Fa the soil reaction on a segment extending from point I 
point 16, Eq. “15 takes, the fo form K or ‘Substituting 
therefore, for in Eq. 32 f 
(Fa cos in = 
— 


relationship existing between I 

— p 
Eq. 18 may be written in the f 
— D1) — 2(P Di) = + => cosa 6 — (82a) 
— 

9 


— 
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«K (Fa cos a, — 


Fs for any given loading, 
: From Eq. 33b, iti is seen that, if the value of Ki is 1 malted ~~ a eran 


ya Mounwr IN THE TO THE Som Compression ConsTANT 


“Satin in other r words, if t the soil compression constant is is ‘increased, it has s the 
same effect as a corresponding increase in the flexibility of the lining. 
Ae : _Itis interesting to note in Table 15 that, for the various values of K for 
~ which analyses were 1 made, the reactive forces F vary only slightly at similar 
points on the r ring for values of K. These small variations, however, 
are sufficient to cause large changes in the moments. = ane a 
«bi is further noted in Table 15 that, as the soil c compression constant is 
decreased, the effects of soil reaction diminish on the upper parts of the tunnel. 


For low values of a“ and a moment of inertia 1 of 136 in in. 4 the reactive. forces are 


ae Contrary to. the usual conception ofa flattening at the invert of the 1 ring ~ 


W with a predominance of downward loads, the results obtained using the author’ Be’. 


; ‘method of f analysis indicate that iti is quite possible for the ¢ curvature to sharpen — ; 
at the invert. ‘The © computed radial deformations of the lining due to com- 
bined active loads soil reactions, corresponding to the assumed values of 

K, have been listed in Table 16. The final displacements are due to three 
factors— “active loads, Soil reactions, and the Actual deformation 
formation means that the point has moved inw and a positive value indi- 
cates an outward movement with respect to the position of the ring before i it 
% is loaded. Table 16 shows s quite clearly: that, for high values of K, points 1 and 


move inw ard with ith respect to joint A. the soil compression constant is 


| 
| 
| 6G 
* 
— 
— 
al | 
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decreased, the soil, —s more ‘yielding, allows the lining to flatten at | Cor 

; —_ near the invert under tt the effects of active loads, until a value of K is 
reached where all points from 1 to. 6, inclusive, move outward. This | is 3 the 
case for an assumed value of | K= 3. a Thus, for high values of K itis possible 


~ to obtain moments at the invert « causing tension on the outer fibers. 7 

TABLE 16. —DEFoRMATION, 1n Incuzs, SHow1NG THE Errect ora Variation 

THE Sort Compression CONSTANT ON THE ‘Rapin 


DerorMATION or THE Lininc (I= 136 ff 


| | | | Kem | 


79.001 1 0033 7 ‘0045 ana 
4 
| 4 to: 0297 “acti 
0 
0331 


70.0014 0.0020 | —0.0022 


With an increase in the compression constant, the soil. 

a the moments in the lower half of the lining, and, with a simultaneous | 
reduction in the amount of the displacement at the ho horizontal diameter as: 

‘shown in Table 16, they c: cause this part of the e ring” to act somewhat as an 
abutment for the upper half. The upper part of of the ring then behaves much 
a as an arch subject mainly to the effects of active loads. 4 ra 
For purposes 0 of comparison, the v writer has analyzed the parr" chosen 1 by -_“lse 
author as an example i in Section 21 (for which I = 136 in.*) by the | to 

tomary method of analysis described in Section 3, assuming (1) the vertical 

4 25 soil reactions to be distributed uniformly over the horizontal diameter (Fig. 1, 
; loading 12); and (2) an induced lateral passive pressure « of an intensity depend- 


increase in the horizontal diameter active loads. The effects 


whi 


Caused by active forces. —230.2| ~173.8 —247.9 sur 

by passive forces. 2 +1 .8| —134.8] —202.5| —134. +134.6| +179.2 pre 

Net moment.......... — 51.0] — 31.9 . .3] — 39.2) — 68.7 

Author’s method (net moment) + 5 3 — 4.5 5 a é — 31.7) — 58.2 | 

of lateral passive pressure are in accordance with the pressure diagram shown str 


_ in Fig. 1, loading 13a, and are based on th the theory developed by Mr. Drucker, a 

using a “le of K = 12. . The results are shown i in Table 17 and the corre-- 

sponding moment diagram has been plotted in Fig. 19. The moments 

de determined by the author are also given in Table 17 for purposes of comparison. - 


| 
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_ Comparing the two sets of values, the dillnenes in . results is apparent at once, 
- especially in the lower parts of the lining. _ The customary analysis yields a 
crown moment 18% greater than that obtained by the author. The invert 


"moment is about 74% of the ‘crown moment obtained by the customary 
whereas the author’ 8 enalyis shows a small moment this is point 


in the lining — from those ¢ by a 


‘saalysi for tunnels is presented in in this paper including the effects of both the — 


‘active loading and the passive soil resistance, or, as the author designates 
the latter, “soil reactions.” ’ It should be of interest to ascertain t the effects — 
of the active loading separately i in order to study the effect, on the final stresses, 7 
of considering the passive pressures. — This would also show to what extent 
reliance is placed on the capability of the passive pressures to keep the stresses 
in the tunnel from exceeding the maximum permissible stress. In the author’s 


however, the two sets of stresses are e inseparably There- 


J 


= 


a2 


to what. appear to be several In pins the | effect of the 
weight of the tunnel ring, 220 lb per sq ft is used. | ‘This weight consists of 
10 lb for the weight of the cast-iron ring having a sectional area of 27.6 sq 
: . and o f 110 lb for the weight of the concrete lining; but, as the concrete _ 
lining is poured i in normal air, a considerable time after the cast-iron ring is 

erected, the v weight of the cant iron alone should have been considered in con- 
- junction \ with the effect of the internal air pressure used to keep out the water — 
- while erecting the cast-iron lining. — In obtaining the moment of inertia and 
section modulus of the tunnel ring, the author correctly disregards the concrete. 


° Another point that may be noted is that the cates apparently neglects the | 
buoyant effect of the water on the soil when obtaining the lateral active pres- 
sures. against the tunnel. ‘Using Eq. 22c he takes a unit active horizontal 
pressure equal to one third the weight of overlying soil—100 lb per cu ft. ‘This ; - a 
corresponds to an active pressure due to a dry soil having : an angle of. — 
of 30°. if the buoyant effect. of the water were considered, the unit. active — 

lateral earth pressure, of course, would be less than 33.3 lb. The variation in ee? ‘ 


stress resulting from this difference in pressure will be considered subsequently. Pi 


‘Tunnel design analysis, for active loading, based on the ordinary method e 
“Tesulted i in a simplification by employing coefficients for bending moments and - 
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DRUCKER ON TUNNEL STRESSES Discussions 
thrusts fc for. the | component parts of the total loading and pressures. x These 
coefficients, for the stresses at point J, are given in Table 18, Cols. 2 and 4. ‘Tes 
In the problem considered by the author (see Section 21), the outside diameter ba 
3 the tunnel is 17 ft 5 in. and the radius at the neutral axis is 8.5 ft. L oaded wit 
tunnel are an earth cover of 36 ft and 2 7 ft of w water. The 


‘TABLE oF STRESSES FOR J (SEE 
‘ON For Benning Momenr AND — 


= 


‘wi = 36 X 100 =3,600| 
ws =27 X25 = 675 0.25 war? 
wm =25 — 0.027 wi 0. 021 wer? | — 38 
wr = 0.009 — 0.021 wi r? 
ws = 3, 600 = = 1,200 0.25 wsr? —+1.00wsr + 10,200 
= 33.3 +.0.208 we r3 0.625 wor? + 
: wo = 27 X 63 = 4, 700 +0.25 wor? 1.00 wor + 14,400 
— 0.208 wn Swurt | + 2,850 
Summation | 20, 570, | + 28,362 
> 


a «The loading numbers cocrenpen te to: ene shown in Fig. 1 1; stealing 6i is omitted because its ef effect i is ine luded 


7 
‘sectional area of the ring is 27.6 sq in., its gross moment of inertia is 136 ins, 
and the intrados section modulus i is 23. 4 in.’, all per foot of tunnel length. 


Applying the coefficients of Table 18 to this case, the stress at ‘point J is 

= 10.53 — 1.03 = 9.50 kips per sq in. | The author finds 


‘avian in. due to the air r pressure, ¥ Ww hich is independent: # the effect a the 
external loading. Therefore, the Temaining stress, without air pressure, would 
be +1. 2 kips per sq in., 80 that, in accordance with the method presented i in 
ee. - the pap paper, the tensile stress reduction due to the p passive ve pressures may be said — 
be 9.5 — 1.2 or 8.3 kips per sq in. 
On the loading diagram, Fig. 1, the author r represents by loading 13a the — 
2S di stribution of passive pressures as considered by the writer. He states (see 
Sect ction y ‘that these | pressures are “ continuous functions of the depth.” ” It 
would | be more correct to describe this m method as one in which the horizontal. :o 
passive pressures are assumed proportional to the lateral deftections caused by 
the active loading. In conjunction with this method, the excess of the down- 
war ard vertical loading over the upward water pressure is assumed to be un 
-forn mly distributed below the horizontal diameter in accordance with general ig 
practice. — _ The author states several times that such assumed vertical pressure 


distribution is ‘the cause “of high calculated stresses. From ‘the following 
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= 


_ To determine the stress reduction due to the passive pressures in accordance 
‘with the writer’s “method, the horizontal deflection at the horizontal ae 
axis due to all the active loading and pressures will first be found from the 
“resulting bending moment at point J—the assumption being that the deflections. 
due to all the active e loadings will be the same as the deflection due to a uniform 
vertical loading ¢ causing the same moment at this point. 


For a uniform vertical loading of w Ib per sq ft, My = 0. 25 wrt. ‘Table 8 

gives as ft-lb. T herefore, 0. 5 wr = 20,570 and w = 


For such uniform load t the e deflection dy at the horizontal | 
“diameter i is 


For the problem under consideration, dix 12,000,000 


“The effect of the passive pressures may 
te in the following manner: p = 1.28 
1.28 x 140 = 1,460 lb . In this com-— 

is the maximum intensity of a 
‘pressure, varying as sin’ as (see Fig. 91), 
‘that. would” be capable of deflecting the 
tunnel back t to ‘its original position. Tt his 
is an extreme pressure and is obtained only 
as a necessary intermediate st step. 

W ith the foregoing information and the : 
given value of 12,000 lb for the soil con 
stant, K, the reduction i in stress, oz, at point 
J , due to. the passive pressures may be found © 
from the following formula derived by the a 
“writer: = 

12 X 0.182 K rp di 0.59Krpd, 


Water Level Below 


= 2:18 X 12,000 x (8.5)? 1,460 x 0.526 
23.4 (12, 000 X 0.526 + 1,460) 
‘ 
1.59 X 12,000 X 8.5 X 1,460 X 


= 8.22 kips per sq in. sieht 
27.6 (12,000 X 0.526 + 1,460 
8  peovtounty | indicated, the author’s method einen a stress at point J of 8 8. 3 i 
- ips per sq in. less than that obtained by the ordinary method for active loading — Re 


only, value \ W vith the foregoing stress reduction of 8. 22 
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ee - the one methods is found to be only about 1%. It is gratifying, indeed, to 
find that the writer’ 8 s method, w hich includes + approximations for the : sake of 
simplicity and w which re requires only - simple slide-rule computations, gives results 
such close agreement w with the author’ more extensive analysis. 


a WwW ith the reduction i in stress of 8.22 kips per. sq in., due to the developed 
a rs per ‘sq. in. w without considering the in due to o air pressure, and 1. 28 + 0.69 
= von = 1.97 kips per sq in. . considering the stress due to air pressure, as against the 


author’ 1.2 kips per sq in. and 1 kips sq ing q 


K, necessary to keep tl the stress fr rom ain ‘the 1 maximum allowable vi wt 
than to begin with an assumed value for this constant to find the resulting 
stress. In the present stage of know ledge of soil mechanics, it would be rather 
‘difficult to defend any ny definite K-value even if it were based on the best possible. , 
available information as to the nature of the soil through which the tunnel is" 
to be constructed. _ Surely, « one could not contend for : a value of K = 12 while 
another argued for a value of, say, 10 or 8. 
Knowing the stress at point J due to the active loading, computed as indi- 
cated in T able 18, and the maximum permissible meine stress, the necessary 


in n which, as s before, oz is the nemeny stress reduction to be provided | by the ‘ 


_ For a maximum allowable tensile stress of 4 kips pe per sq in. for cast iron : and 
taking into account the stress due to air pressure, oz, for the problem under dis 
cussion, would have to be > 9.5 50 + 0. 69 — 4.00 = 6.19 kips per sq_ in. - § Sub- 


this and other known values in Eq. 36, 

26 18 (8.5)? 9.59 x 8.5 6,190 


27.6 1,460 


. +‘Teason 1 for urging a soil constant of K = 12in chetking as to ‘whether the tunnel 
ef “oa is of sufficient strength for a certain location, , when 0 one of only K = 44. 
_ would keep the stress from exceeding 4 4 kips per sq in, 3 ie i 
In this discussion, reference has been made to the stresses at point PE only, 
ee: because it is the point of greatest tensile stress and, therefore, governs fo for cast 

‘Theoretical computations, including those by the author, as well as 
tests on rings used for the London tubes,?° show this to be the case. The 


= oading during these tests (which were conducted under conditions approxi- , 


_ ®**Tunnel Linings with Special Reference to a New Form of Reinforced Concrete Lining,” by G. L. 
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mating those affecting a tunnel after construction) ¥ was ‘increased until cracks 
aa in the. rings—at the | crown of the tunnel. ~ Also, if the stress is known 
at one ne point, the stress at any other point e: can | be determined : readily, if for any 
unusual reason it should be desired to do s 80, , by considering the bending and 
thrust at the point of known stress it in conjunction w with the known or assumed — 
loading, including the passive pressures. 


As previously noted, the author assumed the weight of the tunnel ring as 
220 lb per sq ft whereas he should have used only 1 110 lb for the ring having a + oan 
sectional area of 27.6 sq in. He also neglected to . consider the buoyant effect a 
of the water on the surrounding soil. Making the correction for the weight of | 
- the tunnel ring, the final stress at J for this ring was found to be 1.71 kips per 


in. as indicated in Table 19 for condition Correcting for the effect: of 


— 


‘TABLE 19.—CoMPARISON OF AT Pont 


Author’s specified conditions (see Section 21) and method 

Conditions A, with passive pressures based on the writer’s method 

Conditions B, with corrected tunnel weights> 

Conditions CG, taking buoyant effect of water into account, with angle of 
repose = 30° 

Conditions D, with angle of repose = 45 

Conditions C, except that all active side ‘cudk pressure is disregarded 


‘ 


_ ¢Stresses are in kips per square inch, for section areas of 27.6 sq in. and 11.0 sq in., res ctively. »% For 
_A = 27.6 sq in., the weight was ‘corrected to 110 lb; and, for A = 11.0 sq in., the weight was | corrected 


buoyancy, assuming 40% vi roy and taking the angie of Fepose of the soil as 


= 62.2 per eu on ft and the unit pressure as 62. 2 1 — sin 30° 


20.7 lb. For this s reduced active e earth pressure below 
“water tes stress at J was found to be 2.71 kips per sq in. for condition D. _* 
ee The angle of repose was taken to be 30° as assumed by the author. It may 
not be out of place, however, to give this assumption a little thought. An 
- angle of repose of 30° is considered a reasonable basis for obtaining side pres-_ 
sures on structures, whose safety requires that the highest, ‘Possible active 
_ Pressure be assumed to ) provide sufficient resistance against it. In the case of 
tunnels, however, | as the active side earth pressures reduce the tensile stress _ 
at the g governing point, or at the crown of the tunnel, it becomes 1 necessary to 
decide on an angle of repose that would give : active earth pressures which may — 
be counted on to exist at all times. ba This may be a rather difficult question 
; decide; but the angle of repose for maximum active earth pressure should not 
be used when the minimum pressure | is desired. In this connection, it may be 


noted that Lieutenant- Commander Spangler? assumes no active side earth 


_— ‘The Structural Design of Flexible Pipe Culverts,” by M. G. Spangler, Bulletin No. 153, Eng. Experi- 
j ment Station, lowa State te College, A Ames, 1941, a 
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pressure whatever although he ay analyzes fills at culverts. : herefore, for such : a 
soil as that to which the author assigns a value of of K ines. 2 for the soil constant, | 


which is compact and t undisturbed, an angle of repose of 45° w vould not be too 


large. For such angle of repose, and taking buoyancy into account, the unit 

active earth pressure would be 62.2 X 10.7 Ib. Based on this” 

active earth h pressure the stress. at J was found to 3. 87 ips per in. for: the 


active side earth pressure, oy w was found to be +5. 10 kips per sq. in, r~ con- 
dition the angle of repose of the undisturbed soil would very likely 
- greater than 45°, the final stress, « even for so large a soil constant as K =1 2 
would | be betw een + 3.87 and + 5.10 kips per sq in. 
ra ‘The low stresses determined by Mr. Bull prompted him to state (see Sec- 
27) that w ith a lining as actually designed more than 70% of the 
cast iron is wasted.” ‘That. these low. “stresses a are not inherent in Mr. Bull’s s 
method of analysis is s demonstrated in, Table 19 which shows that, for the same 
tunnel ring, the same low stresses were obtained by the writer’s method under — : 
the conditions assumed by the author. The low stresses. were due to two 
factors: (1) The author neglected the buoyant effect of the water on the soil 
when computing the active horizontal soil pressures; ‘and (2) he assumed a 
high value for the soil constant without assuming correspondingly large angle 
of ‘more correct basis is that represented by condition | E which 
includes the effects of buoyancy and assumes an angle of repose of of 45° Si te 


this condition oy Was found to be 3. 87 gee per sq in. ¥ erga is close shee for 


4 Stress, in Kips per Square Inch 


- To compare the writer’s method with that. of the author for lighter ring» 
sections, a ring having a moment of inertia of 8. 5 in.‘, or one sixteenth of that 
far eonsidered, , was analyzed. « condition B, Table 19, the stress at J 
was found to be 3.70 kips per sq in. as oe 3.18 Kips per sq i in. given by ‘the 


represents a in stress of 62. 3 per sq in. due the passive 
sures. For the other conditions, after making: corrections for the weight of 
the ring, the: stresses increase somewhat more than for the heavier section | 


For the tunnel ring having an area of 27.6 sq in. the s stress: at J, computed 


by the wi riter’ s method, was found to be 0.08 kip higher than that ; given by 
the author. =, For the ring h having an area of 11.0 sq in. the stress obtained d by 
the writer is about 0. 5 kip, higher ‘than that given by the author. © ‘However, 


a me for the area of 17.4 s sq in. the writer ’s value for the stress. was found to be 
0.9 kip less than that given by the author. Ina an to these 


in Fig. be by smooth. curve Lp Ww the 
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tina s stress points, due) to ne stress of 3.14 kips for the area of 17.4 sq in., 
-eannot be so joined except for the stress points for 11.0 sq in. and less, through > 

which a curve Ly could be passed. The irregularity of the other points may 
be explained by the fact that, as .s indicated in Table 14, for areas of 11.0 sq in. 7 
and less the stresses were obtained on the basis of having si six ring segments: 
> in contact with the e soil Ww hereas “J heavier sections sigh ave segments on each 


Kips 


+ 

a 

le 

or 

a "areas between 11.0 and 17.4 sq in. Between these areas there should be « one 


“s for Ww which, if six segments be assumed in contact with the soil, the | pressure on 


a - the sixth segment would prove to be zero or nearly so. . For such sectional” 


‘a area, which may occur about halfway between those for which the stresses are - 
a — & given, and which may be represented by line A- A in Fig. 22, five segments in- 
i. contact with the soil could then be assumed. It would seem reasonable that 
on by extending | line Ly toward the right the stress for this area would be at 
~_ point a or not far fre from it. Similarly, by ahning Te line F toward the left, the 
“ stress for this area of ring would be represented by point b or b; if the | stress 
- line be assumed to curve upward a: as ‘suggested - by the other stress. curves. 
2 Although the computed stress values for the a area at section A-A 1 may be 


7 ‘somewhat different from those indicated in Fig. 22, it would appear that for 
m the same sectional area and loading a consideration of six segments and then 


be &g five segments would give stresses differing by more than a kip for an average 
yo & oftwo stresses of about 3.0 ‘Kips. — - A break i in stress values such as at line A-A 


wd "suggests that the accuracy | of computed stress values, by the author’s method, a 

er, for areas on either side of this line would also be questionable. — If this be so it 
be would indicate that an accurate method of tunnel design, including the effects — 

- it should also be noted that the stresses represented by lines Lg, Lp, a and: Fo 

‘based o on the s loading D shows the stresses that 
n 
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‘Table 19, for the emmanuel a active soil pressures which, in his opinion, 


ee ‘Mr. Bull asserts that the — critical moments obtained by methods other 


‘the horizontal diameter * * *.” this is not the cause of 


P may be noted om the close agreement betw een 1 the stresses for mee same con- © 


tribution the diameter) and those obtained by the author. 
fact, from the vertical components of reactions in Col. ll, Table 11, the p prewvere 
distribution obtained on the horizontal projections of the corresponding seg- 
‘ments is greatest: at the vertical a axis of the tunnel and less toward the sides. 
Such pressure e distribution | causes greater moments at J than does uniform 
upward | pressure for the same total load. | Therefore, the relatively small 
stresses, obtained he either t the author’ s method or the writer’s ‘method, are 
due to the horizontal passive pressures involved in both methods. 7 | 
When a tunnel section is considered or finally decided upon for a project, 
it it may be necessary to to obtain the stresses for various conditions of head of 
water, height of earth : cover, and n nature of soil, for which different soil constant 
values may have to be assigned. To further er simplify the process of obtaining 
the stress at point J for various conditions, : a chart similar to Fig. 23 w wall: 
be useful. This chart was prepared for the ring having a ‘sectional a area of 


sq in. on the basis of the author’s loading assumptions, without 
corrections. 


The notation is defined as fallow oA = bending. stress for active loading, 
including side earth pre pressure, cover = = Ho (see | Fig. 21); = bending stress 


the anther” s assumption), of water = = bending s stress 33 due to all 
active loading (od = 7 TB); Ga = compressive stress due to the thrust 
a caused by active earth ] pressure only; oy = compressive stress in addition t0 oa, 
due to the thrust caused by water, head of water = Hj; and o, = compressive 
stress ina addition to oa, due to the ‘thrust caused by water, head of water = a. 
— ‘Theo oR- -lines, ened on Eq. 36, : are used to obtain the fn al bending | stresses 
_ after taking account of the bending stresses due to the passive pressures. 
" 2 corresponding to the indicated soil constants. _ The other lines in Fig. 23 w ere 
‘ obtained by selecting the proper ‘moment and ‘thrust coefficients from Table 18. 
a From Fig. 23 the final stress can be obtained readily after assuming a 4 


a 


5 Problem 1.—Under the conditions stipulated by Mr. Bull in Section 21 (see 
conditions B, Table 19), let Ha (Fig. 2 21) = 36 ft and H, = (27 ft. Entering 
tao ‘Fig. 23 at H= 36 ft, the vertical dropped to the g,-line in Fig. 23(a) intersects 
: ame at Ta =- 0. A kip per sq in.; and, extended vertically t to the oa-line in Fig. 
intersects at o4 = = 23.7 Kine per sq in. Again: entering ‘Fig. at 
27 ft, the vertical dropped to the in Fig. 23(a) intersects 


‘2 


>= — 0.6 per sq in.; and, extended vertically to the os-line i in Fig. 23(b), 
at og = — 13. 2 kips per sq in. The bending stress op due to all 

; 


arg 


i 


Stresses, in Kips per Sq In 


» nl Stresses, in Kips per Sq In C4 Values, or Tp Values for Use in the Tr Curves) 


a 


i 
— 
ue 4 
— 
q 
— 
— 
— 
— 
~ 
| 
= 


Q 
= 
a 
=) 
Zz 


DRUCKER 


4 


‘February, 1: 


ae 
pue jo sanjer 


oF 
o a 


Une 


a 


HS Hg, Hy OR He (Fig. 21) 


! 


“uy sed sdiy ul 


4 


a | 


y | 


ay) ui asp 10) sanjea 


40 — 


2. 


| 


CAC 


Te 


ry 


TRESS AT 


| bg sad 4g 'sassasys 
jad sdiy ui‘ 4 
Ss 


| 


‘NING THE 


ETERMI 


Soil Constant, K, 


had 


Jomnt J F 


CoMPREHENSIVE CuarT For D: 


‘Tunnet Ring Havina an OF 27.6 Sq IN. 


4 40 DEPTS Me, OR H, (Figg21) = | —— 
| 
a) > 
4 eae AY 4 
se 
it 
— 


JOHANNESSON (ON TUNNEL STRE STRESSES 
active loading is 23.7 13.2 = 10.5 kips per in. In Fig. 23(6), the inter- 
= "section: of a horizontal from op = 10.5 and a vertical from K = 12, gives” 
= 2.3 per sq in 1. The? total of og = — 04 6 =-1 0 kip. 
po sq in., and. os, W without considering air pressure, is 2 2.3 — 1.0 = 1.3 kips per 

ie = Adding the a air-pressure effect, + 0. 69 kip per sq in., the final value of _ 

Cos = 1340. 69 = 1.99 kips per sq in., which checks the stress given in 
Table e 19. - For a total stress oy limited to 4 0 kips per sq in., or could be as 
reat as 4.0 + 1.0 -0. 69 = 4.31 kips" per sq_ in. Entering Fig. 23(b) at 


= 10. 5 kips per sq in. and drawing a line horizontally to or = 4.31 yields 


2 (vertically) K = = 4. 4, which demonstrates how the nec necessary soil constant m: ay. 
be obtained so » that the stress shall not exceed a certain value. 


Stevatp JoHANNeEsson,” M. Am. Soc. C. E.*—The use of cast iron for 
tunnel linings was conceiv ed by Mare Isambard Brunel i in 1818, but it was 

not. until the Tower Subway under the Thames River in London, England, 
_ was built i in 1869 that a tunnel was constructed w with cast iron as lining material. 


Since then cast-iron tunnel | lining has been used extensively for : shield-driven 
tunnels, particularly i in Great Britain and 


method for determining the stresses in linings, the author is be 
ue ‘congratulated for having made a most important step in this development. a 
9 Based on his formulas the author finds that an excessive volume of metal — 
is used in the present cast-iron nm linings. Although this may be true 1 theoret- 
"ically, it is still questionable whether cast-iron linings could be materially 

decreased in weight. _ The author gives two reasons why this may not be done: 


The first is the difficulty of casting large segments of slender shape; and the 
a second is the need for a large cross section to resist the pressure of the shoving 
jacks. Also, cast- -iron linings are cast in segments ts bolted together at the 


Bee, horizontal flanges, | the deflections of w hich, under circumferential tension, con-— 
found the t theoretical conclusion; : and forces « originating in the interior of the 


tunnels, such as an explosion « or a , collision, might cause the tunnel to ‘collapse. 
if the > cast-iron li lining were as thin as the theoretical stress computations w ould 
‘indicate to be safe. - These conditions have been, and must be, considered i in 
proportioning the lining as long as the present type of cast-iron lining is used. 
_ However, with the e development of the art of making, shaping, and fabri- 
¢ cating steel, the: present type of "cast-iron. lining should , and will, become 
—. _ obsolete, and it will b be superseded | by linings of steel (as ‘suggested by the 
gt writer in 1922)”, provided that the natural hesitancy of tunnel designers to 
break with precedent may be overcome. One of the greatest obstacles to 
a has been the uncertain knowledge of the stresses to which a tunnel 
lining may be subjected. ea he author’s contribution to the determination ball 
= stresses will greatly facilitate the solution of this problem. | 


Chf., Highway Planning Bureau, State Highway Dept., Trenton, N. J. 7 
Received by the Secretary January 23, 1945. 
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ESTON GAVETT, M. Am. So _E. 
profiles, obtained from careful observ vations of velocities i ina wide open channel 

; with known bottom roughness, 2 are an important part of this « comprehensive 


‘a The author states (see heading, , “Results: Velocity Distribution”) that “In 
call r runs the velocity distribution fit fits the logarithmic law very well, although 
‘near the bottom there i is a systematic deviation from the law, which is not 
understood. if Figs. 6 and 7 are plotted with semilogarithmic coordinates in 
accordance with the von K4rmén velocity defect (Eq. 14), which the 
author states (see heading, “Previous Work”) «* * -* has been shown by — 
Nikuradse®- to apply to pipes by G. ‘Keulegan” { to apply also to 
open channels” 
_ Theg general | acceptance of the von ] Karman law which the author endorses 
‘in this and in an di discourage interest in the 


{ 


| 
> 
= Note. —This paper by Vito a -Vanoni, Assoc. M. Am. Soc. C. E., was as published i in n June, 1944, Pro- 
Codings. Discussion on this paper has appeared in Proceedings, as follows: December, Ralph Ww. 


Received by the Secretary January My 1945, 


‘“Gesetzmissigkeiten der turbulenten Strémung i in Glatten by Nikuradse, 


1“Strémungsgesetze in Rauhen Rohren,” by J. Nikuradse, Forschungsheft 361, Vereines Deutscher 


“Laws of Turbulent Flow in Open Channels,”’ by Garbis H. Keulegan, J Journal of Research, ‘National 
Bureau of Standards, U. S. Dept. of Commerce, Vol. 21, 1938, pp. 707-741. _ 


Velocity Distribution in Channels,” by Vito A. Vanoni, Civil Engineering, June, 1941, PP. 
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ON SUSPENDED SEDIMENT 


iden ed the Hissar law and that n varied with the Reynolds amex. 
Prandtl‘ has mentioned the seventh- -power law based on Blasius’ experimental 
law of pressure drop and has. _ shown that, although the seventh-power law 
breaks d down for Reyn nolds’ numbers above about 50,000, similar. derivations: 
may be made for other po pow wers. The writer has plotted a number of velocity 
profiles f from various sources on log-log paper and has found that the greater 
number show good agreement: with the exponential formula. =” 
_ The velocity defect law has the advantages of being cndvenesl and of locating - 
the approximate velocity profile when the friction velocity i is known, _— 
The exponential formula also has merits: @ In. giving a more accur: ute 
expression for velocity distribution and (2) : in providing | a simple e expression — 
for the velocity distribution i in a pipe or channel and perhaps also for that i ina 
plate or side of a . ship. ; To apply the exponential formula, it is necessary to 


know the appropriate value of n, which appears to vary in accordance with 


o 


“Tae 


at 


cal 
definite laws that should be determinable from sufficient experimental data. ia 
For smooth pipe, varies with the Reyr nolds number and decreases as wl 
increases. it also varies in smooth pipe with 4, the laminar ‘thickness, showing §§ ind 
an exponential ‘relationship. If nis plotted against 6, the laminar thickness If 
for smooth pipe, and ¢, “the equivalent | sand ‘size for Tough pipe, on log-log cor 
paper, the (n-e)-curves for rough pipe appear to intercept the (n- 5)-curve for = 
smooth pipe at points determined by a value of r/e and to branch off with @ 
slopes greater th than the h pipe curve 
The ‘author states” (see heading, “Results: ‘Velocity Distribution”) that 
“The fact that the velocity profiles fit the logarithmic law so — means that in 
the average ‘velocity will always occur at the depth -= 0. 368. 


If the exponential law is valid, there will be some variation with 


yal 

> replotting the points fen Figs. 6 and 7 on | log- log paper, the w writer te: 

- diana a mean value of n of 0. 145 for runs 1A, 2B, and 3C with a bottom a 

i ‘ roughness 0 of 0.47-mm sand and no suspended load. For bottom sand size of - 


0.88 mm and no suspended load (runs 14A, 14B, and 21C), the value of mwas Bh lop 
0.175. hain these values of n, the mean velocity depths from Eq. 30d : are as | 


Description 
Roughness coefficient n. 0.145 
Relative depth ratio ——=.... 0.393 


Ymax 


3 

““Gesetzmassigkeiten der turbulenten Strémung in Glatten Rohren,” by J. Nituradee, Forachungshat Am 
wd 356, Vereines Deutscher Ingenieure, Vol. 3, 1932, pp. 1- 36.0 “iy a 

48 ‘Applied Hydro- and Aeromechanics,” based on lectures of L. Prandtl, by ( O. Tietjens, 
Hill Book Co., Inc., New York, N. Y., 1934, p. 70. 
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The variation in YavelYmax i is not great, but, for the extreme values of n for 
‘smooth and rough surfaces of about 0.10 and 0. 36, the corresponding \ values of 
would be 0.385 and 0425.0 

_ For pipe, the author®® has g given a constant mean distance of 0.223 ro. The 
following einnomnat are derived from the exponential formula for pipe: 


| ‘For values of n of 0. 1 and 0. 2, the mgeniien' values of Yavg would be 
on 0237 rand 0.2697, 


1a ie The author’s data indicate that the bottom sand size, equivalent to that 
to used | by Nikuradse, is smaller than the actual size of sand used, Possibly 
cause a part of the sand was embedded in the coating. 
ta. bs The average velocities computed by Eq. 30a check closely with the chi 
Rk velocities given. by the author computed from the velocity defect law, 
ing § indicating that the von K4rm4n law may be correct for the average velocity. a 


ess If so, the formula of B. A  Bakhmeteff, An. Soe. C. E., is and 
log. convenient: 


hat in which is value = 0. 4 ‘in Eq. 32, 


-; X — : same of 3. 75 results from substi-_ 


0a) | tuting 0.223 for y yliro an and k = 0.4 in Eq. 14, Prandtl’ s value for this constant 
se B according to Professor Bakhmeteff is 4.07, which agrees better with = values 


J. Wirzia,® Assoc. M. Am. Soc. C. E. service 

ter to: the science of hydraulics | has been rendered by Mr. Vanoni i in enlarging and — 
om completing: his } previous paper on sediment transportation,®® and he is to be 
e of commended for a , careful expos ition and presentation of ‘details leading to. re. j 
was logical conclusions. problem of sediment suspension and transportation 
‘ mii by water involves fundamental principles that are only imperfectly understood, - =4 
4B, and only by painstaking ‘Tesearch will limited existing know ledge of the ‘phe- ae 2 
na Mechanics of Turbulent Flow,” by Boris A. Bakhmeteff, Princeton Univ. Princeton, 
2d Ed., 1941, Eq.116,p.70, | 


4 “Engr. (Civ.), U. S. Engr. Office, Buffalo, N. Y. 
al Received by the Secretart January 15, 1945. 
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me Experiments on the Transportation of Suspended Load, Vito A /Vanoni, Transactions, 
— Am. Geophysical Union, Pt. III, 1941, pp. 608-621, y — 
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Discussions 


it is that the distinction een these modes is 
because of the insurmountable difficulties of treating the sediment load of a a 
stream as a whole, there should be no slighting of the possibility that a ‘much. 
= larger quantity of | sediment in some stre: ums may move as bed load than as 


_ The relations and the results: giv en in this paper should not be applied 
‘indiscriminately to the total ‘suspended load. The exponential vertical 
tribution of ‘sediment applies only to the coarse suspended | load. How ever, 
in some 1e streams, a large part of the total suspended load consists of 1 very fine 
silt, termed the e “wash load,’ *? and i is very nearly uniformly distributed over 
the stream cross section. | “Any conclusions based on a logarithmic or expo- 
nential distribution of the wash load would be misleading or erroneous. — — Appa- 


Tently an appreciable wash load did not ¢ occur in the author’ s experiments, in 


which a . carefully - graded sand was used. 7 It would be of interest, however, to 
know Mr. views on the probable i of distributed 


was 
load, on the pediment transfer 


Although probably, of merely academic interest, it seems. appropri: ite to 

po point out that a transitional mode of sediment transportation exists, known as_ 
“saltation. ” Tn saltation, grains of bed material are lifted momentarily into. 


: the m moving fluid, then bouncing back on to the bed, | possibly dislodging other 
particles When the turbulence of exceeds a certain 


Cc. E.. , concluded that, when velocity is enough cause s ae 
tion, , the turbulence level is also sufficient to obscure the saltation effect en- 


tirely, by placing the material from stream bed into 
paper succinctly 


Ln 


is equal to the momentum transfer Em, and that these 
are constant over the depth. The experimental results described demonstrate” 
= falsity of the former assumption, but it is not apparent whether they either 
prove or disprove the latter. The writer considers it probable that s some of 
the apparent discrepancies i in the published data, such as the large var iation 
in the values of z and a, run 17, Table 3, might | be due to some | unaccountable 
de tion in the ( 
eviation in the y)-relation, 

_ The author explains | that the laboratory fi flume used in his experiments h had 

rubber er sides, with apparently a very small ‘roughness’ coefficient. A narrow 
flume, how ever, is: bound » to exert a pronounced “side- wall « effect,” ” especially 


ii 


ao 


rs 
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Mechanics for Hydraulic Engineers,’”” by Hunter Rouse, McGraw-Hill Book Co., Inc, 

Discussion by Joe W. Johnson of ‘Sedimentation in Reservoirs,’ "by Berard J. Witsig, : 

Am. Soc. C. E., Vol. 109 (1944), p.1072, ly 


&“Criteria for Determining Sand- Transport by Surface and by A. Ka Jinske, 
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is chown in n Fig. 5. Although Fig. 5 demonstrates the effect, 
experimental flume extended only about one sixth» of the flume width| 
from each side, the w riter is of the opinion that such an extensive series of tests 


uch - should have been made in a wider flume. A fundamental assumption for the 
.as J two-dimensional analysis of the problem i is that the theoretical equations for 

a the vertical distribution « of velocity and of sediment are valid only in “wide” 
lied channels. In this case, t the ratio of depth to width was a about 1:5.5, and its 
dis- + recognized effects on the. experimental observations are indicated by the fact 


- that all measurements in Table 3 were confined either to the center line of the 
or to a vertical not than 0.5 ft from the center line. 


in the of width to ‘depth should desirable i in future experiments. 
The effect of the depth-w width ratio in a turbulent. and rapidly flowing 
natural stream ‘is illustrated by "discharge measurements made in 1943 under 
q ‘the v writer’s direction near the head of the Cascades, in in the American channel 
the Ni viagara River at Niagara Falls, 1 > At the time of the measurements, 
the average depth was about 6 ft; the w idth at the | gaging section, about 480 ft; 
the measured discharge, 16,150 cu ft per sec; and the mean velocity, 7.5 ft per 
sec. Vertical velocity traverses made with a current meter showed — 
1 as ~ generally exponential gradient existed, with the maximum velocities occurring 


r, to 


ided 


into atthe surface. These me sasurements also indicated that the mean velocity in 
ther 2 a vertical occurred, on the average, at the 0.55 depth from the surface, instead 
hese of at the 0.6 depth commonly assumed in ‘discharge measurements. ~The 
tion BD logarithmic law places the mean at the 0.632 depth from the surface. It seems 


Soe. | probable that such variation in field observations from theoretical assumptions 
sand laboratory results must place in doubt the validity of generally applying 
the latter to the large-scale phenomena existing in natural streams. a 

_ The writer w vas particularly interested in the analysis of the effect of : 
plied - sediment loads on the roughness coefficient of the flume. : Runs 18, 20, , and 21, 
"Table 4, how ever, apparently it indicate some inconsistencies. For the clear 


sient, . water run in each case, Manning’s n was 0.0118, 0.0121, ,and 0.0128, respectively; Mi 


> 
— 


ents whereas, with a sediment concentration of 1.2 g per liter, n was 0.0107, 0.0116, 
trate and 0.0116, ‘respectively. . Runs 18 and 21 showed a ratio of 1.10 for a 
‘ither flows: compared | with silt flows, but run 20, under practically identical condi- 
ne of tions, gave a ratio of 1. 04. Can this apparent inconsistency be explained by — Pare 
ation bo change e in slope and v velocity, or or should it be ascribed to the difficulties of 

a | . It w would be of interest to find some explanation of the mechanism w hereby ; 
aaa ; the channel roughness is decreased by the presence of sediment in the fluid. 


arrow - May not the lamina of high concentration fill the rough projections of the — 

cially channel bed, thereby acting as a lubricating film and decreasing the apparent 

roughness coefficient? — Merely as a speculation, it is also suggested that an 

increase in ‘silt load —— the apparent Toughness at a diminishing ra rate >i 
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N SUSPENDED SEDIMEA 

~The author suggests that one of the reasons for the ‘difference ‘ides 
the momentum and sediment transfer coefficients is the “slip” between the 
fluid and the ‘sediment particles. _ The writer considers that this ‘ ‘slip”” is due 
in part to the inertia a and the “drag” of the particles. — ~The difference in the 
effect of fine and coarse sediments on €s “compared with Em might reasonably 
be expected. Small particles, with a relatively larger surface area. & per unit 
ae * mass as compared with larger particles, should naturally be more readily 
- moved or kept in suspension _ by the random turbulence; the “slip” effect, also. 

should be less_ because o of the smaller mass mass and drag; and | molecular impact 


should be more pronounced. 


; — It is suggested that discrepancies between z nzand 21, as W well as the ‘consistent 
» variation actually demonstrated, might be caused by secondary, or more 
properly “tertiary,” ” circulation i in the third (horizontal) dimension. mn. If such 
baa. 4 variations appear in carefully controlled laboratory tests, even wider \ variations 
may be expected in natural streams, where com plex circulations are induced 
_by winds at the surface, bends in ‘linement, tributary and surface inflow, and 


irregularities of the sides and bed. 
“The description | of the horizontal distribution of sediment in segregated 
bands and clouds omits reference to any observed vertical segregation. Did 
such vertical separation | exist, and, if so , did the sediment tend to follow the 
‘fluid filaments of highest velocity? A The writer’s observations of silt- laden 
_ natural streams indicate that sediment flows appear to move in “slugs” or 
7 - clouds, but he has not noticed band-like movements. Iti is probable that these 
“may not occur in natural streams because of the large degree | e of turbulence. 
_ With reference to Fig. 15(c),_ it is not readily apparent how al secondary 
circulation i is derived from the velocity and sediment. diagrams. 


Ais whether increase in velocity secompanying silt 


An increase in ‘stream velocity causes increased erosion the 

. ee: stream is not burdened to capacity with its sediment load. When the stream | 
. Boe; 7. is fully charged with sediment, it seems probable that the concentration neat 7 
oe the bottom may be so great that no net increase in load may result, although a rs 
— continual interchange takes place. % Such a stable condition may be due, not so * 
_ much to the theoretical decrease i in shear accompanying the smaller depth and é 
as to the f fact tha at the full capacity is already 
ey. Bs _ previously an anticipated effects of sediment loads on the transfer coefficient, and ; 
‘ eos that his experiments indicate the relative magnitude of these effects for the ’ 
of concentrations and sizes of sediments studied. The apparent effect 
a Be Oe on channel roughness 0 or resistance suggests that further r intensive research on ye 
5 Boy particular aspect of the problem would be of interest economically, as it a 
oe ‘may appreciably affect practical design of r pipe lin lines or flumes for commercial in 
or industrial transportation on of materials i in suspension. — 
The author has shown that a complex problem hydraulics is. 
susceptible, at least i in ‘part, to ‘rational and experimental — Other | 
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EVAPORATIO ON 


5 ad 


Rouwer,! E.4«—Observations on the evaporation 


from free water surfaces in pans are usually made either to discover tk 
mental laws showing the relation between the the meteorological da data and ‘evapora- 7 


when developing evaporation formulas to simplify t the analysis of the effect of 
the influencing factors. — Since the laws governing certain similar physical 
phenomena such a s turbulence, heat conduction, and mass transfer have been’ 
fairly well established, the analogies between t these phenomena and evaporation 
are sometimes used as a basis for attacking the problem of determining | the 
laws of evaporation. This i is the method adopted by the author. However, 
evaporation laws cannot be developed by analogy alone because ‘there are 7 
certain constants involved which must be determined experimentally, 2 and = 
for which carefully controlled evaporation observations are necessary. 
‘The evaporation formula (Eq. 5) based on analogy between heat transfer 
and diffusion developed by the author is represented by the dashed line in : 
Fig. 5. _ Evaporation data obtained | by the author, by means of the optical | rs 
interferometer designed by him , When plotted on the same figure, result i in the | 7 
_ two full lines that meet at the point where the vapor-concentration difference — 
is: 4.2 x 10~ * lb per cu ft. Since the line representing Eq. 5 and the lower — 

- part « of the line representing the experimental data both pass through the origin, 
—-~Eg. 5 could be made to conform to the experimental data by multiplying by a - 
“constant. However, Eq. 5 cannot be made to conform to the remainder of the = 
- experimental data by : a simple transformation. — The fact that there is a break — 
in the line indicates that a different law applies beyond this point. mele ta 


Note.—This paper by G. H. Hickox was published in October, 1944, Proceedings. 
ee ating, C6 Engr., Div. of Irrig., SCS, U. 8. Dept. of Agriculture, Colorado Experiment Station, Fort 
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4 probably ar are never conducting evaporation ¢ even 
are excluded, any ‘difference in temperature or density will cause the air to 
- This i is confirmed by | the experience of the writer in conducting ¢ ev rapora- 
. ae experiments under controlled conditions in the laboratory at Fort Collins 
(16a),*° although the exact significance of the phenomenon was not recognized 
at the time. The evaporation was at a maximum when the water was warmer 
than: the air. When this occurred, the air in contact with the w vater Ww ould be 
warmer than the surrounding air and would rise, “carrying with it a certain 
_ amount of water vapor which Ww ould have had to be removed by diffusion if 
the: air were absolutely still. Ww hen the water was colder than ‘the : air, the 
water surface because it was held there by the rim of the - tank. When this 
‘occurred the water vapor could id escape only by diffusion, which is a slow process. ‘ 
_ evaporation dra the break in the curve (Fig. 5) is about four times that com: 
ima 
puted by Eq. 5. — In fact, the discrepancy beyond the break in the curve is so | 
_ According to the author's analysis the evaporation should vary as the: tha 
0. 75 power of the wind velocity. The results obtained by several investigators. phe 
investigators varies considerably, with one exception the results conform to the Hy 
ndamental 
or whether it is due to the uncertainty of wind measurements at low velocities (se 


> 3 under controlled laboratory | conditions. . Although all external disturbances 
= 
a density gradient was reversed, and the air would remain in contact with the hia 
These . explanations, how ever, -, do not account for the fact that the observed 7 
at that it doubtful whether the analogy of mass transfer is valid. 
are shown plotted i ink Fig. 8. Although the total effect of wind | found by tl these i 
0.75 power at the higher velocities. Ww hether this is a fundamental difference J tion 


is not known. _ The author states that F. Graham Millar was able to approxi- = 
mate actual al evaporation records on the basis of turbulence and eddy diffusivity, 
t but could correlate the writer’s data only o on the assumption n that the ‘anemom- : 


eter was in error. 7 Cup anemometers probably are not the most accurate 


devices measuring wind velocity, but they ‘are almost universally 
and for this reason they were adopted for the experiments at Fort Collins. . 


we 


——— the indicated velocities differ from the true velocities, particularly 
high velocities, similar differences would occur whenever the anemometers. 


"were used and consequently comparable results would be obtained. 
ni ‘The effect of pan diameter on ev yaporation has been plotted by the author 
in | Fig. 9, showing that ; beyor ond 10 ft the diameter has very little effect on the 
vaporation. is the | generally accepted conclusion. However, the 
author’s formula for ev aporation into air in motion, Eq. 156, indicates that the 
- evaporation is inversely proportional to the Seuetli root of the diameter of 
pan. According 1 to this formula, the evaporation from a ft pan ould 
be more than five times as great as that from a pan 9,000 ft in diameter under 
similar conditions instead of about 10% greater as indicated in Fig. 9. 
aie i ig. 10 shows the relation between the barometric pressure and the relative 
evaporation at different altitudes based on data obtained by the writer, 
— in accordance with Eq . 20. - *Fig. 13 shows the same data and in 


addition the relative evaporation as as by the writer (16d). will be 


aA: 4> Numerals in parentheses, thus: (16a), refer to corresponding items in the Bibliography (see Appendix 
of the paper), at th the end of discussion i in issue, 
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observed th: it the points fall much closer to a line wit 
kad 


do those computed by the author. This indicates analysis 


used by the writer gives results seiaian more nearly in accord with the facts 


Relative Evaporatio 


H. Hickox 4 


(b) 12-Fr SuNKEN Pan aT 


1940 | 1941 
0.85 

66 

0.58 

(0.68 

(0.77 

0.72 


Mean 


. 
February 


phenomena. 
tion from a a Weather Bureau pan to that from a reservoir. . The author cites 
| 1939 
0.79 


38 
The Sub-Committee on Evaporation of the Special Committee on Irrigation — 
(see heading, “Introduction’’) that the ratio is probably not constant through- 
TABL 5. —Monruty Ratios FOR THE 
E.sinore, 
0.90 


than the analysis | based on the analogy betwee een evaporation a 
Hydraulics of the Society (2) adopted 0.70 as the Tatio ¢ of the annual evapora- 
& CONVERSION OF THE EVAPORA TION 
WEATHER Bureau Pan To THE EVAPORATION FROM OTHER 
0.69 
0.72 


00 G0 G0 G0 


0.90 

(0.76) 


6. 
0.78 
(0.78) 


Sesss 


(0.78) 


the year. Although this view is proved by the writer’ experiments at 
Fort Collins (16c), it is doubtful whether it is worth while to use & different — 
ratio for each month because conditions vary from ye to year on account o. 


nich 
« 
ven 
| 
this 
rved 
— 
is — 
the: Pa 
hese 
the 
ence 
rity 
— 
im 
ithor 
0.78 | 0.69 | 0.87 0 — 
er of 0.84 | 0.84 | 0.92 | 0 0.71 
-ould O90 0 0.76 
inder Mean (weighted) (0 (0.78) | 
— 
. 
Young, Assoc. M. Am. Soc. C. E. (43)(44), for the Division of Irriga- 


tion, Soil Conservation Service, ‘Uz Ss. Department of Agriculture) shows that 
- the ratio reached a maximum during the cold months i in the fall and winter at 
Elsinore, ‘Calif. whereas at Fullerton, Calif., the maximum ratio occurred 
~ during: the warm months of the year. The 1 mean values for these two ratios 
Z: the entire year are, -Tespectiv ely, 0. 77 and 0.78 instead of 0.70 as recom- 
_ mended by the writer and by the Sub-Committee on Evaporation of the 
_ Special Committee on Irrigation Hydraulics of the Society. _ There was little 
y change in the ratios from year to year. _ These r results § support the author's 
_ contention that the ratios would not be the same under Southern California 
conditions as they would be in areas where the evaporation is very small during 
4 thewinter months, 
a The optical interferometer developed by the author i is an ingenious device 
— ea for measuring very small rates of evaporation with precision. It is unfortunate 
¢ a = the experiments with this device were limited to a narrow range of condi- 


tions. A series of observations under controlled conditions i in the laboratory 
on the | evaporation from pans of different sizes into air moving at definite 


we : velocities would | have made e it possible for 1 Mr. Hickox to check the accuracy 
? of Eq. 15b w ithout making assumptions as to the conditions under which the 


Pe * study of the device raises the questions of Ww vhether the plate glass re- 
-flector F, in Fig. 3, would not reflect the beam of light from both surfaces of 
the plate glass : and of ae as a result, two sets of interference bands would 
appear to the observer. one from the interior face of the glass plate 
_ would be less distinct, but, when the dark bands of one set J w ere superimposed 
7 over the light | band of the other, the difficulty of counting the bands would 
be increased. This possibility is mentioned | because of a similar difficulty 
Pl -- experienced by the writer in making | evaporation observations by the use of 
ay >. - the optical lever. Until the plate glass mirror first used was replaced by 3 
‘ metal: mirror with only one reflecting surface, two distinct images of the 
| which the evaporation was read appeared in the use 
of a metal mirror completely eliminated the second image. 2 
After studying the paper the reader cannot help but feel that ene i is con- 
siderable doubt as ' to whether the analogy used is suitable for interpreting 
a evaporation n phenomena . At least the available data do not fit the theory very 
WwW well. The writer does not wish to leave the impression that he is attempting to to 
ah ‘Mr. Hickox’s work even though the results are not conclusive. 
Much of the work on evaporation has been on an empirical basis and any 
attempt at a rational analysis of the problem should be 
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Like Prof. O. Mohr,!? Mr. Nubar considers the geometric and 
loading factors as forces applied at the joints of the frame, although, regrettably, 
his illustrative example (Fig. 13) fails to demonstrate any superior a 
in practical application. By graphical methods, the moments i in the waseda 
shown i in Fig. 13, can be determined with less effort and time. ee 
Graphically, the writer would proceed in two steps:!* The first step is to 
determine the physical properties of the frame: by locating the characteristic 
points, termed ‘ ‘central points,” ‘ “inflection points,” and ‘ ‘elastic points.” ” 
The relative lances of the characteristic points is itenainat by constructing 
a “three-line polygon,” as described briefly by the writer elsewhere.'® The 
second step is to determine the moments produced by the given loading and _ 
Referring to Fig. 20, w hich is an . analysis of the frame shown in Fig. 13: 
The central moments i in member CD with a load of 40 kips are M; = M; = 80 


If the joints are the ‘moments in the frame are 


moments due to distortion i in member CD are 


bs Nore.- —This | paper by Yves ogy Am. Soc. oc. C.E,, , was published i in April, 1944, Prodeedings. 


Discussion on this aper appeared in Proceedings, as follows: June, 1944, » by Leon Beskin; and September, — 
16 


Associate Bridge ‘Engr., ” ‘Bridge D Dept., Div. of Highways, State State Dept. of of ‘Public ¢ Works, Sacramento 
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“926 Stresses i in Cecil Vivian vo von A Abo, Transactions, Soc. E., Vol. 89 

of Continuous Premes by Distribution of Moments,” A. 
Mento, Calif., 1943 Ed. Published by the author. 
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With a unit displacement at joint C is shown in Fi 


Mpc = N bh = = (668) 


and ky is a distortion coefficient; | a and ba are the the points of 
inflection (show n in the parentheses, Fig. 20(a)); A is transverse d displacement 
member ends (Fig. K (= = I/L) is the basic stiffness of member, 
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Diagram 
shown in the circles in Fig. 20(a). ~The ‘numerical factor 600 Eq. 67 is, 
introduced for convenience computations. . The basic moments in each 
member determined by Eqs. 66 and 67 and distributed, as shown in Fig. 
j i a a condition of static equilibrium, the reaction at joints E, F, and D i “i ZT 
a ‘Using the moments in Figs. 20(c) a and 20(d) to determine the values a He, ‘Hr, es 
‘and Hpin Eq. 68, a value of k, = 0.48 is computed. £ 
Superposing the moments in Fig. 20(c) and the moments in Fig. 20(d), the 
ag final bending moments in the frame are Moc = 161.5 and Mgr = 29.1, w hich 


ees differ slightly from those computed by the algebraic formulas introduced by 
Contrary to Mr. Nubar’s 3 method, the operations in the graphic distribution 2 
of moments may easily be visualized _ 1 Therefore, the moments may readily 
be checked. ~ Furthermore, graphics may be 5 ened i in checking the analy tical 
‘computations of moments. . The constant C used by the author in Eq. 35a 


may be checked by the ‘inflection point in the graphical ‘method. 


Vie 
Ths, for member cD, Fig. 20(a): Ci 0.22, which is the 


that computed by the author. of the inflection point dis 


"tances, a and b, in members with variable sections be simplified by using 
constants given by the writer elsewhere. 
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